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1.0  INTRODUCTION 


1.1  Significance  of  the  Problem 

Advanced,  high  performance  composite  materials  are  really  material  systems.  The  consti¬ 
tuent  materials  interact  in  such  a  way  that  their  collective  response  is  more  than  the  linear  sum  of 
die  response  of  the  constituents.  This  simple  reality  provides  the  technical  community  with  a 
remarkable  opportunity  to  create  composite  material  systems  which  arc  uniquely  suited  to  per¬ 
form  specific  engineering  tasks.  At  the  same  time,  this  systems  aspect  of  composite  materials  is 
a  very  great  challenge  to  the  research  community.  It  introduces  complexity,  nonlinearity,  and 
scaling  problems  (to  name  a  few)  which  require  new  developments  to  represent  geometry  and 
material  behavior,  from  the  standpoint  of  mechanics. 

This  challenge  is  even  more  formidable  when  one  attempts  to  find  modeling  approaches  to 
the  representation  of  the  long-term  response  of  composite  systems  to  cyclic  mechanical,  chemi¬ 
cal,  and  thermal  loading.  The  systems  aspect  in  that  context  creates  the  need  to  represent  defor¬ 
mation,  degradation,  aging,  and  other  processes.  These  process  are  "multidisciplinary”  in  every 
sense,  and  the  mechanics,  chemistry,  thermodynamics,  and  physics  of  their  activity  is  generally 
coupled. 

The  "performance"  of  a  material  system  is  not  a  material  property  like  stiffness.  It  depends 
on  the  manner  in  which  the  environment  of  mechanical,  chemical,  and  thermal  loads  are  applied, 
and  on  the  history  of  that  application.  The  physical  events  that  determine  performance  are  often 
a  "process,"  with  rates  and  interactions  that  must  be  considered  and  characterized  in  order  to 
properly  describe  and  anticipate  the  consequences. 

Recent  advances  in  the  technologies  of  manufacturing  and  materials  have  enhanced  the 
current  application  of  composite  materials  from  being  used  as  secondary  structural  elements  to 
becoming  primary  load-carrying  structural  components.  Consequently  thicker  and  thicker  com- 
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posites  structures  are  being  made  to  cany  higher  loads.  Due  to  the  inherent  inhomogeneity  and 
anisotropy  of  the  materials,  analysis  of  these  composite  structures  imposes  new  challenges  to 
engineers.  A  widespread  and  efficient  application  of  composite  materials  requires  detailed  and 
reliable  knowledge  of  their  physical  properties  and,  in  turn,  of  their  behavior  under  applied 
loads.  There  are  a  number  of  important  technical  problems  associated  with  the  mechanics  of 
composite  materials.  One  such  problem  is  die  effect  of  discontinuities  (holes  and  notches)  on 
the  strength  of  composite  laminates.  This  issue  is  critical  for  the  determination  of  the  load  bear¬ 
ing  capacity  of  composite  laminates;  which  is  directly  applicable  to  the  design  of  composite 
panels  and  the  location  of  fastener  holes.  Indeed,  the  manufacture  and  repair  of  advanced  com¬ 
posite  structures  have  serious  problems  connected  with  die  placement  of  fastener  holes.  This  is 
especially  relevant  to  composite  panel  repair,  both  in  the  field  and  at  the  repair  facility.  At  the 
present  time  all  depots  are  confronted  with  these  problems.  The  lack  of  appropriate  data  has 
resulted  in  new  and  in-service  designs  which  are  often  unnecessarily  conservative  and  expensive 
(both  in  cost  and  tum-around-time).  Another  related  problem  is  the  issue  of  interlaminar 
response  of  composite  materials  which  is  directly  related  to  delamination  and  edge  effects  in 
composites.  In  recent  time,  delamination  has  become  the  most  feared  failure  mode  in  laminated 
composite  structures.  It  can  exhibit  unstable  crack  growth,  and  while  delamination  failure  itself 
is  not  usually  a  catastrophic  event,  it  can  perpetrate  such  a  condition  due  to  its  weakening 
influence  on  a  component  in  its  resistance  to  subsequent  failure  modes.  Study  of  delamination  is 
one  of  the  prominent  topics  in  composite  mechanics  research.  Another  issue  in  die  engineering 
application  of  composite  materials  is  the  modeling  and  study  of  structures  with  curved 
geometries.  Because  of  the  complex  nature  of  these  structures,  present  computational  capabili¬ 
ties  are  far  behind  the  engineering  developments  and  only  very  limited  simulations  of  these  sys¬ 
tems  are  feasible  presently.  All  the  foregoing  problems  share  one  common  deficiency,  namely, 
the  lack  of  an  adequate  and  sound  theory  predicated  upon  principles  of  continuum  mechanics 
that  could  be  implemented  through  an  accurate  and  efficient  numerical  scheme.  Here  in  Berke¬ 
ley  Applied  Science  and  Engineering  Inc.  (BASE)  we  started  this  research  to  address  these  very 
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basic  issues  associated  with  application  of  composites.  In  particular  die  following  objectives 
were  followed  during  the  course  of  this  research. 

1.  Development  of  a  thermomechanical  theory  for  composite  laminates  that  has  a  continuum 
character.  The  theory  should  be  able  to  account  for  the  three-dimensional  responses  of 
laminated  plates  and  shells.  The  theory  should  also  account  for  effects  of  micro-structure, 
anisotropy,  and  geometric  nonlinearities. 

2.  Formulation  of  the  theory  in  the  context  of  finite  element  and  numerical  implementation  of 
the  theory  through  computationally  efficient  algorithms  suitable  for  composite  applications. 

3.  Verification  of  the  theory  through  analysis  of  a  series  of  benchmark  problems. 

1.2  Theories  of  Laminated  Composite  Plates  and  Shells 

There  has  been  an  increasing  amount  of  research  activity  pertaining  to  the  mechanics  of 
composite  laminates  and  multilayered  plate  and  shell  theories.  The  scope  of  multilayered  shell 
theories  encompasses  all  the  methods  used  in  conjunction  with  two-dimensional  treatments  of 
composite  shells.  These  methods  generally  lead  to  a  system  of  partial  differential  equations  in 
two  independent  spatial  variables,  along  with  a  set  of  boundary/initial  conditions  compatible 
with  them.  As  in  the  case  of  single-layer  isotropic  shells,  all  the  different  approaches  for  con¬ 
structing  multilayered  shell  theories  can  be  viewed  as  either  a  single  approximation  or  succes¬ 
sive  approximations  of  three-dimensional  elasticity  models.  For  a  review  and  a  complete  list  of 
references,  the  reader  is  referred  to  Noor,  A.  K.,  and  Burton,  W.  S.  [1990].  The  following  four 
general  approaches  for  constructing  two-dimensional  theories  for  multilayered  shells  were 
identified  in  this  work: 

1.  method  of  hypothesis; 

2.  method  of  expansion; 

3.  asymptotic  integration  technique; 
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4.  iterative  methods  and  methods  of  successive  corrections. 

The  first  approach  is  an  extension  of  the  Kirchhoff-Love  approach  and  is  based  on  intro¬ 
ducing  a  priori  plausible  kinematic  or  static  assumptions  regarding  the  variation  of  displace¬ 
ments,  strains  and/or  stresses  in  the  thickness  direction.  The  simples  of  these  hypotheses  is  the 
linear  variation  of  the  displacement  components  used  in  conjunction  with  first-order  shear  defor¬ 
mation  theories.  Although  the  method  of  hypotheses  has  die  advantages  of  physical  clarity  and 
simplicity  of  applications,  it  has  the  drawback  of  not  providing  an  estimate  of  the  error  in  the 
response  predictions. 

The  second  approach  was  initiated  by  Cauchy  and  Poisson  around  1828,  and  is  based  on  a 
series  expansion,  in  terms  of  the  thickness  coordinate  for  displacement  and/or  stresses.  For  iso¬ 
tropic  and  anisotropic  plates  and  shells,  power  series,  Legendre  polynomials,  and  trigonometric 
functions  have  been  employed.  The  second  approach  also  includes  the  method  of  initial  func¬ 
tions  in  which  the  displacements  and  stresses  are  expanded  in  a  Taylor  series  in  the  thickness 
coordinate.  The  relations  between  die  higher-order  derivatives  of  each  of  the  displacements  and 
stresses  and  their  lower-order  derivatives  are  obtained  by  successive  differentiation  of  the  three- 
dimensional  elasticity  relations. 

In  the  third  approach,  appropriate  length  scales  are  introduced  in  the  three-dimensional 
elasticity  equations  for  the  different  response  quantities,  followed  by  parametric  (asymptotic) 
expansions  of  these  quantities  in  power  series  in  terms  of  a  small  thickness  parameter.  The 
three-dimensional  elasticity  equations  are  thereby  reduced  to  recursive  sets  of  two-dimensional 
equations,  governing  the  interior  and  edge  zone  responses  of  the  shell.  The  edge  zone  (or  boun¬ 
dary  layer)  is  produced  by  self-equilibrated  (in  the  thickness  direction)  boundary  stresses.  The 
lowest-order  system  of  two-dimensional  equations,  depending  on  the  choice  of  the  length  scales, 
corresponds  to  the  thin-shell  approximation.  The  higher-order  systems  introduce  thickness 
correction  effects  in  a  systematic  and  consistent  maimer.  This  approach  was  first  applied  to  iso¬ 
tropic  shells  by  Reissner  [I960].  Later,  it  was  extended  to  anisotropic  shells. 
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The  fourth  approach  includes  various  iterative  approximations  of  the  three-dimensional 
elasticity  equations,  trd  predictor-corrector  procedures  based  on  a  single  or  successive  correc¬ 
tions  of  the  two-dimensional  equations. 

The  following  comments  cm  the  different  approaches  of  constructing  two-dimensional  shell 
theories  and  the  boundary  conditions  to  be  used  in  conjunction  with  these  theories  are  in  order 

1.  The  state  of  stress  in  the  shell  can  be  decomposed  into  an  internal  state  of  stress  and  a 
boundary  layer.  The  first  is  generated  by  external  surface  forces,  and  by  boundary  and 
reactive  stresses,  which  are  not  self-equilibrated.  The  boundary  layer  is  generated  by  self- 
equilibrated  (in  the  thickness  direction)  boundary  stresses.  The  method  of  hypotheses  and 
the  method  of  expansion  can  describe  well  the  internal  state  of  stress,  but  are  not  suited  for 
describing  the  boundary  layer  (because  of  tire  complicated  nature  of  die  displacement 
and/or  stress  distribution  through  die  thickness).  By  contrast,  the  asymptotic  integration 
technique  is  well  suited  for  describing  both  the  internal  state  of  stress  and  the  boundary 
layer  of  the  shell. 

2.  If  the  method  of  expansion  is  contrasted  with  the  asymptotic  integration  approach,  the  fol¬ 
lowing  two  major  differences  can  be  identified: 

(a)  No  a  priori  assumptions  are  made  regarding  the  relative  magnitudes  of  the  different  stress 
components  in  the  method  of  expansion.  By  contrast,  in  the  asymptotic  integration 
approach,  such  assumptions  have  to  be  made  either  explicidy  or  implicitly. 

(b)  Whereas  the  method  of  expansion  leads  to  a  set  of  simultaneous  equations  in  all  the  param¬ 
eters,  the  asymptotic  integration  technique  leads  to  recursive  sets  of  equations  for  both  the 
interior  and  the  edge  zone  (or  boundary  layer)  of  the  shell.  The  lowest-order  equations  for 
the  interior  of  the  shell  correspond  to  the  classical  Kirchhoff-Love  theory. 
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3.  The  aforementioned  four  approaches  are  not  mutually  exclusive.  Some  of  the  theories 
developed  can  be  classified  in  more  than  one  category.  Also,  hybrid  methods,  combining 
more  than  one  approach,  have  been  proposed.  Examples  of  these  hybrid  methods  are:  (a) 
the  use  of  a  three-dimensional  model  for  the  core  and  a  two-dimensional  model  for  the  fac¬ 
ings  of  a  sandwich  shell,  and  (b)  the  two-step  approach  based  on  using  a  two-dimensional 
shell  theory  to  evaluate  the  in-plane  stresses  and  then  applying  the  three-dimensional 
equilibrium  equations  to  evaluate  the  transverse  shear  and  normal  stresses  in  laminated 
composite  shells. 

4.  Although  most  of  the  theories  developed  for  laminated  composite  shells  replace  the  actual 
shell  (or  each  of  its  layers)  by  a  smeared  ordinary  continuum,  some  microstructural  and 
generalized  continuum  shell  theories  have  been  proposed.  In  the  first  class  (microstructural 
theories),  the  shell  is  considered  to  consist  of  alternating  layers  of  relatively  rigid  material 
(with  properties  representative  of  fibers)  interspersed  between  flexible  layers  with  proper¬ 
ties  typical  of  the  matrix  material. 

5.  The  derivation  of  the  correct  boundary  conditions  for  a  particular  shear  deformation  shell 
theory  from  prescribed  data  seems  to  be  important,  even  for  thin  shells.  Some  recent  work 
on  isotropic  plates  and  shells  indicated  that  the  use  of  approximate  boundary  conditions  in 
conjunction  with  higher-order  shell  theories  can  lead  to  significant  errors  in  the  predictions 
of  the  shell  theory.  Therefore,  Saint  Venant’s  principle  needs  to  be  re-examined  when 
applied,  in  conjunction  with  higher-order  shear  deformation  theories,  to  shell  problems. 

Extensive  research  effort  has  been  devoted  to  the  classical  laminated  theory  (C.L.T.)  in  the 
past  and  a  huge  amount  of  literature  is  available  on  this  topic.  The  classical  laminate  theory  is  a 
direct  extension  of  classical  plate  theory  in  which  the  well  known  Kirchhoff-Love  kinematic 
hypothesis  is  enforced.  This  theory  is  adequate  when  the  thickness  (to  side  or  radius  ratio)  is 
small  and  anisotropy  is  not  pronounced.  The  range  of  applicability  of  the  CL.T.  solution  has 
been  well  established  for  laminated  flat  plates.  It  indicates  that  a  theory  which  accounts  for  the 
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transverse  shear  defamation  effects  would  be  adequate  to  predict  only  the  gross  behavior  of  the 
laminate. 

In  order  to  overcome  die  deficiencies  in  C.L.T.,  refined  laminate  theories  have  been  pro¬ 
posed.  These  are  single  layer  theories  in  which  the  transverse  shear  stresses  are  taken  into 
account  They  provide  improved  global  response  estimates  for  deflections,  vibration  frequencies 
and  buckling  loads  of  moderately  thick  composites  when  compared  to  the  classical  laminate 
theory.  A  Mindlin  type  first-order  transverse  shear  deformation  theory  (SD.T.)  was  first 
developed  by  Whitney  and  Pagano  [1970]  for  multilayered  anisotropic  plates,  and  by  Drag  and 
Tso  [1972]  fra  multilayered  anisotropic  shells.  Both  of  these  approaches  (CJL.T.  and  S.D.T.) 
considered  all  layers  as  one  equivalent  single  anisotropic  layer,  thus  these  approaches  are  inade¬ 
quate  to  model  the  warpage  of  cross-sections,  that  is,  the  distortion  of  die  deformed  normal  due 
to  transverse  shear  stresses.  Furthermore,  the  assumption  of  nondeformable  normal  results  in 
incompatible  shearing  stresses  between  every  two  adjacent  layers.  Also  the  later  approach 
requires  the  introduction  of  an  arbitrary  shear  correction  factor  which  is  dependent  on  the  lami¬ 
nation  parameters  for  obtaining  accurate  results. 

The  exact  analyses  performed  by  Pagano  [1989]  on  the  composite  flat  plates  have  indicated 
that  the  distortion  of  the  deformed  normal  is  dependent  not  only  on  the  laminate  thickness,  but 
also  on  the  orientation  and  the  degree  of  orthotropy  of  the  individual  layers.  Therefore  the 
hypothesis  of  nondeformable  normals,  while  acceptable  for  isotropic  plates  and  shells  is  often 
quite  unacceptable  for  multilayered  anisotropic  plates  and  shells  with  very  large  ratio  of  Young’s 
modulus  to  shear  modulus,  even  if  they  are  relatively  thin.  Thus  a  transverse  shear  deformation 
theory  which  also  accounts  for  distortion  of  the  deformed  normal  is  required  for  accurate  predic¬ 
tion  of  the  behavior  of  multilayered  anisotropic  plates  and  shells. 

Along  this  line  the  work  of  M.  Epstein  and  P.  G.  Glocker  [1977,1979],  P.  M.  Pinsky  and  K. 

O.  Kim  [1986],  and  J.  N.  Reddy  [1988,  1993]  can  be  mentioned  where  the  theory  of  multi- 
director  surfaces  was  used  to  model  multi-layered  plates  and  shells.  Pinsky  and  Kim’s  work  was 
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based  on  multi-layered  shell  theories  of  Epstein  and  Glockner  where  the  concept  of  multi¬ 
director  field  defined  over  one  reference  surface  was  employed  for  the  description  of  the  initial 
geometry  and  motion  of  multi-layered  shells. 

Reddy  proposed  a  displacement  based,  layerwise  shear  deformable,  C°  theory  which  also 
accounts  for  the  warping  of  the  composite  cross  section.  In  his  theory,  there  is  a  single  reference 
surface  and  a  director  is  associated  with  this  reference  surface.  The  variable  kinematic  finite  ele¬ 
ment  is  developed  by  superimposing  several  types  of  assumed  displacement  fields  within  the 
finite  element  domain.  The  underlying  foundation  of  die  displacement  field  is  provided  by  the 
assumed  displacement  field  of  any  desired  equivalent-single-layer  theory  and  the  layerwise  dis¬ 
placement  field  is  included  as  an  incremental  enhancement  to  this  underlying  field.  This  work 
has  been  reported  for  linear  analysis  and  for  fiat  geometry  domain. 

Unlike  the  equivalent  single-layer  theories,  the  layerwise  theories  assume  separate  dis¬ 
placement  field  expansions  within  each  material  layer,  thus  providing  a  kinematically  correct 
representation  of  the  strain  field  in  discrete  layer  laminates  and  allowing  accurate  determination 
of  ply  level  stresses.  During  the  course  of  this  research,  we  developed  a  layerwise  shear  deform¬ 
able,  multi-director  theory  which  directly  address  the  technical  drawbacks  present  in  most  of  the 
theories  that  have  been  proposed  for  composite  analysis  to  date.  The  main  features  of  the  theory 
are  summarized  as  follows: 

•  The  displacement  field  proposed  in  this  work  is  continuous  in  3-D  where  as  the  rota¬ 
tion  field  is  layer-wise  continuous  (in  2-D)  and  can  be  discontinuous  across  the  finite 
element  layers  through  the  thickness  direction. 

•  The  displacement  field  fulfills  a  priori  the  static  and  geometric  continuity  conditions 
between  contiguous  layers. 

•  The  novel  idea  in  the  assumed  displacement  field  lies  in  its  capability  to  model  the 
distortion  of  the  deformed  normal,  without  increasing  the  number  and  order  of  the  par- 
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theory. 

•  Another  new  idea  in  the  theory  is  its  3-D  feature,  thereby  modeling  the  interlaminar 
conditions  and  predicting  the  3-D  edge  effects  more  accurately. 

•  A  salient  feature  of  the  proposed  theory  is  that,  at  most,  only  first  derivatives  of  dis¬ 
placement  and  rotation  fields  appear  in  the  variational  equations.  The  practical 
consequence  of  this  fact  is  that  only  C?  continuity  of  finite  element  functions  is 
required  which  is  readily  satisfied  by  the  family  of  Lagrange  elements. 

•  The  number  of  partial  differential  equations  in  the  resulting  system  is  independent  of 
the  number  of  plies  and  their  orientations  in  the  composite. 

•  Another  advantage  of  the  proposed  composite  shell  theory  lies  in  the  greater  flexibility 
in  the  specification  of  the  boundary  conditions. 

•  The  theory  covers  a  wide  range  in  the  sense  that  in  one  limit  case  when  there  is  only 
one  layer  of  proposed  elements  through  die  thickness,  one  recovers  the  features  of  the 
standard  Shear  Deformation  Theories  (S.D.T.).  However  the  added  advantage  in  the 
present  case  lies  in  the  3-D  feature  of  the  theory  which  controls  the  variation  in  the 
thickness  via  the  Poisson  terms  rather  than  ad  hoc  mathematical  tricks  as  done  in  the 
literature. 

•  In  another  'unit  case,  one  can  model  the  composite  with  one  element  per  ply  through 
the  cr  ite  thickness,  a  procedure  that  is  typically  done  while  using  the  standard 
3-D  anisotropic  elasticity  elements.  The  added  advantage  of  the  proposed  theory  in 
this  limit  case  is  that  because  of  the  shear  deformation  capability  of  the  proposed  ele¬ 
ments,  they  model  the  warping  of  the  deformed  normal  more  accurately,  thereby 
improving  the  bending  behavior. 

•  From  a  practical  design  point  of  view  it  provides  the  engineer  the  freedom  to  deter¬ 
mine  the  precision  in  analysis.  If  a  general  response  of  the  composite  structure  is 
required,  the  composite  can  be  modeled  with  one  element  through  the  thickness.  On 
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the  other  hand,  the  designer  can  model  die  thickness  with  as  many  layers  of  the  pro¬ 
posed  element  as  deemed  necessary  to  achieve  the  required  accuracy. 

•  Furthermore,  it  is  feasible  to  employ  this  formulation  for  constructing  plate  and  shell 
finite  elements  via  the  finite  element  displacement  method. 

Details  of  this  work  are  discussed  in  this  report  and  a  summary  of  the  content  of  die  report  is 
presented  next 

1.3  Summary  of  the  Report 

The  results  of  our  efforts  during  the  course  of  this  research  are  presented  in  section  2 
through  section  14  of  this  report.  A  summary  of  the  contents  of  these  sections  is  presented  in  the 
following. 

In  section  2  the  kinematics  of  the  micro-  and  macro-structures  were  examined  and  the  rela¬ 
tionship  between  strain  measures  at  micro-  and  macro-levels  were  derived.  The  field  equations 
for  composite  laminates  were  derived  through  a  direct  integration  of  field  equations  of  classical 
continuum  mechanics.  The  linearized  kinematic  measures  were  derived  in  the  context  of 
infinitesimal  deformation  and  the  relation  of  linear  strain  measures  with  displacement  vector  and 
director  displacement  vector  were  obtained.  The  equations  of  motion  in  the  linear  theory  were 
derived  and  were  presented  for  both  curved  and  flat  geometries. 

Section  3  showed  the  derivation  of  constitutive  relations  for  composite  laminates.  A  pro¬ 
cedure  for  deriving  the  relation  between  composite  quantities  (i.e.,  composite  stress  tensor  and 
composite  couple  stress  tensor )  and  strain  measures  at  macro-level  were  presented.  The  deriva¬ 
tion  was  performed  for  a  bi-constituent  composite  laminate  and  the  constitutive  relations  were 
expressed  in  terms  of  material  constants  associated  with  every  individual  layer. 
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Section  4  presented  the  complete  theory  for  linear  elastic  composite  laminates.  The  rela¬ 
tionship  between  the  displacement  vector  and  the  director  displacement  vector  was  derived 
based  on  the  geometrical  continuity  at  interfaces.  The  field  equations  were  derived  in  terms  of 
displacement  vector  and  it  was  shown  that  classical  continuum  theory  can  be  derived  from  Cos - 
serat  composite  theory  (or  the  case  of  a  single  constituent  The  theory  was  further  simplified  for 
bi-laminate  micro-structure  composed  of  isotropic  constituents.  Finally  the  constitutive  relations 
for  composite  stress  tensor,  composite  couple  stress  tensor  and  interlaminar  stress  vector  were 
derived  in  terms  of  the  displacement  vector,  its  gradients  and  material  constants  of  individual 
constituents. 

Section  5  was  the  extension  of  the  theory  for  multi-constituent  composites.  The  micro¬ 
structure  or  representative  element  was  assumed  to  be  composed  of  several  constituents  which 
repeated  themselves  in  the  layering  direction.  The  development  of  this  section  is  particularly 
suited  for  fiber  reinforced  composites  where  the  fiber  direction  changes  in  the  stacking  sequence 
of  the  plies.  The  theory  was  simplified  for  the  case  of  isotropic  constituents. 

Section  6  presented  the  extension  of  the  theory  from  a  purely  mechanical  theory  to  a  ther¬ 
momechanical  theory.  In  this  section  composite  field  quantities  corresponding  to  the  heat  flux 
vector,  the  heat  supply  and  the  specific  entropy  of  classical  thermo-mechanical  theory  were 
introduced  and  the  equation  of  local  balance  of  energy  and  the  Clausius-Duhem  inequality  were 
derived  in  terms  of  these  composite  field  quantities. 

Section  7  presented  the  constitutive  relations  of  linear  thermoelasticity  for  composite  lam¬ 
inates.  These  constitutive  relations  were  derived  for  the  composite  stress  tensor,  composite  cou¬ 
ple  stress  tensor,  entropy,  heat  flux  vector  and  heat  flux  couple  vector.  The  developments  of  this 
section  were  parallel  to  those  of  section  4  and  a  set  of  coupled  thermomechanical  field  equations 
in  terms  of  the  displacement  vector  and  the  temperature  were  presented. 
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In  section  8  a  linear  theory  for  cylindrical  laminates  was  presented.  Relative  kinematic 
measures  for  cylindrical  geometries  were  discussed  and  linearized  field  equations  along  with 
constitutive  relations  in  cylindrical  coordinate  systems  were  obtained.  The  theory  was  extended 
to  thermoelasticity  and  explicit  thermoelastic  constitutive  relations  for  isotropic  layers  were 
derived. 

Section  9  followed  developments  parallel  to  section  8  but  for  composite  laminates  with 
spherical  geometry.  The  theory  was  extended  to  thermoelasticity. 

Section  10  presented  the  results  of  stress  analysis  of  a  composite  laminate  with  traction  free 
edges.  The  problem  of  a  finite-width  symmetrically  laminated  composite  plate  under  uniform 
one-dimensional  stretch  was  studied  and  it  was  shown  that  the  present  theory  captures  the  three- 
dimensional  response  of  the  laminate  at  the  free  edge  boundaries.  The  predicted  results  were  in 
agreement  with  earlier  studies  of  the  subject. 

Section  1 1  presented  the  wave  equations  in  laminated  fiat  composites.  Expressions  of  wave 
velocities  for  longitudinal  waves,  horizontally  polarized  shear  waves  and  vertically  polarized 
shear  waves  were  derived. 

Section  12  followed  studies  parallel  to  those  of  section  11  but  for  elastic  waves  in  cylindri¬ 
cal  and  spherical  laminates.  Expressions  for  motion  of  rotary  shear  waves,  axial  shear  waves 
and  radial  waves  in  cylindrical  composite  laminates  were  derived  and  general  solutions  in  terms 
of  Hankel  functions  of  first  and  second  kind  were  presented.  Similar  developments  for  spherical 
laminates  were  followed.  It  was  shown  that  for  isotropic  materials  the  displacement  equations  of 
motion  for  waves  with  polar  symmetry  can  be  recovered. 

In  Section  13  the  finite  element  formulation  of  the  theory  was  presented.  The  approach 
proposed  in  this  work  utilized  a  displacement  field  which  fulfilled  a  priori  the  static  and 
geometric  continuity  conditions  between  contiguous  layers.  The  order  of  the  system  was  the 
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same  as  in  the  first-order  shear  deformation  theory.  The  chief  advantage  of  the  assumed  dis¬ 
placement  field  rested  on  its  capability  to  model  the  distortion  of  the  deformed  normal  and  to 
satisfy  the  continuity  requirements  without  increasing  the  number  and  the  order  of  the  partial 
differential  equations  with  respect  to  the  first  order  transverse  shear  deformation  theory.  The 
theory  was  used  to  construct  plate  and  shell  elements  for  composite  laminates  which  accounted 
for  the  3-D  effects,  through-the-thickness  variations  of  stress  and  strain  measures,  and  permitted 
the  warping  of  the  deformed  normal.  These  capabilities,  in  particular  for  curved  geometries,  are 
unique  features  of  the  present  developments.  Based  on  these  developments  accurate  stress 
analysis  of  composite  shell  structures  is  no  longer  a  formidable  task. 

Section  14  presented  the  results  of  several  finite  element  modelings.  These  analyses  were 
performed  for  both  flat  and  curved  geometries.  Various  fiber  orientations  were  considered  and 
different  loading  conditions  were  examined.  The  study  included: 

•  extension  analysis  of  flat  composite  laminates  with  free  edges. 

•  bending  analysis  of  composite  plates  with  different  boundary  conditions. 

•  stress  analysis  of  composite  laminates  with  geometric  discontinuity  in  the  form  of  a 
circular  hole. 

•  bending  analysis  of  cylindrical  shell  composites  with  free  edge  conditions. 

These  analyses  showed  the  main  features  of  the  present  theory.  In  particular,  an  accurate  model¬ 
ing  of  discontinuities  in  composites  and  analysis  of  laminates  with  curved  geometries  was  shown 
through  the  application  of  the  proposed  theory.  These  unique  enhancements  of  mechanics  of 
composite  materials  provide  the  required  computational  capabilities  for  further  application  of 
composites. 


2.0  MICRO-MACRO  CONTINUUM  MODEL  OF  COMPOSITE  LAMINATES 


2.1  Kinematics  of  Micro-  and  Macro-Structures 

Let  the  points  of  a  region  3^. in  a  three  dimensional  Euclidean  space  be  referred  to  a  fixed 
right-handed  rectangular  Cartesian  coordinate  system  x‘  (i  =  1,2,3)  and  let  6*  (i  ■  1,2,3)  be  a  gen¬ 
eral  cmvecud curvilinear  coordinate  system  defined  by  the  transformation  x‘  =  x‘(&0-  We  assume 
this  transformation  is  nonsingular  in  %,  Furthermore,  let  £  represent  die  coordinate  of  a  micro¬ 
structure  in  the  layering  direction  with  £  *  0  corresponding  to  the  bottom  surface  of  the  micro¬ 
structure.  We  recall  that  a  convected  coordinate  system  is  normally  defined  in  relation  to  a  con¬ 
tinuous  body  and  moves  continuously  with  the  body  throughout  the  motion  of  the  body  from  one 
configuration  to  another. 

Throughout  this  work,  all  Latin  indices  (subscripts  or  superscripts)  take  the  values  1,2,3;  all 
Greek  indices  (subscripts  or  superscripts)  take  the  values  1,2  and  the  usual  summation  conven¬ 
tion  is  employed.  We  will  use  a  comma  for  partial  differentiation  with  respect  to  coordinates  9“ 
and  a  superposed  dot  for  material  time  derivative,  i.e.,  differentiation  with  respect  to  time  hold¬ 
ing  the  material  coordinates  fixed.  Also,  we  use  a  vertical  bar  ( I )  for  covariant  differentiation. 
In  what  follows,  when  there  is  a  possibility  of  confusion,  quantities  which  represent  the  same 
physical/geometrical  concepts  will  be  denoted  by  the  same  symbol  but  with  an  added  asterisk  (*) 
for  classical  three  dimensional  continuum  mechanics  and  no  addition  for  composite  laminate 
(macro-structure).  For  example,  the  mass  densities  of  a  body  in  the  contexts  of  the  classical  con¬ 
tinuum  mechanics,  and  the  composite  laminate  (macro-structure)  will  be  denoted  by  p*  and  p, 
respectively. 

The  micro-macro  continuum  model  of  a  composite  laminate  is  illustrated  in  Figures  1  and 
2.  Figure  1  shows  a  typical  composite  laminate  (only  three  micro-structures  are  shown  in  this 
figure).  Figure  2  shows  a  shell-like  micro-structure  with  its  associated  coordinates.  This  micro¬ 
structure  is  composed  of  two  constituents  and  can  be  generalized  for  cases  of  multi-constituent 
composites. 


We  begin  the  development  of  the  kinemadcal  results  by  assuming  that  the  position  vector 
of  a  particle  P*  of  r  representative  element  (k*  micro-structure),  i.e.,  p*(0a,63(k),4»t)  the 
present  configuration  has  the  form 

p* = r^e3®,# + 4d(ea,e3®,t)  (k  =  i,..4i)  (2.i) 

where  r  is  the  position  vector  for  the  surface  4*0  and  d  is  the  director  field.  O3®,  at  this  point, 
is  an  identifier  for  the  k*  micro-structure.  Greek  super-  or  subscripts  will  assume  values  of  1 
and  2  only.  The  dual  of  (2. 1)  in  a  reference  configuration  is  given  by 

p* = tf^e3®) + 4D(ea,e3<k>)  (2.2) 

If  the  reference  configuration  is  taken  to  be  the  initial  configuration  at  time  t  *  0,  we  obtain 

pV.e300,  W)  =  r(0a,e3(k),O)  +  4d(0a,e3®,O) 

*  R(e°,e3<k)) + 4D(0°,e3®) = p**©0,©3®#  (2.3) 

The  velocity  vector  v*  of  the  three-dimensional  shell-like  micro-structure  at  time  t  is  given 
by 


v*  =  =  p*(ea,e3«^,t)  (2.4) 

where  a  superposed  dot  denotes  the  material  time  derivative,  holding  0°  and  4  fixed.  From  (2.1) 
and  (2.4)  we  obtain 

v*  =  v  +  4w  (2.5) 


where 
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v  =  r  ,  w  =  d 


(2.6) 


The  base  vectors  for  the  micro-  and  macro-structures  are  denoted  by  g*  and  gj,  respectively, 
and  we  have 


2* = ifil  e'= ^e! 

80  ae»  ’  83  34 


(2.7) 


3p* 

Bo=^7  _  .  83  =  -^- 


.ifil 

3e°!4=«  ’ 


34  15=0 


Using  (2.1)  and  (2.7)  we  obtain  the  following  relations  between  gj*  and  gj 


8>8a  +  4da 

g3*  =  g3sd 


(2.8) 


where  (  )a  denotes  partial  differentiation  with  respect  to  6°. 

By  a  smoothing  assumption  we  suggest  the  existence  of  continuous  vector  functions 
gjfB0,©3)  for  the  macro-structure  with  the  following  property 


giO0,^)^^  *  gjfe0,^  (2.9) 

where  g(0a,03(k))  are  defined  according  to  (2.7)2.  A  similar  smoothing  assumption  is  also  made 
for  the  director  d  which  we  like  to  attach  to  every  point  of  the  macro-structure.  Based  on  the 
smoothing  assumptions  we  can  write  (2.8)i  as  follows 


8a  ~  go"*-  48k  (3  *  a) 

where  {  }  stands  for  the  Christoffel  symbol  of  the  second  kind  and  is  defined  as 


(3ko)  =  ^t3oj]  =  ygki( 


dg3a 
30i  ' 


(2.10) 
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The  following  relations  can  also  be  derived  between  die  components  of  metric  tensors 
gjsft-g'andgysgi-gj 


gi*fa&  +  5l(3ka)g0k+  bkg)8oatl  +  $2{3kaH3jg}8kj 

&x3  =  ga3  +  $bka}&3 
83*3  =  833 

which  after  simplification  and  linearization  in  terms  of  %  reduce  to 

8i*8op  +  £80*3 
§o3  =  8o3  +  Y  ^633.0 


(2.11) 


(2.12) 


833  *  833 

The  determinants  of  metric  tensors  gj*  and  gjj  are  also  related  according  to  the  following 
relation 


g*  =  g  +  5A  (2.13) 

where 

g*  *  det(gjj)  ,  g  =  det(gy) 


8iu8l2j833,l 

811  812  813 

A  = 

812  822  823 

+ 

81248223833,2 

813  823  833 

813  823  833 

and  the  final  result  has  been  linearized  in  terms  of  £. 

We  recall  that  the  director  d  is  the  same  as  g3  and  therefore  when  referred  to  the  base  vec¬ 
tors  gi  it  has  only  one  non-zero  component,  namely  d3  =  1,  so  we  can  write 
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d  =  d1|i  .  d°*0  .  d^l 
djBgydi  ,  (i=l,2,3) 


(2.15) 


where  dj  and  d1  denote  the  covariant  and  contravariant  components  of  d  referred  to  g1  and  g,, 
respectively.  The  gradient  of  the  director  d  may  be  obtained  as  follows 


The  vertical  bar  ( I )  denotes  covariant  differentiation  with  respect  to  gy.  For  convenience  we 
introduce  the  notations 


Xy  =  gi  •  d  j  —  di  |  j 

fc-W-dj-dSj 

From  (2.17)  it  is  clear  that 

Making  use  of  (2.17)  and  (2.16)  we  have: 

Ag-fc^kj-RW 


(2.17) 


(2.18) 


(2.19) 


Consider  now  the  velocity  vector  v  which  can  be  written  in  the  form 


v  =  v1^  =  vjg1  (2.20) 

Again  we  make  a  smoothing  assumption  for  the  existence  of  the  vector  function  v(0a,03)  such 
that  v(0a,e3)iei=e»)  =  rCG^B3®)  after  which  we  can  define  the  gradient  of  the  velocity  field  and 
we  have 


(2.21) 
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We  now  introduce  the  notations 


vy-gi-Vj*^ 

vij  =  Ci-\j-v\j 


(2.22) 


From  (2.22)  it  is  clear  that 


vJ  =  vugj  =  vJigi 


(2.23) 


We  observe  that  both  Xy  and  Vy  represent  the  covariant  derivative  of  vector  components  and 
hence  transform  as  components  of  second  order  covariant  tensors. 


We  may  decompose  v^  into  its  symmetric  and  skew*  symmetric  parts,  i.e.. 


vii  =  v<ij)  +  vlij)sTlij  +  Q,Sj 

*lij =  v6j) =  ~2  (Vvji)  =  "Hji  (2.24) 

Also  in  view  of  (2.6),  (2.7)2  and  (2.24)!  we  may  express  gj  in  the  form 


go  =  V,a  =  (Tlka  +  <Bko)gk 

g3  =  d  =  w  =  wkgk  =  wkgk 


(2.25) 


The  gradient  of  the  director  velocity  in  8°  direction  is  obtained  by  writing 


= (nup+^kp)^ + 


=  +  X^(T]{-«oJ)  +  X^Wk]gk 

(2.26) 

The  dual  of  the  above  expressions  in  the  reference  configuration  can  be  written  easily  by 

substituting  appropriate  capital  letters  for  small  letters. 

We  now  introduce  relative  kinematic  measures  yjj  and  such  that 

Yij  =  y  C8ij  —  Gij)  —  Yji 

(2.27) 

^ij  =  ^ij  “  ^ij 

(2.28) 

where 

Gij  =  GiGj 

(2.29) 

Aij  =  [3j,i]  =  y  (G3ij  +  Gjij-Gsj.i) 

(2.30) 

Making  use  of  (2.12)  and  similar  expressions  for  the  reference  configuration  we  can  relate 
relative  kinematic  measures  y*}  of  the  micro-structure  as  follows 

YaP  =  Ypa  =  y  (&P  ”  Gap)  =  y  t(fiap  +  5  Sap, 3)  “  (Gap  +  £  Gap  3)) 

=  YaP  "*■  ~2  £(&xp,3  —  Gap,3) 

=  Yap+y4(Kap  +  ,cPo) 

Ya3  =  Ykx  =  y  CCga3  +  y  4  833,a)  -  (Ga3  +  y  4  G33.a)l 


(2.31) 


=  Y«3+“-$*3a 


(2.32) 


YJ3  =  Y33 


(2.33) 


In  obtaining  the  above  results  we  have  noted  that 


*  [3(3, a]  =  -J  (g3a^  -  g3fra)  +  J  gap3 


Vjo  =  [3 =  y  833.0 


(2.34) 


+  Kpa  =  gap3  “  °a33 


and  we  have  linearized  the  result  in  terms  of  4- 
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2.2  Basic  Field  Equations  for  Micro-  and  Macro-Structures 

The  three-dimensional  equations  of  motion  in  classical  continuum  mechanics  are  recorded 
here  for  the  k*  representative  element  (micro-structure)  in  the  present  configuration 

TWbY^pVg** 

g’xT^O 

where 

t*  =  gM/2T\  .  T*1  =  g*1/2t*^gj*  (2.38) 

The  argument  of  all  starred  functions  recorded  above  is  (©“.O3^,^)  and  the  equations  are  writ¬ 
ten  for  each  and  every  representative  element  (k  =  l,2,...,n)  which  is  assumed  to  repeat  itself  in 
the  present  model. 

Now  introduce  the  following  quantities  for  each  micro-structure: 

Composite  Stress  Vector  T: 


(2.35) 

(2.36) 

(2.37) 


Ti(0a,e3(k>,t)  ~  J  T^G0,©3^,^  (2.39) 

Composite  Stress  Couple  Vector  Sa: 

sa(e°,e3(k),t)=Y-  f  5T*°(ea,e3(k),5,t)d5  (2.40) 

?2  p 

Composite  Mass  Density  p: 


j  pVU2<i5 


o 


(2.41) 


BASE 
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pgw(z°)i^-  )  C&Vf**  <®»U>  (2.42) 

Composite  Body  Force  Density  b: 

1  J  p’b’g"^  (2.43) 

W  • 

Composite  Body  Couple  Density  c: 

M  ci¥  I  pVg*w^  llM) 

The  quantities  on  die  left-hand  side  of  equations  (2.39H2.44)  are  discrete  in  terms  of  the  vari¬ 
able  B3^  which  are  made  continuous  by  smoothing  assumptions.  The  composite  mass  density 
p0  in  the  reference  configuration  is  also  defined  as  follows: 

PoGW-^-/p.‘G‘10d4  (2.45) 

where  Po  is  the  mass  density  of  the  micro-structure  in  the  reference  configuration.  Since 
p*g*t/2  _  PoV1'2,  the  continuity  equation  for  the  macro-structure  is  readily  seen  to  be 

Pgm  =  PoGIfl  (2.46) 

Now  consider  equation  (2.36)  and  first  divide  it  by  £2  and  then  integrate  with  respect  to  £ 
from  0  to  $2  to  obtain  the  equation  for  balance  of  linear  momentum  for  the  macro-structure 

4  k 

=  T~  j  P*(v  +  £w)g*1/2d£  (2.47) 

0 
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Each  term  in  the  above  equation  can  be  represented  in  terms  of  the  quantities  defined  in  (2.39)- 
(2.44)  except  the  second  term  which  is  the  difference  between  interlaminar  stresses  above  and 
below  die  representative  element  divided  by  its  thickness  as 


9T*3 


=  rr^e®.©3^1^)  -  T*3(9a,e3(k),t)] 

S2 


(2.43) 


Now  we  postulate  the  existence  of  the  continuous  vector  function  a(8“,93,t)  whose  values  at 
e3  =  O3®  are  the  same  as  interlaminar  stresses  T*3(9a,03^\t)  and  further  approximate  (2.43)  as 
the  gradient  of  this  function  in  the  O3  direction.  With  this  in  mind  we  write  (2.47)  as 


T°ftt  +  -J—  +  pbg1/2  =  pg1/2(v  +  z*w)  (2.44) 

To  obtain  the  equation  for  balance  of  director  momentum,  (2.36)  is  multiplied  by 
integrated  from  0  to  %2  ^  divided  by  to  get 

.  «* 

■  r-  /  pVwft* + 52*)d5  (2.50) 

o 

Again  the  second  term  in  the  above  equation  can  be  written  as 

T-  1 5  TT  =  J-  RT’3)o1  -  T-  J  Tnd5  =  o  -  T3  (2.51) 

S2  o  S2  ?2  o 

As  a  result  we  have 


S® 0  +  O  -  T3  +  pg1/2c  =  pg1/2(z!v  +  z*w) 


(2.52) 


which  is  the  equation  for  balance  of  director  momentum. 
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Next,  we  consider  (2.37),  divide  it  by  \2  ***4  integrate  with  respect  to  %  from  £=o  to 
and  making  use  of  (2.8)li2  we  get 

■r  J(faxT^  +  bxT^«=0 
o 

or  (2.53) 

4>  4a 

-^-J(ga  +  $dta)x  T*°dlj  +  -^  J  d  x  T*3  d^  =  0 

and  substituting  from  (2.39)  and  (2.40)  we  obtain 

8axTa  +  diOxS°  +  dxT3  =  0  (2.54) 

which  can  also  be  written  as 

glxr  +  dAxSa  =  0  (2.55) 

This  is  the  balance  of  moment  of  momentum  for  the  macro-structure. 

Now  we  proceed  to  obtain  an  expression  for  the  specific  mechanical  energy.  Such  an 
expression  for  each  micro-structure  can  be  written  as 

p*g*1/2E*  =T#i  •  vj  (2.56) 

First,  using  (2.5)  we  write  this  equation  as 

p*g*1/2e*  =  T*°  •  (v  +  £w)a  +  T*3  •  —  (v  +  £w) 

dq 

=  T*°  •  v  a  +  £T*a  •  +  T*3  •  w  (2.57) 

Dividing  (2.57)  by  \2  and  integrating  with  respect  to  £  from  \  -  0  to  ^  =  \2  will  result  in 

4a  4a  4a  4a 

±  j  p*g*1/2e*d^  =  j-  J  T*°d!j  va+-i-|  ^°d^ •  wA +  ] T*^  •  w  (2.58) 
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We  now  define  composite  specific  mechanical  energy  for  the  representative  element  by 


pg“e=f  JpY'Vd*  (2.59) 

«  ; 

From  this  definition,  the  equation  of  continuity  and  other  definitions  (2.39)  through  (2.44)  for 
composite  quantities,  (2.58)  can  be  rewritten  as 


pgi/2e_T°.  v^  +  s°- wa  +  T3  ■  w 

Since  v  =  r,  =  (f)^  =  yja  =  ga  and  w  =  d  =  g3,  we  can  further  reduce  (2.60)  to 


(2.60) 


pg1/2e  =  T5  •  ft  +  S°  • 


(2.61) 


which  is  the  appropriate  expression  for  the  specific  mechanical  energy  of  the  macro-structure. 


2.3  Field  Equations  in  Component  Form 


We  obtained  the  following  field  equations  for  the  macro-structure  (balance  of  mass  is  not 
recorded  since  it  is  a  scalar  equation) 


T°a  +  +  pg^b  =  pgl/2(v  +  z1  w)  (2.62) 

8%  +  a  -  T3  +  pg1/2c  =  pg1/2(z*v  +  z3*)  (2.63) 

gixTi  +  diaxSo*0  (2.64) 

And  also  the  following  expression  was  derived  for  the  specific  mechanical  energy 

pg1/2e  =  V  •  fe  +  S“  •  w^  (2.65) 


By  referring  various  vector  quantities  to  the  base  gj  we  would  like  to  write  the  above  equa¬ 
tions  in  component  form.  First  write 


Ti  =  g1/2xi%i 

(2.66) 

o  =  ojgj 

(2.67) 

S°  =  g1/2  S“jgj 

(2.68) 

b  =  bjgj  ,  c  =  c*gj 

(2.69) 

where  t1^  and  S0^  are  contravariant  components  of  composite  stress  tensor  and  composite  couple 
stress  tensor ,  respectively,  &  is  the  interlaminar  stress.  Now  substitute  in  (2.62)  and  obtain 

(g1/2  *°jgj),a  +  -^3  (o’gj)  +  Pg1/2b»gj  =  pg^v*  +  z1  w^gj 


(g1/2  O^agj  +  gm  T®j  {j k  0}gk  +  +  & {jk  3}gk  +  pg1/2fcpgj 


=  pgl/2(vi  +  z!wi)gj 


or 


( g1/2  +  g1/2  j  a)  +  0»,3  +  O11  {k  j  3 }  +  p  g1/2  b» 


=  p  g1/2  (v*  +  z!w*) 


(2.70) 


Equation  (2.63)  reduces  to 


(g1/2  S°%j),0 + -  gm  T3^  +  pg1/2cjgj 

=  pg1^1^  +  Z^Jgj 


or 


(gt/2  s°j)^  +  g1/2  S0*  (k  *  a)  +  ^  *“  g1/2 +  pg1/2c^  =  p  g1/2(z!v*  +  z2^)  (2.71) 


Equation  (2.64)  can  be  rewritten  as 


gj  x  (g1/2  x'jgj)  +  X'ogj  x  {gm  Sojgj)  =  0 
or 


g1/2(xij  +  V0S“j)gi  x  gj=0  (2.72) 

since  g  *  0,  gixgj  =  £ijkgk  and  is  skew-symmetric  we  conclude  that  the  quantity  in 
parentheses  in  (2.72)  must  be  symmetric  in  i  and  j.  As  a  result,  the  conservation  of  angular 
momentum  in  component  form  is  the  symmetry  of  T®  defined  by 


rj^j+o00 

(2.73) 

T»j  -fi* 

(2.74) 

The  expression  (2.65)  for  the  specific  mechanical  energy  can  also  be  written  as 


pg^e  =  %m  •  i + g1/2  S*^gj  •  w*  icjj 

=  *  8a  +  S^Wjla  +  T3^  •  fo) 

=  •  V.a  +  T^gj  •  W  +  S«%„) 

-  g1/2(t“jvj|0  + 13^  +  S°jwj  |o) 

We  have  now  the  component  form  of  the  expression  for  mechanical  power 

Pipe  =  x0^  1o  +  t*Wj  +  S“jwjlo  (2.75) 

An  alternative  form  for  mechanical  energy  expression  is  derived  in  which  the  rates  of  rela¬ 
tive  kinematic  measures  will  appear.  Using  (2.25)i  and  (2.26),  we  rewrite  (2.75)  as 

P  =  pe  =  ^(Tlja  +  CDja)  +  T^Wj  +  SW[Xjp  +  +  Wjo)  +  *0  Wj] 

=  (t“j  +  +  (t°j  +  SWXJ)<Dj „  + 

+  (T3j  +  X|sPj)wj  (2.76) 

Recalling  (2.73)  and  using  symmetry  of  and  skew-symmetry  of  Wy  we  can  write  (2.76)  as 

P  =  ’^jnja  +  SpjXjp  +  t3jwj  +  Ta30)3O  (2.77) 

By  (2.25)!  we  have 

ga-gp  =  Tlpa  +  <«>pa  .  8p  '  8a  =  “Hap  +  ©a?  (2.78) 

Therefore 


8ap  =  8a’gp  +  go,gp  =  2rtap 


(2.79) 


In  the  last  result  we  have  used  the  definition  of  from  (2.27).  By  (2.25)lt2  we  have 

fa  83  =  1130  + ©3a  (2.81) 

fofa  =  wa  (2.82) 

Therefore 

faa  =  fag3  +  fo-fa:=ll3a  +  0>3o  +  Wa  (2-83) 

ll3o  =  fa3“(t^o  +  wo)a:2Ya3-(Cl)30  +  Wa)  (2.84) 

Again  by  (2.25)2  and  (2.27) 


o>3  =  83  -  83  =  Y  g33  (2.85) 

Vi3  -  833  =  w3  (2.86) 


and  by  (2.28) 


Xjp  =  KjP  (2.87) 

Substituting  from  (2.80),  (2.84),  (2.86)  and  (2.87)  in  (2.77)  we  get 

P  =  +  rx3{2ya3-(0)3a  +  wa)}+  sWiCjp  +  T3aw0t  +  T33y33  +  T0*^  (2.88) 

which  is  simplified  to 

P  *  +  2Ta3ya3  +  T33yj3  +  SPjK-p  (2.89) 

If  symmetry  of  T*  and  y^  is  considered  we  can  further  simplify  (2.89)  and  obtain 


P-T4  +  S»fc, 


(2.90) 
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(tij  +  X^^y+S0^ 


(2.91) 
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2.4  GeneraJ  Constitutive  Assumption  for  Elastic  Composite 


At  this  point  we  postulate  the  existence  of  specific  internal  energy  in  purely  mechanical 
theory  which  depends  cm  relative  kinematic  measures  Yy  and  as  defined  in  (2.27)  and  (2.28) 

V=¥(Yy.*ja)  (2.92) 

P  =  py  (2.93) 

By  usual  procedures  we  obtain  from  (2.91),  (2.92)  and  (2.93) 

*•«%-*•&  <2-94) 

sai  =  p^  (2.95) 

Now  the  composite  stress  vector  T1  and  the  composite  couple  stress  vector  S°  from  (2.66)  and 
(2.68)  will  be 


s“=pein<^)gJ 


(2.96) 

(2.97) 


The  coefficient  pg1/2  can  be  replaced  by  p0G1/2  by  taking  advantage  of  the  continuity  equation. 
Note  that  by  these  constitutive  relations  for  T‘  and  S°  the  balance  of  moment  of  momentum  is 
identically  satisfied. 
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2.5  Linearized  Kinematics 


For  linearized  kinematics  let 

r(0a,03(k),t)  =  R(0a,03<k))  +  eu(0a,03<k),t) 

(2.98) 

d(ea,e3(k),t) = D(0°,e3(k)) + aqp&Mj) 

(2.99) 

•  • 

v  =  f  =  eu  ,w  =  d  =  e8 

(2.100) 

where  e  is  a  non-dimensional  parameter.  The  motion  of  the  macro-structure  describes 
infinitesimal  deformation  if  the  magnitude  of  the  gradient  of  the  displacement  vector  eu  and  the 
magnitude  of  the  director  displacement  vector  eS  are  of  the  order  of  E  «  1  such  that  in  the  fol¬ 
lowing  developments  we  can  only  retain  terms  which  are  linear  in  e.  The  base  vectors  gj  are 
found  from  (2.7)2  as: 


§a~R.0+eu0  (2.101) 

g3  =  d  =  D  +  e5  (2.102) 

The  corresponding  vectors  in  reference  configuration  are: 

G0  =  R0  ,  G3  =  D  (2.103) 

We  now  proceed  to  obtain  the  relative  kinematic  measures  yy  and  k^.  Using  (2.103)2  and 
(2.101)  together  with  the  definition  of  g^  and  G^  we  write 

goj)  =  (G0  +  eu  „)  •  (Gp  +  eu  p)  =  Gop  +  e(G0  •  u^p  +  •  Gp)  +  0(e2)  (2. 104) 

where  0(e2)  denotes  terms  of  order  e2  in  displacement  gradient,  where 

Ga  •  up  +  u  a  •  Gp  =  Ga  •  iPipgj  +  iPiagj  •  Gp 

=  Ga  •  uT,p»r+  G„  •  u’,^3  +  u’log,  •  Gp  +  u5 10g3  •  Gp 
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=  u^pGj,  •  (Gy+ eu,y)  +  u3ipGa  •  (D  +  e8)  +  uT|a(GY+  euy)  •  Gp 

+  u3  ia(D  +  e5)  •  Gp  (2. 105) 

Retaining  terms  which  are  of  the  order  of  unity  in  (2.105)  and  substituting  the  result  in  (2.104) 
we  find 

Yap  =  J  (gafi  -  Gap)  =  J  («aip  +  “pia)  +  *J  (u?pDa  +  U?aDp)  (2.106) 

Here  covariant  differentiation  is  supposed  to  be  performed  with  respect  to  the  metric  Gy  of  the 
reference  configuration  and  instead  of  eu  we  have  used  u  with  the  same  assumptions  made  for 
linearization.  Similarly  we  can  write 


go3  =  (Ga  +  eua)  •  (G3  +  e8)  =  Goj  +  c(Ga  •  8  +  ua  •  G3)  +  0(e2) 

=  0,3,3  +  £(5a+  u31a)  +  CKe2)  (2.107) 

Again  using  8  instead  of  e8  with  the  same  interpretation  we  obtain 

Ya3  =  Y3a=y  (&x3“Go3)  =  Y(8a  +  U3io)  (2.108) 

To  find  Y33  we  write 

g33  =  (G3  +  e8)  •  (G3  +  e8)  —  G33  +  2e53  +  CKe2)  (2. 109) 

Yj3  =  y  (833  -  G33)  =  83  (2. 1 10) 

As  for  the  measures  Kj„  we  proceed  as  follows 

^ap  =  8a  ’  d,p  =  (Ga+£u  o)  •  (Dp+eSp)  =  Aop  +  e(iP  1  agj 1  Dp  +  Ga  •  5*  |  pgj)  +  CKe2)  (2. 1 1 1 ) 

where: 
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«*!«& *  D.p  *  »Tio(GT-t-  eu.Y)  *  Dj  +  u3|0(G3  +  e8)  • 

=  UY)oAyp  +  U3  Ia^3p  +  0(e) 

Ga  •  6*ipgj  =  5^1  pGa  •  (GY+ eutT)  +  S3ipGa  •  (G3  +  e8) 

=  8aip  +  83ipDa  +  0(e) 

Substituting  these  results  in  (2.1 1 1)  and  using  the  definition  of  we  get 

*ap  =  ^t*p  ”  ^op  *  ^la^jP  +  ^oip  +  83,pDa 
Now  we  obtain  an  expression  for 

X3 a  =  g3 '  d,a  =  (G3  +  e5)  •  (D^  +  e5t0l) 

Xqb  —  +  e(G3  •  5a  +  8  •  Dj0t)  +  CKe2) 

We  simplify  each  term  separately 

G3  •  8  0  =  G3  •  (8j(agj)  =  57iaG3  •  (Gy -f-  eUy) 

+  83|0G3  •  (D  +  eS)  =  ST|0tDy+  83)aD3  +  0(e) 

=  8iloDj  +  0 (e) 

8  •  D  a  =  (8>gj)  •  (D*  loGk)  =  (5igj)  •  D3  laG3 

=  A3(5YgT+83g3)-G3 

=  A3[8T(Gy+ euy)  +  83(G3  +  e8)]  •  D 

=  0(e)  +  A^STDy + A383D3  =  A3^  +  0(e) 

However,  since  D°  =  0  and  D3  =  1 


(2.112) 


(2.113) 


(2.114) 


(2.115) 
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5t>j  =  SjD1  =  53  (2.116) 

Substituting  from  (2.114),  (2.115)  and  (2.116)  in  (2.113)  and  using  previous  notation  for  5  we 
obtain 


K3o  =  X3a-A3o  =  8*|0Dj  +  A|83  =  8310  +  ^83  (2.117) 

To  recapitulate  the  relative  kinematic  measures  in  linear  theory  are: 

Yap  =  J  (“alp  +  Uplo  +  U3loDp  +  U^pDJ 

Ya3 ~  Y3a  =  (u3la+  &a) 

Yj3  =  «3  (2.118) 

M“APujla+8alp+5  tpDa 
K3a  =  83la  +  A^53 

For  a  composite  with  initially  flat  plates  we  can  always  choose  our  base  vectors  Gj  such  that 
Gy  =  G,J  =  5y  and  as  a  result  D°  =  Da  =  0  and  if  we  confine  ourselves  to  small  deformations,  then 
all  Christoffel  symbols  vanish  and  covariant  differentiations  reduce  to  partial  differentiations  and 
equations  (2.1 18)  for  relative  kinematic  measures  will  further  reduce  to 


Yap  =  y  (Ua,p  +  Up,a) 

Ya3  =  Y(U3,a  +  5cx) 

(2.119) 

YJ3  =  53 

Kap  =  8a,p  *  K3a  =  83,a 

In  writing  these  relations  it  has  been  noted  that  Da  =  0,  D3  =  1  and  s  0. 
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It  is  also  desirable  to  find  a  relation  between  g  and  G,  determinants  of  metric  tensors  in 
present  and  reference  configurations,  in  the  linear  theory.  We  recall  the  following  relations 

g1/2  =  8i*g2*S3  >  G^-GjxGj-Gj 

Si  x  g2  =  (Gj  +  eutl)  x  (G2  +  eu^) 

=  Gj  x  G2  +  £(uj  x  G2  +  Gj  x  +  0(c^) 

[818283]  =  [Gi  x  G2 + e(u!  xGj  +  GjX  u^]  •  [G3  +  e5]  +  0(e2) 

=  G1/2  +  c[G]  x  G2  •  8  +  Uj  •  G2  x  G3  +  G3  x  Gj  *  +  0(e^) 

=  Gm  +  eG1/2[G3  •  5  +  u j  •  G1  +  G2 •  uj  +  CKe2) 

Retaining  terms  which  are  of  the  order  of  e  and  using  the  previous  notations  for  u  and  5  we  get 

(-£-)W=l  +  83  +  u°,a  (2.120) 

G 

Now  the  equation  for  balance  of  mass  will  readily  reduce  to 

Po  =  P(^-)ln  =  P(l  +  S3  +  ““ia)  (2.121) 

and  since  in  linear  theory  displacement  vector  u  and  director  displacement  S  satisfy  linearity 
assumptions  we  obtain 


p  =  p0(l-83-uala) 


(2.122) 
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2.6  Linearized  Field  Equations 

We  use  pertinent  results  from  linear  kinematics  and  usual  procedures  for  linearization  to 
write  the  field  equations  in  linear  theory.  It  should  be  recalled  in  such  analysis  that  g  is  replaced 
by  G,  p  by  p0  and  Christoffel  symbols  are  calculated  with  respect  to  Gy.  By  omitting  the  details, 
the  linear  version  of  the  equations  of  modem  are  recorded  here 


(G1/2  +  Gm  f*  (k  J  «}  +  ai  +  o*  {k  J  3)  +  p0  Gm  b*  =  p0  Gm  (tf  +  z1^)  (2. 123) 

(G172  S°j).0  +  Gm  S'*  {k* «}  +  oi  - Gw  4-  p0  G1/2  c*  -  =  p0 Gm  (z1*?  +  z2#)  (2.124) 

T*j  =  x‘j  +  =  P*  (2.125) 

For  a  composite  with  initially  flat  plies  further  simplification  can  be  made.  As  mentioned 

earlier,  G  =  1  and  all  Christoffel  symbols  vanish  identically.  The  resulting  balance  equations  for 
such  a  situation  will  be 


^.0  +  0*3  +  Pob*  =  Po^  +  z1^) 

(2.126) 

S“i  0  +  d  -  t*  +  p0oi  =  Po(z1iP  +  z2#) 

(2.127) 

T*i  =  tij  =  t* 

(2.128) 
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3.0  CONSTITUTIVE  RELATIONS  FOR  LINEAR  ELASTICITY 


For  the  representative  micro- structure  let 


a  =  U 


where 


(3.1) 


0<S<Si 


02) 


and  (a  =  1,2)  are  material  constants  in  associated  layers.  Now  we  proceed  to  calculate  T1 
and  S°  defined  in  (2.39)  and  (2.40).  First  we  recall  that  T*J  =  gm  g*1/2  =  gl/2  (l  +  -^  — ), 

2  g 

Sy  -  By  +  £d>r  g3  -  g3  and  for  brevity  we  omit  the  index  a  in  relations  (3.1)  and  (3.2) 


T  J  TX  =  £  J  g-1" 

=  Ti  l  8”l/2  CfSUy*/g&  =  IT  i  +  2c*"i**5  +  c‘i33'ft3)SJ‘d^ 

Substitute  from  (2.31).  (2.32)  and  (2.33)  in  the  above  relations  and  get 


, 

V  =  f-  j  gm  (^op + 2cija3ya3  +  cij33%3)$  d£ 

S2  o 


$2 


+  f  J  8*lfl  f  Kj^Kg-d? 
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+  *3a^  Jc*>“V1/2g*$d$  (3-3) 

We  calculate  each  integral  separately,  noting  that 

+  (3.4) 

g* lfl  g-  ,  gl/2  ^  +  A  fe)  =  g'fl  (,  +  A  5)fe  (3.5) 

The  first  integral  in  (3.3)  is 

*;<*=£/  (c*vi/2  g; + g->«  0.6) 

The  first  term  of  (3.6)  from  (3.4)  is  equal  to 

g,B  ^  i !  «** + g,s  (£  * + d-^i  I  ^ + if  { *’  <3-7) 

and  its  second  term  can  be  written  as,  from  (3.5), 

g'°  g3(-£  J  c°“d|  ^J3Udy  (3.8) 

Combining  (3.7)  and  (3.8)  we  rewrite  (3.6)  as 

im  ij  J  «"*  +  gj  J  5cljud5  +  gw  d,T  J 

+^d4{^  <3-9) 

The  second  integral  in  (3.3)  is 
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£  J  <?**  i  w  %;  \  <* = £  f  *•»  *• + 

The  first  term  of  (3.10)  from  (3.4)  is 

*' 1/2  (r  £  J  5<=w<15 +  *r + ^  J 

+  2?»fe{^Cn°>^  <311> 

and  its  second  term  from  (3.5)  is 

g1/2  fc(-£  J  ^  “  J  ^)  (3.12) 

Combining  (3.1 1)  and  (3.12)  we  rewrite  (3.10)  as 


E'n  8i  ¥ !  ?cWld4  +  gj  ¥ ! 42  c‘Wd4 + gW  ¥  { 42ci,nPd4 


+ 


A 

2gm 


(3.13) 


The  third  integral  in  (3.3)  is 


-jT  J  £clja3g*1/2  g*d|  =  g1/2  gj  ~  J  £cija3d£  +  gj  J  52  cij0t3d^ 

+  im dY ^  J  4W*  +  ^  d.T -1-  J  (3.14) 

The  last  result  was  written  by  noting  the  development  in  (2.13).  The  results  in  (3.9),  (3.13)  and 
(3.14)  can  further  be  simplified  by  recalling  (2.9)3,  namely  d  .f  =  X^gj  and  using  the  following 


definitions  and  results. 
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Define 


mss5i/^2<1 

f  0<5<$i 

*={  »2  tl<5<42 

Then 

-5*  *k  . 

^ /^k®d^  =  -^r[inai  +  (~-m)a2]  k#-l 


Now  (3.9)  is  equal  to 

gw  %(mc}^  +  (l-m)^  + -^  [mc!i“ + (— -  m)c^ 

4g  m 

+  M  -y  [loci’"*  (■—  -m)c||t'|  +  (-^  -m)c|>ul) 

=  tm  gj(0  +  ^•>'0=!*'+  yl{(l  + 

(3.13)  can  also  be  written  as 

g1/2  gj{[nac|j°^  +  (^-  -  m)c|aP]  +  A  ll_  [me!*3#  +  (A  - 
m  2  2g  3  mz 

t2  k  3 

+  y  [mcfl“P  +  (A  _  m^hapj  +  A  \i  A  [mc|iap  +  (J__  m)^?] } 


(3.15) 

(3.16) 

(3.17) 


(3.18) 


6g 


x  4id-m)  ,,  A 
+  ^-(1  +  m 


+  ij£  > . 

3g  m  3  m 


m 


+  2  M(-L-m)]cJW])  (j.19) 

8  g  nr 

The  expression  for  (3.14)  is  exactly  the  same  as  (3.19)  except  that  (a,p)  in  (3.19)  should  be 
replaced  by  (a,3).  These  results  should  be  incorporated  in  (3.3)  to  find  an  expression  for  T*. 
Due  to  the  presence  of  die  factor  g112  gj  in  all  these  expressions  and  also  the  equality  (2.66),  we 
can  find  the  constitutive  relation  for  t^.  However  before  doing  so  we  exploit  the  symmetry  of 
&V  to  further  simplify  (3.3).  Since  c1*^  »  c^°  we  can  write 


Ci^1Cpa  =  C«PaKpa  =  Ci^Kap  (3.20) 

Therefore, 


=  y  C^OCap  +  KpJ  (3.21) 

In  view  of  (3.21),  now  we  write  (3.3)  as 

.  5.  .  ,  «■  . 

T1  =  Tu  J  C,jk/  g*1/2gj*  d^  +  Kta  -T-  J  cijfa  gm  gj*  4  (3.22) 

*•2  0  0 

With  the  explanation  presented  above,  now  we  write  the  constitutive  relation  for 


t“  =  {<1  +  -^■)mc}ik'  +y^(1  + 

+  (l-m)(l  +  ^  — X^+  ^'(1—  X|  (1  +  m  +  ^ 

4m  g  2m  '  3g  m 
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2m  3g  m  3  m 

+  *r  (-V  -  mHc^lK*,  (3.23) 

8  g  nr 

If  we  further  restrict  ourselves  to  small  deformations  of  a  composite  with  initially  flat  plies, 
equation  (3.23)  can  be  simplified  to 

T«=  (mcJ»,+(l-tti)c^,)-&/+  {ooc}>**  +  -l=si  (324) 

2  m 

Of  course  in  the  above  equations  no  distinction  should  be  made  between  covariant  and  contra- 
variant  components  of  tensors.  Using  (2.114),  equation  (3.24)  can  be  written  in  terms  of  dis¬ 
placement  vector  u  and  director  displacement  vector  S,  hence 

%  =  (meg  +  (l-mfcgjTk  +  ~  {mcg>  +  <$>  Jk*, 

=  {mcigjj  +  (l-mjcg^)  j  (u^p  +  uPt0)  +  2(mc^  +  (l-m)cg^)(u3>0+  8^/2 

+  {mc$  +  (l-m^)^  +  y  (mc$>  +  —p  cg>}8/>0  (3.25) 

Using  the  symmetry  of  Cjjy’s  (3.25)  can  be  written  as 

*ij  =  +  (l-mjcg^luo.p  +  {mcija3  +  (l-nOcg^u^ 

+  {mc^  +  (l-mJc^JSo  +  {mCj^  +  (l-m)c@}83 

+  i.(mcXS  +  il^tc®)Sla 

=  {mciji  +  (l-mJcgiJu^a  +  {mc^  +  (l-mjc^)^ 


(3.26) 


The  same  steps  can  be  followed  to  calculate  S°  and  we  record  the  procedure  here  KS 

i  4i 

5,  5a 

]  tfW^dt  +  n*  ± (3.27) 

This  is  basically  the  same  result  as  (3.22)  except  that  i  has  been  replaced  by  a  and  the  integrand 
has  been  multiplied  by  %.  The  first  integral  can  be  reduced  to  the  following  by  referring  to  (3. 1 8) 

t'°  *jt%  [mef* + (—  -  m)  c?i“]  +  [mcf^  +  ( ±r  -  m*?**)  + 
z  m  og 

fi  2  3 

H  +  (~Y  -  ^  [mcf**  +  (-~  m)cf^ }  (3.28) 

J  m  og  m 

Similarly  the  second  integral  is  simplified  and  by  reference  to  (3.19)  the  result  is 


^  2  j 

g1/2  gj{ ~~  [mcf^  +  (~-  m)c|^]  +  [mcfjip  +  (4r  -  m)cf^] 

J  m  og  nr 

+  »4'y[,nc1“^  +  (-^-  -  m)cf#]  +  XiM[mc  «)«  +  (•_  m)ca^] , 

=  tf  «*  gj((f  ♦  — )cf^ + Xj4,(i  +  M)nKl^ 


Again  it  is  seen  that  because  of  the  factor  g1/2  gj  and  the  relation  S°  =  Saj  g1/2  gj  we  can  readily 
write  the  constitutive  relation  for  S°J,  the  result  is 


S“i  =  5,  ( i  <mc?*)(l  +  +  + 


+  ~ 


1  , 1-m2  .  1-m3 


1-m3  .  1-m4 


m 


3  8g 


)cf^)Tk i 


Hf  ((}  ♦  +<1^  +  1T 


♦«**#+  10g 


$1*  l-m5 


m 


(3.30) 


This  is  the  general  constitutive  relation  for  S°j  in  linear  elasticity.  If,  as  before,  we  confine  our¬ 
selves  to  small  deformations  of  a  composite  with  initially  fiat  plies  (3.30)  can  be  simplified  to 


1  5,  (me?*  +  c?*)*,  +  i  tfoncfl*  +  (3.31) 

2  m3  m 

with  no  distinction  between  covariant  and  contravariant  tensors.  Once  written  in  terms  of  dis¬ 
placement  vector  and  director  displacement  and  simplified  as  done  in  obtaining  (3.26)  we  get 


Soi  =  i  4i<mc?)f>k  +  ±zt  +  i  5,(mc?*  +  c?^ 

L  m  r  i  m 

+  \  +  -^1  (3.32) 

3  rrr 

As  the  results  of  this  section  indicate,  even  in  the  simplest  cases  of  small  deformations  of 
an  initially  flat  composite  (composites  with  flat  plies)  higher  gradients  of  displacement  vector 
become  significant  and  they  appear  in  the  constitutive  relations  for  composite  stress  and  compo¬ 
site  stress  couple.  As  defined  by  (3.13)  m  is  of  the  order  of  unity,  and  £2  are  usually  small 
lengths;  however  their  products  with  components  of  Cyu  and  even  the  product  of  their  higher 
powers  with  elastic  constants  may  be  indeed  significant  quantities,  in  which  case  we  get 
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contributions  to  ty  and  Saj.  In  the  trivial  case  m  =  1 ,  2^  =  £2  -» 0  we  get  x,j  =  c^y^,  Saj  ■  0  and 
the  equations  of  linear  elasticity  are  recovered. 
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4.0  COMPLETE  THEORY  FOR  LINEAR  ELASTIC  COMPOSITE  LAMINATES 

The  results  of  sections  (2)  and  (3)  are  combined  to  obtain  the  complete  equations  for  a 
linear  elastic  composite  laminate.  However,  before  doing  so  we  should  derive  appropriate 
expression  for  p0,  z1  and  z2.  As  before,  we  assume  that  the  representative  micro-structure  is 
composed  of  two  homogeneous  layers  with  respective  densities  and  p2  in  the  reference 
configuration  (pj  and  p2  are  constants).  Recalling  equations  (2.45)  and  (2.46)  we  write 

pg1/2  =  peG1/2  =  /  p*G*1/2d£  (4.1) 

Now  by  (2.13)  we  have 


G*1/2  =  G1/2(1  +  ^-)  (4.2)* 

where  A  is  understood  to  be  the  sum  of  two  determinants  similar  to  those  expressed  in  (2.14) 
except  for  substituting  gy  by  Gy.  We  have 


Po 


Pi 

P2 


0<5<5i 

$1<$<$2 


Substituting  (4.2)  and  (4.3)  in  (4.1)  and  using  (3.7)  we  set 


(4.3) 


tiA  1+m 

Po = a  +  “-)mp!  +  (l  +  ^“7  4iA)(l-m)p2 
Of  course  the  composite  mass  density  p  is  related  to  p0  through  the  equation  (2.122). 


(4.4) 


We  proceed  similarly  to  calculate  z1  and  z2  using  their  definitions  in  (2.42) 


S2  o 


(4.5) 


*  Here  again  the  result  has  been  linearized  in  terms  of 
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ta 


p*«¥ = p.0^  -  ^  /  5W 

After  substituting  from  (4.2)  and  (4.3)  in  (4.5)  and  (4.6)  and  using  (3.17)  we  get 


Po* 


l  r/i  .  x  .  /  1-m2  .  l-m: 

=  T[(1  +  lo)niP>  +  (—  +  1g^ 


■)pJ 


.  _j  5?  .  .  ,1-m3  .  3^A  , 

Po^"T[<l  lG-)mpl  +  (-rf-+-SG- 

For  a  composite  with  initially  flat  plies  (4.4),  (4.7)  and  (4.8)  ait  reduced  respectively  to 


(4.6) 


(4.7) 

(4.8) 


p0  =  mpj  +  (l-m)p2 


(4.9) 


i  §1  ,  ,  1-m2  . 

Poz  =  ~T  (mpi  +  — —  P2) 
JL  m 

_2  5?  ,  .  1-m3  N 

P0r  =  —  (mp!  +  — —  pj) 
j  m 


(4-10) 


(4.11) 


To  formulate  the  complete  theory  it  is  also  worthwhile  to  derive  a  relation  between  the 
director  displacement  5  and  the  gradient  of  displacement  vector  u  in  the  (^-direction.  In  order  to 
derive  such  a  relation  we  enforce  the  continuity  of  the  position  vectors  p*  and  P*  between  two 
adjacent  micro-structures.  Recalling  (2.1)  and  (2.2),  we  have  the  following  relations  for  the  k* 
micro-structure: 


p*(e°,e3(k>,4,t)  =  KG'se^.t)  +  ^dce^e^t)  (4.12) 

P*(Ga,03(k),4)  =  R(G°,03(k))  +  SDCG^G300)  (4.13) 

Now  in  order  that  position  vectors  p*  and  P*  be  continuous  on  the  common  surface  between  k* 
and  (k+l)*  micro-structures  we  should  have 

P*(G“,  e30^  =  P'^.e3^1^)  (4.14) 


BASE 


P 'co3®**!)  -  p,(e“,e*fc*1>AO 


(4.15) 


Using  (4.12)  and  (4.13)  we  can  write  (4.14)  and  (4.13)  as 

R(e°,e3<k+1))  *  R(e°,e3<k)) + yxe0,©30^  (4.i6) 

r(e»te^+1).t)  =  r(0°,03<k),t)  +  (4.17) 

Recalling  (2.98)  and  (2.99)  and  identifying  eu  and  e5  with  u  and  5  as  before  we  conclude  from 
(4.16)  and  (4.17)  the  following 

u(0°,e3<k+1),t)  =  u(0°,03<k),t)  +  ^^.t)  (4.18) 

or 

^.OW.t)  =  ~~  {u(0a,e3<k+1),t)  -  u(0°,03(k),t)]  (4.19) 

S2 

By  smoothing  assumptions  and  approximating  the  right-hand  side  of  equation  (4.19)  as  the  gra¬ 
dient  of  the  displacement  vector  in  the  O3  direction  we  have 

8(0g,e (4-2°) 

In  component  form  we  have 

6  =  5igj  =  (Uigj)^  =  ui|3gj  (4.21) 

or 

8»  =  ui|3  ,  6j  =  Uj,3  (4.22) 

For  a  composite  with  initially  flat  plies  the  equation  (4.22)  reduces  to 

5i  =  UjJ 


(4.23) 


With  this  simplification,  equations  (2. 1 19)  reduce  to 


‘ftj*y(uiu+uii)  (4*24) 

Kja  =  Uj,a3  (4.25) 

Using  (4.22),  equations  (2.121)  and  (2.125)  are  also  reduced  to 

P<,  =  (X^)m-Pa+rf|i>  (4.26) 

p»p0(l-u>,P  (4.27) 


The  constitutive  relations  (3.26)  and  (3.32)  for  ty  and  Saj  for  a  flat  composite  are  also  further 
simplified  by  using  (4.23) 

*y  =  (meg  +  (l-m)cg)uu  +  {megj  +  eg,}  uta3  (4.28) 

XU  / 


and 


=  T  c®)uu+  4  +  -Mf-  c&)utfl3  (4.29) 

Now  we  can  write  the  field  equations  (2.113)  and  (2.124)  in  tenns  of  the  displacement  vector  u 
and  its  gradients.  It  should  be  recalled  that  the  resulting  equations  are  the  linearized  field  equa¬ 
tion  for  small  deformations  of  a  composite  with  initially  flat  plies.  These  are  the  counterpart  of 
the  Navier-Cauchy  equations  in  linear  elasticity.  The  appropriate  equations  for  a  general  com¬ 
posite  will  be  derived  in  later  chapters.  Using  (4.9)-(4.11),  (4.23),  (4.28)  and  (4.29)  we  write 
(2.123)  and  (2.124)  as 
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(mc$/  +  +  (rocjffip  +  cffip)  -y  utpy. 


+  0j3  +  [mpi  +  (l-mjpjbj  =  [mpj  +  (l-mjpjiij 
2 

+  y  5i(mpi  +  Pjiiy  (4.30) 

and 

^  2 

\  5i  {“*$/+  c^Ji}uk.ta+  y  [wcSfe  +  yf~  4&)Uk.|Ba 

t  2 

+  Oj  -  [mc$  +  (l-m)cg,}uw -  y  {mc|^  +  -yy  cg^}uta3 

2 

+  [mpj  +  (l-m)pJCj  =  \  ^(mp!  +  p^Uj 

z  m 

1  1  m3 

+  i^(mPl+-L:T-P2)uJ.3  (4.31) 

j  m 

At  this  point  we  may  notice  that  an  ordinary  continuum  (a  single  material  continuum)  can  be 
regarded  as  the  limiting  case  of  a  composite  laminate  when  =  £2  -*  0-  Therefore,  we  may 
anticipate  to  derive  the  equations  of  linear  elasticity  by  letting  m  =  1  and  -*  0  in  equations 
(4.30)  and  (4.31).  Doing  so,  equation  (4.30)  reduces  to 


cajkfikja  +  °>3  +  Pbj  =  Puj  (4.32) 

where  subscript  and  superscript  1  are  dropped  because  we  have  only  one  material.  To  simplify 
equation  (4.31),  first  we  recall  the  definition  of  c  in  equation  (2.44)  and  notice  that  by  the  mean- 
value  theorem,  c  ->  0  as  %2  “*  0.  hence 


°j  “  C3jkiuk,/ = 0 


(4.33) 
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Substituting  for  Oj  from  (4.33)  in  (4.32)  we  get 

cct*iU*Ja + c3jk^kj3  +  pbj  *  piij  (4.34) 

and  combining  the  first  and  the  second  terms  we  get 

CijuUk.fi  +  pbj  =  piij  (4.35) 

For  a  completely  isotropic  continuum 

CjjVf  as  +  |l(5jk5y  +  5y5jk)  (4.36) 

where  X  and  p  are  Lame  constants  and  (4.35)  reduces  to: 

puy  +  (X+iOujjj  +  pbj  =  piij  (4.37) 

which  are  the  equations  of  motion  for  an  isotropic  media. 

For  the  case  of  the  composite  laminate  we  can  also  eliminate  dj  between  equations  (4.30) 
and  (4.31)  to  obtain  the  appropriate  equations  for  displacement  vector  u.  First  we  do  this  for  a 
static  problem  with  no  body  force.  For  such  a  case  we  let 

b  =  c  =  ii'=  0  (4.38) 

in  equations  (4.30)  and  (4.31),  hence 

(mc^  +  (l-m)c$/}uUa  +  {mcgp  +  c$p)  y  u^p^  +  0^  =  0 
and 
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y  {“C$/  +  yy  Cofi/Kto*  Y  (mcaj^  +  '^-  C®p)»»tP3o 

+  0j-  {mcjiJi  +  O-inJc^lUkj-—  (mcffig  +  -1~— -  c^a}ukta3  =  0 

Z  111 

Eliminating  oj  between  these  equations,  we  get 
{mc^  +  (l-mjcg/ju^ta  +  {mc$p  +  -yy  c^p}  y 

+  {mcj$+  (l-m)c|^/}uk>i3  +  y  {mc$«+  1~n*  c^Ju^s 

Z  111 

-  y  {mCaj^i  +  YY~  c®  lUk.faa  -  y  {mci^  +  cgp}utPa33  =  0  (4.39) 

By  combining  the  first  and  third  and  also  the  second  and  fourth  terms  of  the  equation  (4.39),  we 
get 

{mc$+  (l-m^Ju^  +  y  {mc^  +  yy-  c^Ju^ 

-  y  {mc^+  -yy  c^}utJa3  -  y  {mc^  +  yy  c^)utpa33  =  0  (4.40) 

This  is  a  fourth  order  partial  differential  equation  for  displacement  vector  u.  Now  we  apply  this 
equation  to  a  composite  laminate  whose  micro-structure  is  composed  of  two  isotropic  layers. 
For  such  a  case  we  can  write 

ciji)  =  +  StfSjk)  (4.41) 

cijk)  =  ^2^ij8k/  +  +  (4.42) 

where  Xj  and  m  (i  =  1,2)  are  Lame’s  constants  for  the  respective  layers.  Introducing  equations 
(4.41)  and  (4.42)  in  (4.40)  we  get  the  following  for  each  term: 
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{mc$j + (l-mjcijjjut/i  =  +  u jjM 

+  (l-inMXjUkj^  +  ji^Ukjc +  jd) 

-  {mCXj  +  m)  +  (l-mKXjfHijJJuytj  +  {m^  + 

(mcyka  +  ^  c§n)ukja3  *  m{Xi^j5jBB  +  |ii(^^ja  +  \s5jk)}Ukjo3 

+  ~  +  M&iSjQ  +  ^aSjk))t^ia3 

=  ^{^l^jSpaU^jaj}  +  mjIifS^SjaU^^  +  SflaS^Uj^j} 

+  m  ^2(^j5poup,ia3)  +  m  ^(S^SjoU^o^  +  S^S^U^^) 
-  {“^1  +  ^2)uoJa3  +  (“Vl  +  —  ™  *  ^}5jaukjtn3 

+  ltrqn  +  42)Uj,aa3 

{mc^j^/  +  m  ca^)utta3  =  m{^lS0j5k/ +  Hl(Sak8j/  +  Sq/I^Ju^^ 

+  ”  {X^SajSy  +  ^(S^Sj/  +  8a/8^) }utta3 

=  mX.i(8ajUk^a3)  +  +  uj,aa3) 

.  (1-m2)  ,  /C  x  .  (1-m2)  „  , 

+  m  ^2^ojUWc«3^ +  m  ^2(uaja3 +  uj,aa3) 

=  (mXj  +  --^  ^  ^2>5ajUtka3  +  (m^!  +  *  H2>(u0ja3  +  ujaa3) 


(4.43) 


(4.44) 


(4.45) 


BASE 


(®Co&  +  cajJp)uMa33  =  m*Xl5aAp  +  M5*5*  +  fiapS^Ju^j 


1— 

+  -“2— '  {^ajSkP  +  M2(5ak8jp  +  6^6jk)}utpa33 


*  mXl80jUp<pa33  +  D3Ji1(8jpUa,pa33  +  Uj.0033) 


1-nr 


m‘ 


1-m 


^2®ajup,pa33  + - -  M2(8j3ua.Pa33  +  ujJ«*33) 


m 


1-m3 


1-m3 


-  (mA-!  +  X2)80jUp  pa33  +  (mm  +  — 5-  li2X8jpUoupa33  +  Uj,aa33) 


Substituting  (4.43)-(4.46)  in  (4.40)  wc  obtain 

~Y  (“^1  +  +  +  "1  ™  ^  ^K8jaUkJ«3  +  u*aa3)  Y 


+  {m(X!  +  Hj )  +  (l-m)(X2  +  +  {mPi  +  (l-m^u^ 


£i  .  (1-m2)  ^  xC  £i  .  (1-m2)  ,  x 

“  Y  +  m —  ”  T  +  „  M2)(“<*ja3  +  uj.aa3) 


“T  (“^1  +  T  ^2^ajup,pa33 


m* 


1-m3 


~  (m^!  + - —  H2)(8jpUa>po33  +  ^0033)  -  0 

j  m 


(4.46) 


Now  we  introduce  the  following  definitions  in  the  last  equation: 


X12  =  mX|  +  (l^njXj 


4i2  =  mfii  + 


Xj2  -  noXj  +  — — —  X2 

m 


-  1-m2 

M12  =  “Mi  +  — — —  \h. 
m 


^12  =  +  -Mp  ^2 


m 


=  1-m3 

M12  =  m^i  +  — T-^2 
nr 


The  result  will  be 

(X12  +  Ml2)ukjy  +  Ml2ujJI +  Y  {^12uojo3  +  ^12(8jaukjHl3  +  uj.aa3)) 


Y  ^12  ®ajutka3  “  Y  ^12^U°tOct3  +  uj.«*3) 


£2  £2  _ 

Y  ^12^ojuP.Pa33  "  Y  Pl2(8jpua.pa33  +  uj,aa33  =  0 


(4.47) 


(4.48) 


The  above  equations  are  counterparts  of  the  classical  equations  for  linear  elasto-static  problems 
in  the  absence  of  body  forces.  We  need  to  write  these  equations  in  the  expanded  form.  The 
result  would  be  three  scalar  equations  as  follows: 
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9  dui  du2  3u3  d^Ui  9^,  d^i 

(^12  +  Hi2)^(— r  +  -^-  +  -^')  +  lii2(-r-2-  +  -r-2-  +  -r-y) 

CFv  C7U  CnT  Crl T  Owj  (a?2  OO3 


.  $1  -  d2  ,^1  .  *“2*  .  $1  -  d2  >1  ,  3U2  .  *>3, 
+  T  a12  '3a“  *  32T' *  “7"  to  Ha  aa  *35~  +  35T  +  3a' 

Z  C/vjC/V]  C/uj  Ovj  Z  OyjOvj  Ovj  vuj  OO3 


JU-  a  ^ui  ,  a^ii,  $i_r  92  jui  8u2  ,  au3, 
2  412  863 1  aef  ael  2  12  ae^  ^  ae2  ae3} 


li"  a2  ,  ^2,  Si-  a  ,jV  a2^ 

2  ^12  ae^  *ao,  a©2  2  ^12  303 1  aej  ae22 


Ilt  a3  au2  4?=  a3  au,  au2 

3  12  aejdel  ^  a02J  3  *2  ae^2  ^ae,  ae2' 


=  92  a2^  3% 

“Tp“Mr<a5r+aP’)-° 

J  UU3  OUj  OU2 


(4.49) 


This  is  the  first  equation  for  j  =  1  which  can  further  be  simplified  as 

a  3ui  3u2  du3  aV  a^,  a2^ 

a‘2 +Mae7(ae7+de;  +  a%,+ ^ai? + lef  +  ae?’ 


-  -  a3u3 


£?  _  _  a3  du,  3u2 

T(Xl2+Pl2)^F(aer+3^) 


52  =  a2  .a^i  aV 

3  ^12  ae2  ae2  +  90 2 

J  OU3  OUj  uv/j 


(4.50) 


The  second  equation  will  be 
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„  a  dui  du,  dui  d^u?  d2^ 

(Xl2' +l **  *; ‘ ae7 +  a?' t  atf' *' ‘aif +  mT +  *?> 


4i  -  -  83u3 

+T°‘u"x,j)a ^ 


-  -  3ui  du2 

^"^aS^aST**^ 


-  a2  d2u2  d2u2 

^  OO3  Ovj  OUj 


And  the  third  equation  is 


3  dui  3ui  du3  d^ui  d^u3  3^ui 

<X‘2 + M,2)  3% ( M7 +^+aer)+ Ml2(^f +  ^f+W> 

^1  ~  a2  du!  du2  -  a  ^3  ^3 

+  T  ,x'2  <397  +  V  + 1112  3«r  V  +  mF” 


-  a2  ^1  d  ^“3  ^3 

T 11,2  ( aeT  +  W)  +  *7  (3ef  +  lef  ” 


S2  =  a2  ^3  ^3 

T  "“aif<aiPW)"0 


or 


(4.51) 
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d^Uj  d^Uj  d2^ 

(X12  +  H12)  aS"  +  3S“  +  -S") +  +  +  "^T) 

Owj  Ov2  CTO  3 


3  du]  0u2  du3 
5fl“  (  v>  ^  ao  ^  in 

CTO3  CTO  j  CTOj  CTO3 


£1  -  -  32  9ui  9u2 

+  yfti2-M^7(3^+j^) 


-  — E  —(—  +  —1=0 

3  ^12  dfl2  Sfi2  +  a2  '  ° 

J  CTO3  CTOj  »2 


(4.52) 


For  future  reference  we  will  also  calculate  Gj  for  static  problem  in  the  absence  of  body 
forces.  The  preceding  results  are  used  in  conjunction  with  equation  (4.31)  which  has  been 
simplified  for  such  a  case  and  reproduced  before  the  equation  (4.39).  From  (4.45)  and  (4.46)  and 
(4.47)  we  have 


{mc$j  +  (1~^L-  c^}utta  =  ln  SajUtko  +  SuKtja  +  Uj.aa)  (4  53) 


l_m3  —  — 

{mc$p  +  — —  c$p}utp3a  =  Xn  5ajupiPa3  +  P12tyliaja3  +  Ujtaa3)  (4.54) 

m 


Using  (4.36)  we  get 


{mc$,+  (l-mJcgiJu^^  [mXj  +  (l-m)X2]53j5wuu 


+  [mpj  +  (l-mjpj^fij/  +  63/SjkJu  iM  =  X^jU^  +  Mi2(u3j  +  Up)  (4.55) 


c  jjlca )  uk,a3  =  ^i253jSkauk,a3  +  +  53a8jk)uk,a3 


=  ^■I2&3jua,a3  +  Ml2$jau3,a3  (4.56) 

Substituting  (4.53)-(4.56)  in  relation  for  Oj  we  obtain 
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0j  =  ^I2^3juk^  +  Ml2(u3J  +  °p) +  ~2~  I^'l2^3ju<*.a3  +  ^12^jau3,a3] 

2*  ^2®ajukjux  +  ^12^UaJa  +  Uj^)] 

-  y  S?  I^125ajup.pa3  +  Pl2(5jpUopa3  +  Uj.aa3)l  <4-57) 

Writing  down  the  components  of  Oj  separately  we  get 


du3 

de^ 


°i  -  M-i2(  +  35“^ +  T  ^12 


duj 

ae^ 


Si  -  ** 

2 


iir  a  aui  au2  fa3 

de^  "  2  12  30!  ^30!  d02  d03' 


Si  -  a  3u2  -  a2^ 

<5  M-12  vj  '  v\  3n  '  f\  1^12  '  vj  2 

Z  (7Uj  OTj  Ot>2  Z  TO] 


3% 

w> 


if,  _jt_ 

3  12  ae,ae, 


dui 

<ae: 


Si2  =  a2  3u2x  S?  =  a  ,3^1  ,  a2ui x 

3  ^12  30^03  ( 30!  +  302  }  3  ^12  303  {  30!2  302  ' 


dll!  dtl3  Si  —  r- 

=  M-12  (~xT  +  3S“) +  T  ^12  _  ^12) 
(W3  </\/j  z 


d2^  E,  -  -  3  3ui  3u2 

- y  (^12 + ^12) 

C7vj(/\73  Z  (a/j  (a/j  OV2 


_  !l  JL  ^2Ui  .  ^2ui  .= 
3  ae3  c  a©2  aol 


(4.58) 


du2  du3  El  -  r 

°2  =  Hl2(^~  +  ~jxr)  +  T  0^12  “  ^12) 

(A/3  OU2  z 


362883  2  (Xl2^12  de2  (  362  ae2 


(4.59) 


d2^  4?  jb  -  d2  dui  du2 

■  t  t,i2('srf  i^)  ■  t  ^ + ^  ( »: + 3«r) 

5?  a  aa,= 

3  ae3  aef  39? )til2 


£  _  £  _  £3 

°3  *  ^12ukjc +  2M-12u33  +  y  “  y  PufrnJa  +  u3.aa)  “  y  Pl2u3.aa3 


=  Xi2( 


duj 

de7 


dlH 

ay) +  2^i2 

OU3 


^3  Sir  9  ,dui  d«2 

V\  ^  ^  *12  V\  V  V\  ^  V\ 

OU3  Z  (/v/3  wvj  OO 


Si  -  a  /3u1  9u2^  ^1  -  ,4 , 9^3, 
2  ^12  ae3  (ae2  ae2)_  2  ^12( aef  +  ae2 } 


3^3 

39?  + 


.  3ui  3u2  ^  3ua  Ei  _  _  a  du,  du2 

°3  =  ^12(30" +  ) +  (^12 +  2^n)  +  y  (^12 ~ ^12)  t^y  +  ^y ) 

OU 1  002  OO3  Z  OO3  OO  J  OO2 


_  Si-  ^3  .  =  _JL ^3  ^u3 

2  *2  (  39?  3?  }  3  1112  393  (  36?  39? 


(4.60) 


The  constitutive  relations  (4.28)  and  (4.29)  for  tjj  and  Saj  for  a  flat  composite  are  also 
simplified  here  for  a  composite  laminate  whose  micro-structure  is  composed  of  two  isotropic 
layers.  Using  (4.36)  first  we  simplify  different  terms  of  the  expressions  (4.28)  and  (4.29) 


{mc$)  +  (l-mjc^ju^/  =  {mXj  +  (l-m^lSySuUu 

+  {m^ij  +  (l-m)n2}(5fc8j/  +  8j/8j/  +  u^/ 


=  A.^SjjUkjc  +  M-i2(uio  +  uj.i) 


(4.61) 
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1 

1 

where  (30.47)1>2  are  used. 

■ 

1 

2 

{nx^  +  ^  ^  Ca^/lut/ =  ^ifiafikPkJ  +  ^12(Sak®j/ +  $a£jk)uk ./ 

1 

^^AjUkJt  +  MUaj  +  V) 

(4.62) 

1 

1 

2 

{  DSCyka  +  m  cijka )  uk.a3  ~  ^12^ij^kauk.a3  ^  Ml2(^k®ja  +  ®ia^jk)uk,a3 

1 

1 

=  ^125ijUa,a3  +  Ml2(8ja«i,a3  +  8iauj,a3) 

(4.63) 

1 

«&+  m2  c^)uw=X1AA»»W  +  MI2(5ak8#  +  ^u^ 

1 

-  ^12^ajuP,P3  +  Pl2(8jpUa,P3  +  uj.a3) 

(4.64) 

1 

Substituting  (4.61)-(4.64)  in  (4.28)  and  (4.29)  we  obtain 

I 

Tjj  =  XjzSyUtt  +  n12(uy  +  uj4) 

■ 

+  ~Y  [^12^0,03  +  Hi2(Sjaui.a3  +  8kxuj,a3)] 

(4.65) 

Saj  =  Y  [^nSajUtk  +  ?12(Uaj  +  Ujt0)] 

1 

1 

+  J  [^125ajup.p3  +  Pl2(5jpUo,p3  +  uj,a3)l 

(4.66) 

I 

1 

BASE 

5.0  LINEAR  CONSTITUTIVE  RELATIONS  FOR  A  MULTI-CONSTITUENT 
COMPOSITE 

In  this  section  we  assume  that  the  representative  micro-structure  is  composed  of  n  layers 
with  different  constituents.  For  such  a  micro-structure  we  let 

t(aj  =  c$S£  (a=U,...,n)  (5.1) 

where  c$  (a  =  l,...,n)  are  material  constants  in  the  associated  layers.  As  before  the  variable  4  is 
designated  to  change  across  the  micro-structure  whose  thickness  is  assumed  to  be  5„.  It  should 
be  noted  that  although  the  micro-structure  is  composed  of  n  layers.  £n  is  still  supposed  to  be  a 
very  small  number.  The  range  of  variation  of  5  in  the  Ith  layer  of  the  microstructure  is  from  5m 
to  5/  where  /  =  l,...,n  and  5o  =  0-  This  convention  is  adopted  due  to  its  agreement  with  the  spe¬ 
cial  case  of  a  two-layered  micro-structure  which  was  studied  before.  We  further  define  (n-1) 
constants  mlt . . .  .nvi  according  to  the  following  relations 

5/  =  m£n  (/  =  l,...,n-l)  (5.2) 

As  a  result  of  this  definition  the  thickness  of  the  /*  layer  of  the  micro-structure  is  equal  to 
(mf-m^)^  where  /  =  l,...,n  and  =  0,  m„  =  1.  The  composite  stress  -vector  T1  and  the  compo¬ 
site  couple  stress  S®  and  other  quantities  are  obviously  defined  over  the  whole  thickness  of  the 
micro-structure.  For  example  we  have 


(5.3) 


(5.4) 


In  order  to  derive  appropriate  constitutive  relations  we  make  the  following  definition.  Let 
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a  =  a*  for  4m <5 <5/  (/=  l,2,...,n)  (5.5) 

where  a/s  are  constants  and  4o  *  0  as  noted  earlier.  The  function  a  is  piecewise  continuous  for  4 

€  (0, 4n)  and  we  can  evaluate  the  following  integral 

<k*~1)  (56) 

Using  definitions  (5.2)  we  further  simplify  (5.6) 

J  I  1  •Km,k*1-mJiV)  (k*-l)  (5.7) 

where  m*,  =  0  and  m„  =  1.  To  simplify  the  final  results  in  constitutive  relations  we  first  notice 
that  the  integrals  which  appear  in  these  equations  are  the  weighted  averages  of  the  constitutive 
coefficients.  So  we  adopt  the  following  definition 

i  ^ 

!<k)pqr»A  1  j  ^pqr.^  (5 .8) 

Sn  o 

which  by  (5.7)  is  seen  to  be  equal  to 


^licr(m,ktl-m“)  <59) 

We  use  the  same  contravariant  or  covariant  index  notations  for  I  and  c.  However,  the  weighting 
number  k  is  always  written  as  a  superscript  in  parentheses.  Whenever  the  covariant  components 
of  constitutive  coefficients  are  used,  the  layer  index  (/)  is  also  written  as  a  superscript  in 
parentheses.  Recall  (3.22)  which  for  the  present  situation  can  be  written  as 


g*!«  g;  dk + Kkt4~  S  *;<* 

Sn  o  »n  o 


Combining  (3.4)  and  (2.9)3  we  obtain 


(5.10) 
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(‘“gg'-^O  +  ^Xlg  +  SMlj)  (5.11) 

and  (3.5)  reads 

gm  fc  -  i'n  (1  +  ^-)fc  (512) 

The  first  integral  in  relation  (5.10)  is  simplified  using  (5.1 1)  and  (5.12).  We  have 

f  J  c*1  g*w  t*  d4 = f  J  («»"  sln  ti + 0°“  *,,fl  fcM5 

Sn  o  in  o 

=  t'a  eg  £  /  (1  +  f^d!;  +  g>«  X&  i-  /«!  +  ^-JcWdS 
+  g‘«g,^J<l  +  (|V“d4 

-  g1/J  J  0  +  jj*s*(l4  +  ^  -i- 1 4(1  +  -^->^4)  (5.13) 

Now  we  use  definition  (5.8)  to  simplify  (5. 13).  The  following  expression  would  be  the  result 


1“  term  in  (5.10)  =  g1/2  gj{I<0)i^+  £-  +  xja(1)ipk/+  A I™^)]  (5.14) 

2g  2g 

The  second  integral  in  (5.10)  is  also  simplified  similarly 


1  j  ^cij/a  g*i/2  gj*  d£  =  j-  j  %  g*1/2  (c*tagj  +  c0tag3)d4 

7n  o  Sn  o 

■  i ! 5  gW  a + + £  j  gis  a + ^ 


+  -jr  {  4  gm(l  +  •|j)cB*‘fod4 
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-*1®  *](•£•  J  to + |j)c|',x+>4  a  j  + 1^)0^} 

(5.15) 

which  again  by  using  (5.8)  is  simplified  to 

2nd  term  in  (5.10)  =  g1/2  g^1^  +  A  I0***  +  Xil*2^  +  xl  A  I0***)  (5.16) 

J  2g  v  y  2g 

Substituting  (5.14)  and  (5.16)  in  (5.10)  we  get  the  following  constitutive  relation  for  T1 


r  =  { +  A  jOMju + xjaoww  +  A  I®**)]*, 

2g  ~  2g  '  ** 


+  +  A  J(2WI»  +  Xja(2)ipta  +  A  iW^lK^JgW  gj  (5.17) 


The  expression  inside  the  bracket  is  obviously  the  constitutive  relation  for  x1-’.  Hence 


Vi = p®**' + A  ioxh + xtoawu + A  fMUffa 

2g  K  2g 

+  pw + A  i®»*> + jJaa»'“+  A 

2g  H  2g 


(5.18) 


The  same  steps  are  followed  to  derive  the  constitutive  relation  for  the  composite  couple 
stress  S“.  By  (3.27)  we  have 

sa  =  Tk  T-  J  5  g*1/2  cajkig/d^  +  Kjg  A  J  ^2caj/p  g*1/2  g* d£  (5.19) 

Sn  o  w  o 

which  can  be  reduced  to  the  following  form  by  exactly  using  the  same  procedure 

S“  =  (l'lw*t'+  A  ,<2)aiki+^I(2x.1ii+  ± , 


+ p®"1^  +  A  +  A  i<4>ot^)]K/p)  g1/2  gj 

^g  *g 


(5.20) 


BASE 


5-5 


Subsequently  the  constitutive  relation  for  S°^  would  be 


scg  =  pomp  +  A  pw + +  A  rPWyfo 

2g  ~  2g 


+  pOJajfli  +  A  jO*#  +  xiflO)*#  +  A 

o 


4 

2g 


(5.21) 


For  small  deformations  of  a  composite  with  initially  flat  plies,  the  foregoing  equations  are 
further  simplified.  The  resulting  relations  are  recorded  here: 


saj  =  i<1)ajki^  +  i(2)ajiplc^ 


(5.22) 


(5.23) 


with  no  distinction  between  contravariant  and  covariant  components.  In  terms  of  displacement 
vector  u  and  its  gradients  these  equations  can  be  written  as  follows: 


*ij  =  ^K/+I$Xa3 

Saj  =  Jojd/Uk./ +  Iqjku4ft3 


(5.24) 

(5.25) 


Using  (5.9)  the  constitutive  coefficients  are  written  in  the  expanded  form 


I$)  =  CiS&m,  -  m^) 

1$=}^ c^mr  ~  mr-i)  (5.26) 

' =  J  $n2  C$W  -  “r-l) 

To  recapitulate,  d^  (r  =  l,...,n)  are  the  constitutive  coefficients  of  the  micro-structure  layers  and 
m/s  (r  =  l,...,n-l)  are  dimensionless  constants  related  to  the  thicknesses  of  different  layers  with 
iHq  =  0  and  m„  =  1. 
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If  the  micro-structure  is  composed  of  n  isotropic  layers,  we  can  write  (5.26)  in  terms  of  the 
Lame’s  constants  of  various  layers.  In  fact,  we  have 


effir  —  +  p<r)(5jk8j/  +  8y5j^)  r  =  l,...,n  (5.27) 

For  such  a  case,  the  relations  (5.24)  and  (5.25)  are  written  in  expanded  forms  as  follows 


n  n 

xij  =  Uk*$ij  £  ^r)K-nVi)  +  ( W  2^  M<r)(m1-mr_1) 
E_  n 

+  ~2  I  ^tr)(®r2-n£l)  +  (Mja'^kiSj/) 


x  u/.a3  E  MW-mr-i)] 


xij  =  ukjc  S|j  *  \r)  Amx  +  (Ujj+Uj^i)  2:  Z  X^Am2 

£_  n 

+  Y  (ui.a35ja+Uj.a35ia)  £  H^jAm2  (5.28) 

Saj  =  4-  USojSuUu  +  (8ak8j/+ So^jOuti  I  li^yAm2] 

•  I^l  pi 

+  J  ^nt8aj8^U/p3  XfrjAm3  +  (S^p  +  5ap8j/)uiip3  l^rjAm3] 

1  n  .  n 

Soj  =  -Z  USajUk*  \  +  (Uaj+U^a)  Z  ll^Am2] 

“  Pi  i=l 

1  .  n  n 

+  J  U5«jUp,p3  VrW  +  (8jpUa.p3  +  Uj,a3)  S  H^Am3] 


(5.29) 
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where  for  brevity  we  have  introduced 

AmP  =  -  m£i  r=l,...,n  (5.30) 

and  in  the  above  relations  p  =  1,2,3- 

In  order  to  obtain  the  complete  field  equations  for  a  linear  elastic  composite  whose  micro¬ 
structures  comprise  n  layers  we  should  substitute  the  constitutive  relations  (5.24)  and  (5.25)  in 
the  equations  of  motion  (2.123)  and  (2.124).  However  before  that  we  should  obtain  appropriate 
expressions  for  pQ,  p0z1  and  p0z?.  We  assume  that  the  micro-structure  layers  are  homogeneous 
with  densities  (r  =  l,...,n)  in  the  reference  configuration.  Recalling  (2.45)  and  (2.46)  we  can 
write 

pgi*  =  p0Gw  =  j-  J  P;G-'«d4  (5.31) 

Sn  o 

where 

Po=Por)  ft*  £1-1  <\<Ki  (r=l . n) 

and  (5.32) 

^,  =  0 

Using  (4.2),  (5.32)  and  (5.7)  in  (5.31)  we  conclude 

i  .  Pa  n  aL,  n 

P°  =  /  (1  +  f|)Pod$  =  X  Por)Amr  +  X  p^Am,2  (5.33) 

We  will  proceed  similarly  to  calculate  PqZ1  and  ppZ2  using  their  definitions  for  the  present  situa¬ 
tion.  By  (2.42)  we  have 

pg^z1  =  p^z1  =  j-  j  !jp*G*1/2d£  (5.34) 

Sn  o 
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(5.35) 


Similar  to  what  was  done  in  the  derivation  of  (5.33)  we  write 

P°I’  =  t'I5(1+  2G)P;d5=2  5oJ1P?)W  +  |  PoWAm,5  (5.36) 

and 

Po*2  =  +  H>Po‘d!; - i «• f  i  p,«Am?  +-^-|  (5.37) 

For  a  composite  with  initially  flat  plies,  equations  (5.35),  (5.36)  and  (5.37)  are  reduced 
respectively  to 


Po=  I  po^Am, 

(5.38) 

r=l 

PoZ1  =  7  5n  S  Por>Amr2 

l  pel 

(5.39) 

PoZ2  -  y  Sj  pir>Amr3 

(5.40) 

Now  we  can  substitute  (5.24),  (5.25),  (5.38M5.40)  and  (4.23)  in  (2.123)  and  (2.124)  to 
derive  the  linear  equations  of  motion  for  a  flat  composite.  The  resulting  equations  are  recorded 
below: 


lo$/Uk.ta  +  Ia]iKa03  +  bi ,?!  P»r)Am*  +  aP 


(5.41) 


Iojiuk.ta  +  4jtula|i3  +  °j-  I^k./  “  lij&utP3 


+  Cj  X  p “Am,  =  i  yi,  E X  p“Am,3  (5.42) 

i»l  Z  r=l  J  i=l 

The  appropriate  differential  equation  for  the  displacement  vector  u  is  obtained  by  eliminat¬ 
ing  Oj  from  equations  (5.41)  and  (5.42).  For  static  problems  with  no  body  force  these  equations 
reduce  to 


(5.43) 

(5.44) 

Ia$7uk.ta  +  !ajl[puk,aP3  +  +  lijlcPuk.P33  “  ~  Ia?Puk.aP33  =  0 

or 


+  lA&ufca3P  +  °p=  0 

+  I§pUk.ap3  +  <*j  “  tij&Uk./  “  !$pukp3  *  0 
Eliminating  Oj  between  these  equations,  we  get 


Ii>k)uk.fi  +  Iijkjluk4p3  “  Ioj^uk.ta3  ”  ^puk.aP33  =  0  (5-45) 

This  is  a  fourth  order  partial  differential  equation  for  the  displacement  vector  u.  The  constitutive 
coefficients  lfy  (r  =  0,1,2)  ate  already  written  in  expanded  forms  in  Eqs.  (5.26). 

For  a  composite  laminate  whose  micro-structure  is  composed  of  n  isotropic  layers  we  use 
(5.27)  and  (5.26)  to  rewrite  (5.45).  The  result  is 


n  n 

Sy5yUk.fi  J-j  *>(i)Amr  +  (SfcSjj  +  6y5jk)uyfi  ^Z  l^rjAnij 

1  n  ®  , 
+  2  ^nl5ij6kpuk<ip3  I  +  (6*5^  +  6ip5jk)uyip3  2)  ^jAm/] 

1  n  n 

-  y  ^n[5aj5k/UtJa3  Z  X^jAm/  +  (5^5^  +  5CJ6jk)uk<jQ[3  ^Z  p^jAnij  ] 
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4n2[5aj5k0Uta033  *<r)M3  +  +  fiapSpJUkjapB  ^  M<r)M3]  =  0 

1=1  I»1 


or 


uyi(  I  X(r)AmT  +  X  HfrjAm,)  +  X  m^Am, 

r=l  i»l  r=l 

1  n  » 

+  2  ^>^UPJ33  *<r)M2  +  (SjP«UP3  +  Uj.pp3>  Sj  H<r)M2] 

1  n  ,  n  - 

-  J  U&tfVfa*  £  +  (Uaja3  +  uj.aa3)  £  ^rM  l 

1  n  n 

-  J  ^[SajUp.aP33  Z  +  (5jPua,a333  +  Ujjua)  jE  ^jAm,3] 

n  n 

=  “iji  £  (\r)  +  MM  +  ujai  Z  MM 

pi  r=l 

+  J  ^n{(Uaja3  "  ^  l\z)  ~  MM2) 

~  4*  4n2{SajUp,ap33  2  [*«  +  MH3  +  uj,aa33  £  MM3)  =  0  (5-46) 

J  1=1  1=1 

There  are  three  partial  differential  equations  of  fourth  order  for  displacement  vector  u  and  we 
can  write  them  in  two  separate  sets,  one  for  j  =  y  and  the  other  for  j  =  3.  For  j  =  y  we  obtain 


n  n 

“iii  £  (h) +  MM  +  Uroi  J*  MM 


l  n 

-  u3.-y33  ^  $-(r)  “  M<r))Mr2 

1  n  n 

-  T  ^n  (uP.lP33  2  (X<r)  +  MM*  +  »y.aa33  2  M<r)M  )  =  0 

J  1=1  1=1 

and  for  j  =  3  we  get 


(5.47) 
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“M3  2  OHj) +  Ni))^  +  u3.ii  +  T  ^“0.033  2  (\r)  ~ 

1=1  1=1  *•  Is  1 

~  y  ^3,0033  MtryAm,3  =  0  (5.48) 

For  the  special  case  where  n  =  2  these  equations  reduce  to  the  equations  (4.50)-(4.52)  derived 
before  for  a  bi-laminate  composite.  From  (5.44)  we  can  also  calculate  Oj 

°j  =  ~  Uj^kjo  +  l$Puk.P3  "  IojicPuk.oP3  (5.49) 

which  for  the  isotropic  case  reduces  to 


n  n 

Oj  =  (53j5yUkj)  2  ^ry^r  +  (^3 k^jl  +  2 


-  —  [SafikPkJa  J*  ^ry^mr  +  (5^8^+  Sa/Sjk)uk'ja  1  1 

”  2  “  2 
+  ~z~  [53j5kpuk.p3  ^  \xfont  +  (83p5jk  +  53kSjp)u^p3  L  ti<i)Amr  ] 

L  r=l  i=I 

-  y  4n[^oj®kP^c,oP3  ^  \'pssi  +  Kk^jp  +  ^apS^u^  2^  l^rjAm3] 
The  equations  (5.50)  are  written  in  two  separate  sets.  For  j  =  ywe  get 

n  4n  "  2 

<JT  =  (Up  +  u3(T)  I  H^An^  +  —  U3iY3  I  [M<r)  - 
i=l  *  r=l 

-  y  Kerr  Jj  K)  +  M**"?  +  \aa  J-  ^r)^2] 

1  11  n 

-  y  ^Kor^  ^  (k{r)  +  li^Am3  +  11^,003 


(5.50) 


(5.51) 


and  for  j  =  3  we  obtain 
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O3  —  Uir.it  S  X^yAm*  +  2U3J  £ 
i^l  r»l 

”  l  “  2 

-  Y  (UaJa  +  u3>aa)  mr)AmI2  +  Y  “<w*3  ^  ^(r^V 

-  y  ^3.003  ^  (5.52) 

For  the  special  case  of  a  laminate  with  two  layer  micro-structure  (n=2)  these  equations  reduce  to 
(4.58M4.60). 
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6.0  THERMOMECHANICAL  THEORY  OF  COMPOSITE  LAMINATES 


In  order  to  develop  a  thermo- mechanical  theory  for  the  composite  laminates,  we  begin  by 
writing  down  the  local  balance  of  energy  and  the  Clausius-Duhem  inequality  for  the  k* 
representative  micro-structure.  First  we  introduce  the  following  additional  five  quantities  which 
we  associate  with  a  motion  of  the  micro-structure: 

The  specific  internal  energy  e*  = 

The  heat  flux  vector  q*  =  q*(0a,03(k),^,t) 

The  heat  supply  or  heat  absorption  r*  =  r*(0ot,03*k\5»t) 

The  specific  entropy  q*  =  T|*(0'\03(k),4,t)  and. 

The  local  temperature  0*  =  which  is  assumed  to  be  always  positive.  The 

equation  for  the  local  balance  of  energy  —  the  first  law  of  thermodynamics  —  can  be  written  in 
the  following  form 


*  *  *ii.’ *  *V  _  , 

p  r  -  p  e  + 1  tyj  -  q  I*  =  0  (6.1) 

£  at#  *  0 

where  p  is  the  density  of  the  micro-structure,  q  *  and  Yy  are  defined  by 

q  =q^k  .  7ij  =  ygij  (6.21 

and  covariant  differentiation  is  performed  with  rep  sect  to  the  metric  tensor  g*  of  the  micro¬ 
structure.  Recalling  the  relations 


.♦  » 

8i  =V4 

gij=g*-gj*  (6.3) 

tWVV 
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We  can  write 


**  1.*  1/*  •  .  *  *v 

Yy  =  7«y  =  7<v.i'*j  +Vfc) 
Using  (6.3)3  and  the  symmetry  of  x*lK  we  can  write 


(6.4) 


=  y  (vj  •  x  +  vj  •  t* V) 

= i  (v-  •  g f-«v> + v*  •  g'-1^) .  ,*-«v  • »; 


As  for  the  divergence  of  the  heat  flux  vector,  we  have 


(6.5) 


div  q*  =  q*ktk  =  (g*1/2  q\  (6.6) 

8 

Introducing  the  results  (6.5)  and  (6.6)  in  (6.1),  we  can  write  the  local  energy  equation  in  the  fol¬ 
lowing  alternative  form 


P  v  -  P  v + g’-'V*  ■ »•  -  (g*‘Vy = 0  (6.7) 

The  energy  equation  can  also  be  written  in  terms  of  the  Helmholtz  free  energy  function  defined 
by 


y*  =  e*  -  0*ii*  (6.8) 

The  Clausius-Duhem  inequality  as  a  statement  for  second  law  of  thermodynamics  has  the  fol¬ 
lowing  local  form  for  the  representative  micro-structure 

•1/2  *k 

pW  -  pv  +  e,g’-ln(£— 2  0  (6.9) 

U 

By  combining  (6.9)  and  (6.1)  and  using  (6.6)  we  have  the  inequality 


which  in  terms  of  the  Helmholtz  free  energy  y*  defined  in  (6.8)  becomes 


-pV+n'eV^'-i-qX*®  («•»> 

u 

Now  for  elastic  materials  the  constitutive  relations  for  Helmholtz  free  energy,  the  specific 
entropy  and  die  stress  tensor  can  be  expressed  in  the  following  forms 


y*  =  y*(Yy,e*) 

(6.12) 

•  fry* 

tl  =~  “T 

(6.13) 

de 

(6.14) 

where  the  partial  derivative  with  respect  to  the  symmetric  tensor  yj*  is  understood  to  have  the  fol¬ 
lowing  symmetric  form 


i(i£.+i£) 
2  drj  H 


The  constitutive  relation  for  the  heat  flux  vector  has  the  form 


q'MXI.eX)  <6.i5) 

and  the  response  function  Q*k  in  the  light  of  the  Clausius-Duhem  inequality  is  seen  to  be  res¬ 
tricted  by  the  inequality 


(6.16) 


With  the  help  of  (6.13)  and  (6.14)  the  energy  equation  (6.1)  is  reduced  to  the  following 
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4  4  *V  $A*  *  •  A 

pr  -q  Ik-P  M  =0  (6.17) 

where  we  have  used  die  definition  (6.8)  in  order  to  calculate  e*  in  terms  of  y*  and  then  used  the 
relations  (6.13)  and  (6.14)  to  further  simplify  the  energy  equation.  It  should  be  recalled  that  the 
argument  of  different  functions  in  the  energy  equation  (6.17)  is  (0°,93<k),^,t)  and  this  equation  is 
written  for  each  and  every  representative  element  (k  *  l,2,...,n)  which  repeats  itself  in  our  model 
*nd  n  — >  •«.  For  a  bi- laminate  representative  micro-structure  with  thickness  wc  introduce  the 
following  composite  quantities.  These  relations  can  be  generalized  for  a  multi-constituent 
micro-structure  without  any  difficulty  (see  definitions  5.3  and  5.4) 

PS,nr=f  JpWdt 

^2  o 

(6.18) 

pglnr,ij  PV1/2r‘tdt 
0 

8,nq“=j-Jg‘,Vidt 

^2  o 

(6.19) 

glVA=^Jg‘1/2q^^ 

,  & 

pg1/2n(m)i-^  J  pV^irtW  (m  =  0,1,2)  (6.20) 

We  further  assume  that  the  variation  of  temperature  0*  across  the  micro-structure  is  a  linear 
function  of  hence 

e,(e»,e3<k>.t.) = ♦0(e°,e3<kV) + ^1(e».e3<k>.i)  (6.21) 

In  order  to  derive  the  appropriate  form  of  the  energy  equation  for  the  composite  laminate,  first 
we  write  (6.17)  in  the  following  form 

BASE 
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pV  “V  -  (g’V)*  =  pW  (6.22) 

which  after  using  (6.21)  reduces  to 

pV'V  -  u“Vu  -  p  W  (♦„ + «*.)  (6.23) 

Now  divide  (6.23)  by  %2  and  integrate  with  respect  to  £  from  0  to  £2,  die  result  is 

£  J  pWo*  -  j-2  J  (g’>V%« -i/A  (g“ 

-^JpW^+J/pY1^'^  (6.24) 

Each  term  in  the  above  equation  can  be  written  in  terms  of  the  composite  quantities  introduced 
in  (6.18)-(6.20)  except  the  third  term  which  is  the  difference  between  the  values  of  g*1/2q*3 
above  and  below  the  representative  element  divided  by  its  thickness  £2,  namely 

7-  J  4-  (s’1' Y5** = 7-  te,wq*W**'>,t)  -  g,1V(e“,e3<l>,0] 

Now  we  assume  the  existence  of  the  continuous  function  MB'^.t)  which  coincides  with 
q*3(0“,e3(k),t)  at  B3  =  B30^,  and  further  approximate  the  right  side  of  the  above  equation  as  the 
gradient  of  this  function  multiplied  by  g1/2  in  the  B3  direction,  i.e.,  B(g1/2  h)/9B3.  As  a  result, 
(6.24)  can  be  written  as 

pg,/2r  -  (gInq%  -  rp-  <gW  »  =  PgW«Mo  +  Mi)  (6.25) 

In  writing  (6.25)  we  have  also  made  use  of  the  balance  of  mass  equation. 

Next  we  multiply  (6.23)  by  £,  integrate  with  respect  to  £  from  0  to  £2  and  divide  it  by  £2  to 
get 


BASE 


(6.26) 


£iJpY,'V5  «  -  J  te-VXa  -  /  5  £  <*’■  Vwy 

.  ^ 

“  -jr  I  P*s'1V$(lPo+&M 

Using  integration  by  part,  the  third  term  on  the  left-hand  side  of  (6.26)  can  be  written  as 

£  /  5  £  (*•■  Vx 15 = { YV<* = g,n  <M’>  (6.27) 

which  in  writing  the  last  term  we  have  used  (6.19)j  and  the  definition  of  h  given  above.  Using 
this  result  together  with  the  relations  (6.18)2,  (6.19)2  and  (6.20)  we  can  write  (6.26)  in  the  fol¬ 
lowing  form 


pg1/2ri  ~  (g1/2q“).o  "  g^Oi-q3)  =  pg1/2(<Mi  +  hVz)  (628) 

To  determine  the  appropriate  form  of  constraints  on  the  composite  heat  flux  vectors,  first 
we  write  the  Clausius-Duhem  inequality  (6.16)  in  the  following  form 


g-^qXso 

which  by  (6.21)  reduces  to 


t,r-  +  $4>i,o)  +  im  Ai  S  0  (6.29) 

Next  we  divide  (6.29)  by  £2  and  integrate  with  respect  to  %  from  0  to  £2  which  after  using  (6.19) 
can  be  written  as 


gli2q\a  +  g 


m  qf4>i.a+g1/2q3<MO 


or 


qXo+qf^u+q^i^o 


(6.30) 
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which  is  the  appropriate  form  of  Qausius-Duhem  inequality  for  the  elastic  composite  laminates. 

When  the  rate  of  heat  supply  or  absorption  is  zero  (r  =  r1  =  0)  the  energy  equations  far  the 
composite  reduce  to 


<gw  <1%  +  ^3  (*in  W  +  P  gw  ».TU  +  Pi'll)  =  0  (6.31) 

<i'n  qf)fl  +  tln  (M3)  +  pgw  (*0T|,  +  PiTij)  »  0  (6.32) 

Far  small  deformations  of  composites  with  initially  flat  plies  the  energy  equations  (6.31)  and 
(6.32)  further  reduce  to 


q.S+-^3-  +  P(<t>o’lo  +  Mi)  =  0 


(6.33) 


q“a  +  h  -  q3  +  +  <>^2)  =  0 

where  obviously  no  distinction  should  be  made  between  contravariant  and  covariant  components 
of  heat  flux  vectors. 


The  derivation  of  energy  equations  (6.25)-(6.28)  and  the  Clausius-Duhem  inequality  (6.30) 
for  the  composite  laminates  is  not  affected  by  the  number  of  layers  (or  constituents)  in  the 
representative  micro-structure.  The  only  necessary  modification  in  the  case  of  a  multi- 
constituent  composite  is  the  replacement  of  £2  by  in  definitions  (6.18)-(6.20).  Here  ^  is  the 
thickness  of  the  representative  element  which  is  supposed  to  consist  of  n  layers.  Of  course  ^  is 
still  supposed  to  be  a  very  small  number. 

To  recapitulate,  for  a  composite  whose  micro-structure  is  composed  of  n  layers  with  a  total 
thickness  of  4n  we  have  the  following  relations  for  energy  balance  and  the  Clausius-Duhem  ine¬ 
quality 


pr  -  g"1/2t(g1/5q“).o + (gw  WJ = pWoiio + Pi’ll) 


BASE 


6-8 


pr,  +  q3  -  h  -  g"1/2(g1/2<ha).o  *  p(Mi  +  M2) 

qa4>0.a + q“4>i^ + q3$i  *  0 

where  the  composite  quantities  are  defined  as  follows 

pgl/2r  =  y-  J  p*g*1/2r*d^ 

Sn  0 

pg1^  =  _L  J  p’g’Wr*^  d^ 

Sn  e 

pg1,2n(m)=T-  /  pWVMt  (m=o,u> 

Sn  0 

gi/2qi  =  1  j 
Sn  0 

g'V-f  Jg-'V^ 
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(6.34) 


(6.35) 


7.0  CONSTITUTIVE  RELATIONS  FOR  LINEAR  THERMO-ELASTICITY 

For  a  composite  laminate  whose  micro-structure  is  composed  of  n  layers  with  different 
linear  thermo-elastic  constituents,  we  recall  the  following  constitutive  equations  for  the  stress 
tensor  t*‘j,  entropy  q*  and  the  heat  flux  vector  q* 

tS'cSSii-c&S*  (7.1) 

(P  T|  )(a) =  c(a)Yij  +  (P  c)(a>®  (7-2) 

q£)  =  -4)0‘  (7.3) 

where  c(a)/’  c(i)’  c(a)  and  (a  =  l,2,...,n)  are  constants  in  the  associated  layers.  Moreover,  we 
have  the  following  symmetries 


Now  we  proceed  to  calculate  the  appropriate  constitutive  relations  for  composite  stress  vec¬ 
tor  T1,  composite  couple  stress  S“,  composite  entropy  T|(m)  (m  =  0,1,2),  and  composite  heat  flux 
vectors  q1  and  q“.  The  contribution  of  the  first  part  of  (7.1)  to  the  constitutive  relations  for  T1 
and  S°  ( and  consequently  and  Saj)  has  already  been  calculated  (see  section  5).  Therefore  we 
need  to  find  out  the  effect  of  the  second  part  of  (7.1)  in  the  constitutive  relations  for  T1  and  S“. 
Similar  to  what  was  done  in  section  5  we  adopt  the  following  definitions  for  the  weighted  aver¬ 
ages  of  various  quantities 


(7.6) 

(7.7) 


Now  recalling  (5.1 1),  (5.12),  (6.21)  and  (7.6)  we  can  write 


*'■"*;<* 

Sn  o  Sn  0 

+  "r-  J  4c(i>  %m  g*  =  —•  J  (c*gp  +  c°fo)  g*1/2  d£ 

to  o  5n  o 

+j-)Ks'm(c%+c*t;xik 

TO  0 

= ♦„sin  eP  J  (i  +  + ♦„  gw  /  4<  i  +  ^*^4 

4b  .  1/2  4b 

+  *.  gW  *3  ^  /  (1  +  gg  J  5d  + 

+ ♦, im +  J 4(1  +  ^K°d4 

=  Ko  g,n  g,(-^  J  (1  +  I|)cijd4  +  xji|4(l  +  ||)c»d4) 

+ *1  t'n  gj(  ^  J  4(1  +  |j)cljd4 + xj  -i-  J  42(1  + 

=  ♦„  im  gjl  j<°*  +  A  + ^(JW  +  A  J0*1")) 

+ 0,  gln  + A  jWj + xjo®* +  A  jW,, 
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By  combining  the  results  (7.9)  and  (5.17)  we  can  obtain  the  response  function  for  the  com¬ 
posite  stress  vector  T‘.  As  before  if  we  disregard  the  factor  g1^  gj,  what  remains  is  the  constitu¬ 
tive  relation  for  which  is  recorded  below 


t(i,(I«»»'+AI(iWU  +  xJaci»u  +  A 

©  “'O 

+  JlOWto  +  A  lC*jto+XJ(l<2)4U«+  A 

-♦.IJ®*'+Aj<iw+ + A  jc»)) 

2g  K  2g 

-  4>!  (J<1)ij  +  jOW  +  xi(J<2^  +  A  J^P)}  (7.10) 

2g  H  2g 

Similar  steps  are  followed  to  find  the  constitutive  equations  for  S°  and  S°j.  The  contiibution  of 
die  thermal  term  is 

9*  g*1'2  *;<15 = g>«  tj  (J«w +Ajow+ xfei®* +A.jO)»i>)) 


+  tog“M^W  +  ^,W  +  ^W+4-J<4)<‘’))  (7.11) 

Again  combining  (7.11)  and  (5.20)  we  find  the  constitutive  relation  for  S°.  Dropping  the  com¬ 
mon  factor  g1/2  gj  would  result  in  the  constitutive  relation  for  S°^  which  is  recorded  below. 


S“i = TUI0**"  +  A  + xtfP**  +  A 

2g  2g 

+  +  ~  +  A  I<4)°n*P)} 


_  ^{jOW  +  A  J<2)aj  +  jjgOX#  +  A  pw)} 

2g  K  2g 
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-  ♦1{j<2)oj + 4~  p** + + 4~ 

2g  F  2g 


a.12) 


The  next  step  is  to  find  the  constitutive  relation  for  the  specific  entropies  (m  =  0,1,2). 
We  calculate  the  contribution  of  each  term  of  die  relation  (7.2)  separately.  By  (6.35)3  we  have 


pA^4  JpVg'^-dt 

Sn  o 

=  T-  J  «*weMt  *4-t  »Vi" “H-at 

Sn  o  Sn 

Since  g*1/2  =  g1/2  (1  +  4^-)  we  can  write  the  first  part  of  (7.13)  as 

2g 


(7. 13) 


j-  J  ‘(Vi!  g',/2  5”d5  =  ^  J  5m0  + 

Sn  o  Sn  o 


which  by  (2.31)-(2.34)  reduces  to 


-  [  /  t"a + IrWivsj +[^-J^"Hl(i+ ^>cg)d5]^ 


(7.14) 


=  g1/2  { Yy[J(m)ij  +  ~  J<m+1)ij]  +  K-p[J(m+1)iP  +  Aj(n»+2)ip]  j 
The  second  part  of  (7.13),  by  (6.21)  and  (7.7),  is  written  as 


(7.15) 


•  ^  1/2  t  * 

T*  J  (P*c)(«,9'  g-lfl  5”d5  =  -t-  /  5”(1  +  ^)»„+$4>iKP*cW« 


+ ^ir !  ’r'a  +  27)<p‘c)<“^ 
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=  *0gm  {K(m) 


4  A  K^1*}  4  gm  4  A  K(m+2)}  (7.16) 


Substituting  (7-15)  and  (7.16)  in  (7.13)  we  obtain 


pg1/2n(m)  =  g1/2  {Ylj(J<rn)ij  +  A  J^1*)}  +  gm  iqa^1^  +  4-  J°n+2)ip) 

2g  K  2g 

4  <J>0  g1/2  (K(m)  +  A  K(nHl>)  +  +!  g1/2  (K(ro+1)  +  A  K(m+2>) 
2g  2g 


pTl(m)  =  7ij(J<m)ij  4  A  J<m+1)ij)  4  Kfc(jfr»  W  4  A 
2g  2g 

4  <|>0(K(m>  4  A  K(m+1))  4  ^(K(mfl)  4  A  K(m+2>)  (m  =  0,1.2) 
2g  2g 


(7.17) 


Wc  start  with  (6.35)4  to  find  the  constitutive  relation  for  q\  Substituting  from  (7.3)  in  (6.35)4  we 


«“V = f  S  8-,Vd4 =-f  jY‘«  k&e-dS 

«i  o  Sn  o 


(7.18) 


By  (6.21) 


e,J  =  ^o,|J  +  ^l,p 


(7.19) 


®3  =  ^1 


(7.20) 


BASE 


7-6 


Noting  these  results  we  simplify  (7.18).  Hence 


*1/2  q* = -  1  (1  +  !j)ki"(*o.|l + 

-£jo+a*%* 

= -  gin  J  (i  +  + ^6-  /  sa  + 

+il(1+f)kBdU 


(7.21) 


Using  (7.8)  we  get 


gm  q4  =  -  gm  (6o.plLW  +  ~  LW]  +  4>i.p[L(1)*  +  A 


-Hj>,[L(0)i3  +  v-L(1)0]} 
2g 


qi  =  _  (L(0)iP  +  A-  L(1»fy|>0,p  -  (L(1)iP  +  -A.  L®^1>p 
-(L^  +  AlO^ 


(722) 


Finally  we  use  (6.35)5  to  derive  the  constitutive  equation  for  q“.  Similar  to  the  above  develop¬ 


ment,  we  write 


gln  tf  -f  J  g'“ VS  d$  =  -  f  J  5  g''0  k“^d5 


J  ^l  +  ^Kk^eJ+k^ej^ 


BASE 


7-7 


which  by  (7.8)  reduces  to 


(7.23) 


q«  =  -(L(1)QP  +  A  L«2)°P»0.p  -  (L*2**  +  A  L0^1>p  -  (L(1)a3  +  A  (7.24) 

This  concludes  our  derivation  of  linear  thermo-elastic  constitutive  relations  for  composite  lam¬ 
inates. 


For  small  deformations  of  a  composite  with  initially  flat  plies  the  foregoing  equations  are 


simplified  to  the  following  constitutive  relations: 

025) 

^aj =  ^aj^ut{13  —  ^aj^o  —  7aj*9l  (726) 

P’lta)  =  JfVj  +  7i^M'1)“i.p3  +  K<m)»o  +  K<-%  (n>  =  0,1,2)  (727) 

ql  =  -4V6-l41)*l,S-4%  (728) 

q“  =  -  LO^oj,  -  -  L«^,  (729) 


The  constitutive  coefficients  I*0),  I*1*  and  I*2)  have  already  been  calculated  and  recorded  in  equa¬ 
tions  (5.26).  As  for  the  other  constitutive  coefficients  we  use  the  results  of  section  5.  Compar¬ 
ing  definitions  (7.6)-(7.8)  with  (5.8)  and  using  the  results  (5.4)  and  (5.30)  we  can  write 


(7.30) 

K<k>  =  fir  ^  1  (p'c>«Am.k*1 

(7.31) 

(7.32) 

BASE 

It  should  be  noted  that  in  equation  (7.32)  are  the  coefficients  of  thermal  conductivity  of  dif¬ 
ferent  layers  of  the  representative  micro-structure  and  are  not  to  be  confused  with  the  superscript 
k  which  assumes  non-negative  integer  values. 


If  the  micro-structure  is  composed  of  isotropic  layers,  die  coefficients  of  thermal  stress  c$ 
and  thermal  conductivity  can  be  written  in  terms  of  only  one  constant  for  each  layer.  For 
such  case  we  write 

=  07-33) 

kf-Mij  (7-34) 

Taking  note  of  these  relations  and  relations  (7.30)  and  (7.32)  we  obtain 

If  =  SljJifS(I>Anv 

If  Mm,3  (7-35) 

Lf  =  Sy  !  k^Am, 

(7.36) 

Lf  = 

Consequently  the  constitutive  equations  (7.25)- (7 .29)  reduce 
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a 

Ti;  =  5:;  £  (Xrf»\Uh  v  ■ 


+  2  Si  ~  P<r)$i)Ainr2 

4n  a 

+  Y  (uM*3Sjo  +  “ia35ta)  2 

“  T5al 


1  n  n 

Soj  =  y  ^n^oj  J ;(*fr)uk*  ~  Pfr^o)^2  +  (“oj  +  uj,o)  £  p^jAm2} 
^  r=l 

+  "J  ^(8aj  J.  &<r)ufcp3  ~  PfrjMAm3 
a 

+  (Sjpu^pj  +  Ujta3  I  p(r)Amr3} 


n 

P^(°) =  ^  iPw^i  +  (p^^JAm, 


+ 1  ^  li  [P«UW3  +  (P^V^ilAm2 


PH(1)  =  [P(r)Ui4  +  (pV^JAm2 

+  J  ^WUP*P3  +  (P’c^ilAm3 

P%) =  J  &  Zj  tP(r)«i4  +  (p*c)w(|) JAm3 

+  7  &  |,  IP(DUP.P3  +  (P*c)(l)<h]Am^ 


(7.37) 


(7-38) 


(7.39) 


(7.40) 


(7.41) 
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®  i  a  , 

Qa  =  "  $o.a  ^(oAnij  -  —  >a  2^  k<T)Anv 

(7.42) 

n 

q3  =  ~  <Pl  2 

i*l 

(7.43) 

qf  — j  U.  |  M®.*  -  y  |  WW 

(7.44) 

In  relation  (7.44),  a  is  written  as  a  superscript  only  for  convenience  and  does  not  signify  the  con- 
travariant  position. 

The  linear  equations  of  motion  and  balance  of  energy  for  a  composite  with  initially  flat 
plies  are  derived  by  substituting  (7.25)-(7.29)  in  (2.123),  (2.124)  and  (6.33).  Using  the  results 
(5.41)  and  (5.42)  in  conjunction  with  (7.25)  and  (7.26)  we  have  the  following  equations  of 
motion  in  linear  thermo-elastic  theory 

Ia£/Uk.ta  +  *ajfeuAa03  ”  JqjVa  "  +  bj  ^  Po^  +  °j3 

=  iij £  Pir)Amr  P^Am?  a.45) 

- Iijiu/JB  +  jf  ♦.  +  +  Cj  E  p«Am, 


=  7  ^nUj  jc  Po>Ani?  +  J  ^iij3  2^  p^Am3  (7.46) 

The  energy  equations  when  the  rate  of  heat  supply  or  absorption  is  zero  are  recorded  in  relations 
(6.33)  for  small  deformations  of  thermo-elastic  composites  with  initially  fiat  plies.  Substituting 
the  constitutive  relations  (7.27)-(7.29)  in  (6.33)  we  find  the  following  coupled  differential  equa¬ 
tions  for  displacement  and  temperature  fields 


BASE 
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<t40)  +  Hi  +  WH'  *  +  (♦«Km  +  ♦,K<,,)+0 

+  (Wt,,) + ♦iK°>)i,  +  J|-  LSSto#  -  iat*M  -  LS>*,a  -  0  (7.47) 

(♦4° + ♦ilfKj  +  (♦<$>  +  ♦iJfl’K*  +  (♦^0)  +  ♦l*®#. 

+  (+^C® + ♦1K<3>)4,1  +  h  +  L^p  +  L|J>fc  #  +  Lg^, 

-  L^W  -  L$4l.a(>  -  =  0  (7-48) 

For  static  problems  in  the  absence  of  body  force  and  heat  supply,  the  foregoing  equations 
are  further  reduced  to 


ii'K*  +  I^p*  -  l«j\a  -  +  Oy  =  0  (7-49) 

Iaj^uk.ta  +  lajfeu4aP3  “  JajV.a  ”  loj^U  +  0j”  IJjEluU 

-  I$»4|I3  +  jf ♦„ + Jf -  0  (7.50) 

|p-  -  L^9o.ap  -  -  LS't’U  =  0  (7.51) 

h  +  L$V|> +  MJ’fl.f  +  1^08  -  “  0  (7-57) 


Similar  to  what  was  done  previously  in  order  to  find  a  relation  between  the  director  dis¬ 
placement  and  die  gradient  of  displacement  vector,  we  enforce  the  continuity  of  the  temperature 
field  across  two  adjacent  micro-structures  to  derive  an  analogous  relation  between  <j>0  and  <j>] 
defined  in  equation  (6.21).  In  order  that  the  temperature  field  be  continuous  on  the  common  sur¬ 
face  between  k*  and  (k+l)“  micro-structures  we  should  have 

©•(e^^.O.t)  =  eV^tt)  (7.53) 

Now  by  (6.21)  we  have 
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©•(©“.©^ao = ^(e0,©30*1^) 

©•(©“.e300,^) = ♦o(0°.03<k>.t) + ^(©“.e300.*) 

Substituting  from  (7.54)  and  (7.55)  in  (7.53)  we  get 

^(e0,©30*1^)  *  ^(©“.e300.!) + ^(©“.©^.t) 

or 

♦!(©“, ©^V)  =  j-  {^(©“.©^.t)  -  ^(fl0, e30^)}  (7.56) 

By  smoothing  assumptions  and  noting  the  smallness  of  4n  we  approximate  the  right-hand  side  of 


(7.56)  as  the  gradient  of  <j>0  in  the  0*  direction.  So  we  obtain 

♦i*©”,©^)  =  6o(0a.©3.t)  (7.57) 

This  conclusion  is  used  in  various  field  equations.  In  particular,  equation  (7.49)-(7.52)  reduce  to 

Ioj£/uk.ta  +  ^o]^utap3  “  Jaj^O.a  “  Jap0O,a3  +  °j3  =  ©  (7.58) 

Iqjll/Uk.Ja  +  I«jfeuJ.a03  “  ^oj^a  ~  Jqfto.ctf  +  aj~ 

- 1^4)3  +  jf Vo +  J|j’*oj  =  0  (7.59) 

hj  -  L^o*.  "  -  L$Vo3  -  0  (7.60) 

■>  +  +  M3W3  -  La^9o.ap  -  ^♦o.ap3  “  La3Wa3  -  0  (7-6* ) 


Eliminating  Oj  between  (7.58)  and  (7.59),  and  h  between  (7.60)  and  (7.61)  we  find  the  following 
coupled  differential  equations  for  displacement  and  temperature  fields.  Since  we  are  investigat¬ 
ing  the  static  problems  in  the  present  derivation  the  equation  for  temperature,  i.e.,  die  equation 
resulting  from  the  energy  equations  is  independent  of  the  displacement  field.  Recalling  (5.45), 


(7.54) 

(7.55) 


BASE 


7-13 


the  displacement  equation  becomes 

I$s)uk.fi  +  I$&«k4P3  ~  lSiUk.ta3  “ 

-  -  ■»£' Va  +  MS  Ws  +  Mj  Wj3  =  0  (7.62) 

and  the  temperature  equation  by  (7.60)  and  (7.61)  is 

LfVi,  +  (Lji>  - L®##*,,  -  I^o.<*33  =  0  (7«) 

Having  determined  the  displacement  and  the  temperature  fields,  the  interlaminar  stresses  Oj  and 
heat  flux  h  can  be  determined  from  (7.59)  and  (7.61),  respectively.  The  results  are 

°j =  Jjjk^kjc  +  iJ&W#  “  _  !q$jutap3 

+  OM) 

h  “  ^'(’O.aP  +  +  LfiVid  “  Ma$0,a  “  “  Lgfyoj 

+  (765) 
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8.0  LINEAR  THEORY  OF  INITIALLY  CYLINDRICAL  LAMINATES 

For  an  initially  cylindrical  laminate,  we  choose  the  usual  cylindrical  coordinates  (r,0,z) 
which  are  related  to  the  coordinates  (Q1,©2,©3)  according  to  the  following  relations 

e^e  ,  e2=z  ,  e3=r  (8.i) 

The  choice  of  coordinate  r  for  O3  direction  is  natural,  since  the  cylindrical  laminates  are  piled  up 
in  r  direction.  Hie  metric  tensor  Gy  and  its  conjugate  G1^  and  G  *  det(Gy)  are  given  by  the  fol¬ 
lowing  relations 


V  o  o 

l/r2  0  O' 

0  1  0 

0  1  0 

rH 

o 

o 

0  0  1 

G  =  det(Gy)  =  i2 


(8.2) 


(8.3) 


Since  we  are  working  within  the  realm  of  linear  theory,  the  reference  metric  tensor  Gy  is 
appropriate  for  calculating  the  covariant  derivatives  of  various  quantities.  We  further  recall  if  vj 
and  Ty  are  covariant  components  of  arbitrary  vectors  and  tensors,  their  physical  components 
and  are  given  by 


V<i>  = 


(no  sum) 


and  in  terms  of  covariant  components 


V<i>  =  ViVG^  ,  T<ij>  =  TijV^VGjJ(nosum) 


(8.4) 


(8.5) 


For  the  cylindrical  coordinates  defined  in  (8.1),  the  non-vanishing  Christoffel  symbols  of  first 
and  second  kind  are  as  follows: 


[13,1]  =  [31,1]  *r  ,  [1 1,3]  =  -r 


(8.6) 
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{l13>  =  b1l>SB7  ’  {l3l}=_r  (8-7) 

For  the  subsequent  analysis  we  also  need  an  expression  for  the  quantity  A  defined  in  (2.14): 


‘2r  0  O' 

V  o  o‘ 

A  = 

0  1  0 

+ 

0  0  0 

0  r  1 

0  0  1, 

(8.8) 


This  quantity  relates  the  determinants  of  the  metric  tensors  of  micro-  and  macro-structures 
through  relation  (2.13). 


8.1  Relative  Kinematic  Measures 


The  relative  kinematic  measures  y-j  and  were  given  by  (2.118)  which  for  cylindrical 
coordinates  since  Da  =  =  0,  by  (2.103)  and  (8.2),  are  simplified  to 


Y.j=y(uiij  +  Uj,i)  (8.9) 

iCap  =  A|juj|0  +  80ip  ,  ^  =  53,0  (8.10) 

In  writing  the  above  relations,  we  have  also  used  the  resuts  (4.22)  and  (2.30)  while  noting  in 
cylindrical  coordinates  we  can  write 


Aj»(33a>.0 

By  straightforward  calculations,  covariant  derivatives  of  displacement  vector  u  and  director  dis¬ 
placement  8  are  found  and  substituted  in  (8.9)  and  (8.10)  in  order  to  obtain  the  covariant  com¬ 
ponents  of  the  relative  kinematic  measures.  These  results  together  with  physical  components  of 
each  tensor  are  recorded  below: 


du] 

Yll  =  le~+ru3 


1  ^  Uj 
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y  _  i  ,3°!  ,  v  _  i  ,  i  du, 

Yl2‘  2  (  01  30  J  ’  *•"  2  (  3z  r  30  ^ 


ui  „  i  due  t  du,  ue 
Yir2(3r  ae)_r  ’  V2(3r  r  ae‘r) 


du2  du, 

•  *-sr 


1  du2  du3  1  du,  du. 

Y23=I("^r  +  lT)  *  Yra='2  (1T  +  1T) 


du3 


dUf 


‘ft3=  dr  *  Yn  =  dr 


”=i(^)+i(ra3)  •  ^=7i(r^)+H(rUr) 


1  d  ,  ^uev  .  1  d 


d  y^ui  «i 


d2!!® 


Kl1  3z  (  dr  T)  '  ^Mt 


1  duj  d2^ 


1  ,du®  d*u. 


*21“  r  d*  +  dGdi  ’  **"  r  (  dz  +  d0^) 


d^2  d2^ 

<22  dzdr  ’  *“  dzdr 


3^3  Uj  *  d\  due 

K31_  dedT“T  '  K'®=7(dedT“ir) 


d2^ 


d2^ 


K32=  dzdr  ’  Kra~  dzdr 


(8.12) 


(8.13) 


The  equation  for  balance  of  mass  is  also  obtained  by  using  (2.121)  and  (2.122).  The  result  is  as 
follows 
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1  d“e  du,  du,  u. 
P  Po  r  80  dz  dr  r 


(8.14) 


82  Linearized  field  equations 

These  equations  have  already  been  derived  for  a  general  composite  laminate  and  are 
recorded  in  (2.123)-(2.125).  In  cylindrical  coordinates  and  in  terms  of  physical  components  of 
different  tensors,  these  equations  are  reduced  to  the  following  forms 


1  .  d**  1  309,,  ^  ,  ijj  , 

TIT  T  TIT  Pob,=Po01*  ^ 


and 


and 


i  dt^  8t„  i  do,  ... 

7 17  +  17  +  7  ir+Pobz=Po(Uz+z  8z) 

1  dt*  diB  “'ee  ,  1  ^  ...  ft . 

7  le- + 17  “  T" + 7  + p A = p»(Uz + z  ^ 


1  dSee  BSjq  Se t  oe  ... 

T'ae-  "aT  ~  p°Ce=Po<z“®  ^ 

7  ■&" + IT + T  '  ^ + PA = p»(z  + ^ 
715"  "ai  7  T“%  p<’c,=Po<z"’  ^ 


Tae  ~  =  “  sez 


(8.15) 


(8.16) 


(8.17) 
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tn  =  *«r 

As  noted  earlier  in  (1.125)  the  composite  stress  is  not  in  general  symmetric  and  in  particular 
in  cylindrical  coordinates  the  asymmetry  is  represented  by  equations  (8.17). 

8.3  Constitutive  relations 

We  assume  that  die  representative  micro-structure  of  the  initially  cylindrical  composite  is 
composed  of  n  layers  with  different  linear  constitution.  The  results  developed  in  section  5  are 
used  to  derive  the  appropriate  constitutive  relations.  For  the  present  case  it  is  noted  that 


A  _  1 
2G  r 


and  the  only  non-vanishing  component  of  AjJ  by  noting  (8.7)  is 


(8.18) 


A{  «  l3 1  !>  =  7  (8.19) 

Using  these  results,  we  simplify  (5.18)  and  (5.21)  and  derive  the  following  contravariant  forms 
for  the  constitutive  relations  in  cylindrical  coordinates 


fl  _  (]<0)ilk/  +  2  j(l)ilk/  +  +  qdHWo  +  1  j(2)iUa  + 

t2 = 0®““  +  j  + 1  j?v*)Kia 

t°  =  (pwu  +  1  +  (I<l)i3/a  +  1 

Sal  _  (j<l)alk/  +  2  j(2)olk/  +  ^  +  (jtfjal/p  +  2  j<3)al/p  +  ±_  j(4)al/p)|C/p 

so2  .  tfXVM  +  1  +  q(2)a2tp  +  1 


(8.20) 

(8.21) 

(8.22) 

(8.23) 

(8.24) 
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S03  =  (I<1)a3w  + 1 (8.25) 

where  the  covariant  components  of  the  relative  kinematic  measures  Yjj  and  tc^  have  already  been 
recorded  in  (8.12).  The  composite  constitutive  coefficients  I<k)PT*  depend  on  the  material 
behavior  of  every  constituent  of  the  composite  laminte.  For  die  special  case  where  die  micro¬ 
structure  is  composed  of  n  isotropic  layers  with  different  elastic  constants,  we  can  write 

qff*  =  A^GWG”  +  +  G^G*)  (8.26) 

where  and  (/ «  l,...,n)  are  the  Lame's  constants  of  each  layer  in  the  micro-structure.  By 
(8.2)  the  non-vanishing  constitutive  coefficients  are  as  follows: 


C1111  *  -y  (X  +  2p) 
r 


rti33  _ 

*”  r2 

C2222  =  X  +  2p 

C2233  =  X  (8.27) 

C3333  =  X  +  2p 

C2323  =  p 


ri3i3_  JL 
'r2 


ri2i2_iL 

'r2 


where  subscript  (/)  is  dropped  for  simplicity.  It  should  be  recalled  that  following  symmetries  of 
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the  constitutive  coefficients  must  be  remembered  when  the  expressions  involving  the  constitu¬ 
tive  coefficients  are  expanded 


(8.28) 


By  (5.9)  the  non-zero  components  of  the  composite  constitutive  coefficients  have  the  same 
superscript  as  (8.27).  Substituting  from  (8.27)  in  (5.9)  and  using  the  results  in  (8.20H8.25)  we 
find  the  following  constitutive  relations  for  the  physical  components  of  the  composite  stresses 
and  composite  stress  couples.  The  summations  in  these  relations  extend  over  the  micro-structure 
from  /  *  1  to  /  *  n 


L  £2 

*tee  =  YeeZ  (*-  +  ^(Am,  +  -y-  Am,2  +  Am/3) 

p  £2 

+  (Yzz  +  Tn)Z  +  —  Am/2  +  Am/3) 


+  *«  Z  +  24;)(  y  Am,5  +  %  W  +  ^  Am,4) 

+  K=ZMYAm/5+-^Am,3+-^-4m ,4)  (8.29) 

g  £2 

XaB  =  2Yz9  Z  I^KAm/  +  —  Am,2  +  Am,3)  +  (k*  (8.30) 

+  « Y  Am,’  +  -^  Am,3  +  -^  Am,4) 

^ie  =  2Yi«Z4/(Am/+y  Am,2+-~Am,3)  (8.31) 

+  Z  M/' Y  Am,3 +  Am,3 +  £  Am,4) 
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P  P  ?2 

^z  =  2YteX  H/CAm,  +  —■  W)  +  (Krf  +  I MK  y  W  +  Am/)  (8.32) 

1,  L, 

=  (Yee+Yn>Z  VAn^f  W>  +  Y«  Z  0*  +  2n/XAm/  +  ^  Am/)  (8.33) 

+  Kee  I M  y  ^  + 17  W)  +  k„  £  (X, + 2^X  y  W  +  y  ^ 


P  P  P* 

V  =  2y«  I  UKAm,  +  -g-  Am,2)  +  K„  £  M  Y  +  Y  Am,3)  (8.34) 

p  P  1} 

%  =  2Yer  Z  )i/(Am/  +  y  An^  +  Kr«Z  UK  Y  Am'2  +  37  (8,35) 

F  F  p2 

^  =  2YzrI  Mi(Am/  +  —  Am/3)  +  K*  £  ^Y  W  +  ^  *“/*>  (8.36) 

V  *  (Yee+Yiz+Yu)  Z  +  —■  Am/2)  +  2y„  £  M/(Am/  +  ~  Am/)  (8.37) 


F  Pf 

+  (K0e+K22)  Z  M  y  W + 17  W) 

See  =  $n  Yee  Z  (*■/  +  Wy  Am/  +  ^  Am,3  +  ^  Am/)  (8.38) 

+  5n(Yz7+Yrr)Z  *Ky  Am,2  +  Amz3  +  ^2  Am/4) 

+  ^ee  Z  (V2M/)(  y  Am/  +  ^  ^  Am/  +  -^  Am/) 

+  Z  ^/  ( J  Am/3  +  ~  ^  Am/4  +  -p-  ^2  Am/5) 

Sje  =  2^  2  H/(y  Am/2  +  Am/  +  yy  Am/)  +  ^(k^  (8.39) 
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+  Am,®  + -J- Am,4  + -p- Am,5) 

S«z  *  21b,  £  nKy  W  +  —•  Am,3)  +  £?(*«*  (840) 

+  «i»>Zt*<<yAmj,  +  -|-Am,4) 

SD  =  4j=  I  ft,  +  2H,)(|  W  +  ^  Am,3)  (8.41) 

+  4?  X  0,  +  Am,3  +  ^  Am,4) 

Sfc  =  ^  Ife  2  (IK  Y  W  +  W)  +  ?.2  «i«  I MK J  Am,3  +  -  J  Am,4)  (8.42) 

S„  =  24.  T„  Z  Am,2  +  ^-  to,3)  +  ?,J  2  (JKj  Am,3  +  ^-  Am,4)  (8.43) 

Using  (5.33),  (5.36)  and  (5.37),  die  composite  mass  density  p0  and  the  composite  mass 
moments  PqZ1  and  p0z2  which  appear  in  the  equations  of  motion  are  also  calculated  for  an  ini¬ 
tially  cylindrical  laminate 

Po  =  I  Po(/)Ani/ +  ^  S  Po(/)Am2  (8.44) 

1  E2 

PoZ1  -  j  Sn  1  PofljAm,2  +  jp  x  Po(oAm?  (8.45) 

1  t3 

pc^2  =  y  42  X  PcftAm,3  +  2  Po®Amf  (8.46) 

BASE 
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8.4  Energy  equations  and  constitutive  relations  in  linear  thennodastidty 

Energy  equations  (6.34)^  in  the  absence  of  heat  supply  or  heat  absorption  reduce  to  the 
following  forms  when  written  in  the  cylindrical  coordinates  (8.1) 

7^  +  ir  +  T^+p°(^+*lT,,)'0  (fU7) 

7^  +  ^i+i,-q,+p<,<Mi+»,%)=°  (S.48> 

where  q^q^q,  are  the  physical  components  of  heat  flux  vector  for  the  component  laminate; 
qie,qiz  are  the  physical  components  of  the  composite  heat  flux  moment,  and  other  quantities 
have  the  same  meaning  as  section  6.  The  constitutive  relations  for  various  composite  quantities 
in  thermoelasticity  can  be  written  with  reference  to  the  development  in  section  7.  The  mechani¬ 
cal  parts  of  the  constitutive  equations  for  composite  stress  and  stress  moments  were  derived  in 
section  8.3.  Therefore  in  what  follows  we  record  the  thermal  contributions  to  these  quantities. 

The  complete  constitutive  relations  in  thermo- mechanical  theory  is  obviously  obtained  by  super¬ 
imposing  these  distinct  parts. 

From  (7.10)  we  have 

^(thermal)  =  +  i.  jO»  +  +  i  J®»))  (g.49) 

-  ( jOW  +  I  jGW  +  (jWf  +  I 

Using  (8.19)  this  reduces  to 

1“  =  +  -Lj<Wi|_ ♦,{  J*1*1  +  -  J®»  +  4  J®*1 )  <8.50) 

r  r  r  r 

f2  =  -♦o[I<0)£!  +  -1-  +  -L  J«»]  (8.51) 
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t*  -  -+0[fl>»  +  }  J0*3)  -  ♦it*1*3  +  }  (8.52) 

where  the  word  thermal  is  dropped  for  brevity.  Similarly  from  (7.12)  we  find  the  thennal  contri¬ 
butions  to  stress  moments  as  follows: 

S»>  =  -♦„{  J<‘»'  +  2  + 1.  jom )  -  4,{  J<J»1  +  2  J0»>  +  ^  J<4»» )  (833) 

s02 = -♦.(J<1W + 2  j®02]  -  + 2  J®"2] 

S"3  =  +  2  J®o3]  _  +  I jO*3] 

r  r 

The  constitutive  relations  for  composite  entropy  and  its  moments  qm,  heat  flux  vector  q  and  its 
moment  qx  are  also  derived  from  equations  (7.17),  (7.22)  and  (7.24).  Using  (8.18)  we  have 

p0Tim  *  +  k^"*^  +  j  J<m+2)iP] 

+  <|>0[K(m)  + 1  K(ra+1>]  +  ^>j[K(m+1)  +  j  (m  =  0,1,2)  (8.56) 

q‘  =  -  [Lw  +  j  L(1)*]4>0>p  -  [L(1)ip  +  j  -  [L(0)i3  +  j  L(1>3]<|>,  (8.57) 

qo  =  ^(DoP  + 1  -  [L<2*#  +  j 

-  (L(1)a3  +  y  L*2**3]^  (8.58) 

where  the  thennal  constitutive  coefficients  and  for  the  composite  laminate  are 

given  in  (7.6)-(7.2)  and  calculated  in  expanded  form  in  (7.30)-(7.32). 

If  the  micro-structure  is  composed  of  isotropic  layers  the  coefficients  of  thermal  stress  and 
thermal  conductivity  of  each  layer  can  be  represented  in  terms  of  only  one  constant.  For  the 

BASE 


(8.54) 

(8.55) 


present  case  of  cylindrical  laminates  we  can  write 


C$=G%  (/=  1,...*) 

(8.59) 

(8.60) 

Substituting  from  (8.59)  and  (8.60)  in  (7.30)  and  (7.32)  we  have 

(8.61) 

J0022  _  j<k»3  _  _i_ 

(8.62) 

(8.63) 

L(k»2  =  L«»  ,  -L. 

(8.64) 

These  are  the  only  non-vanishing  components  of  die  composite  thermal  coefficients  and  K^’s 
are  scalar  quantities  independent  of  the  coordinate  system.  Substituting  from  (8.61)  and  (8.62) 
in  (8.50H8.55)  we  get  the  following  contributions  to  the  thermal  parts  of  the  composite  stress 
tensors  and  moments: 

P 

^(thermal)  =-$0£  ~  Am f  +  Am,3) 

"♦lZPiCy  Am,2  +  %-Am’  +  p-4“i‘)  (8-65) 

^(thermal)  =  ^(thermal)  =  0  (8.66) 

^(thermal)  *  ^(thermal)  =  0  (8.67) 


^(thermal)  =  x^thermal)  =  0 


(8.68) 
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^(thermal)  =  t„(tbennal)  *  ■ -  +  £  PKAm,  +  y  Am?) 

F  F? 

-♦ilWfW+“  Am?)  (869) 

Seedhomal)  =  -  ♦«,  2  PK  y  W  +  yr  Am?  +  Am?) 

F?  F?  M 

-*i2  PKy  Am?  +  y-  Am?  +  —■  Am?)  (8.70) 

F  F2 

S2z(thcnnal)  =  -<|>0X  p,(y  Am?  +  y  Am?) 

£2  ?3 

2  My  Am,3 +^- Am,4)  (8.71) 

The  thermal  contributions  to  other  components  of  S  is  zero.  As  for  the  entropy  and  its  moments, 
heat  flux  and  its  moment,  we  substitute  from  (8.61M8.64)  in  (8.56H8.58)  and  get  the  following 
results 

PoTl(o)  =  (Tee + Y«  + 1„)  2  PKAm,  +  —  Am?) 

F  S2 

+  <««  +  «=>  2  My  Am,2 +-=-An1,3) 

F  F  F? 

+  4>o2  (peMAm,  +  y  Am?)  +  £  (pc)/(y  Am,2  +  y  Am?)  (8.72) 

F  F? 

PoTl(i)  =  (‘fee+Yiz+ Yn)2  PKy  Am?  +  —  Am?) 

F2  f3 

+  <*ee+ K=)S  My  Am,3  +  y  Am,4) 

BASE 


+ ♦.  2  (pc¥y  W + y  Am?) 

+*iE(P«>i(y  W+^W)  (873) 

M  *3 

Po^<2)  -  (%e + Y«  +  Y»)  X  fiK-j*  ^  +  ^  W) 

M 

+  (Kee  +  ^o)  X  PK“j*  Am  *  +  —  Am,5) 

£2  t3 

+  *.  2<PC«-|-  Am,3  +  •=-  An v4) 

M  M 

+  Cl  I<pcW-J-  Am,*  +  Am,5)  (8.74) 

<b=-i^-£kl(Am,+  ^Am,J)-i^-2kKYAm,I+^-Aml3)  (8.75) 

<b  =  -'^IkI(Ani,  +  ^.Am,2)-^-XkKy  W  +  -|-Am,3)  (8.76) 

<k=“8lX8|(4n,l+y  W)  (8.77) 

qie=-  7  %  X¥y  Am,2+  y.  Am,3)-  j  ^  I^y  Am,3+  ^  Am,4)  (8.78) 

<Ii,—  W  +  yAmft-^JWy  W  +  'yAm,4)  (8.79) 


It  is  worthwhile  to  recall  that  in  the  above  relations,  is  the  gradient  of  in  S3  =  r-direction. 


9.0  LINEAR  THEORY  OF  INITIALLY  SPHERICAL  LAMINATES 


For  an  initially  spherical  laminate,  we  choose  the  usual  spherical  coordinates  (r,+,8)  which 
are  related  to  the  coordinates  (G1,©2,©3)  according  to  the  following  relations 

eJ=0  .  e2=e  .  e3=r  (9.1) 

The  choice  of  coordinate  r  for  03  direction  is  natural,  since  here  similar  to  the  cylindrical  case, 
the  spherical  laminates  are  piled  up  in  the  r  direction.  The  metric  tensor  Gy  and  its  conjugate  Gij 
and  G  «  det(Gy)  are  given  by  the  following  relations 
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ii’i)*-#*-* 

(l  22)  =  (22l} 'COte^COt* 

(3  2  2>  =  (22  3}  =  -A"  =  ~ 

y  ^ 

(2 1 2)  =  ‘-sin  9  0’cos  01 »  -sin  $  cos  $ 
{2  3  2}  =  — ^sin2©1  =  -r  sin2^ 

From  (9.2)  and  (2.14)  we  find  the  following  expression  for  A 


2r  0  0  i2  0  n 

A  =  0  rsin2^  0  +  0  2r  sin2$  0  =  4r3sin2<fc 

0  0  1  0  0  1 

9.1  Relative  Kinematic  Measures 


(9.5) 


(9.6) 


The  relative  kinematic  measures  y,j  were  given  by  (2.118)1A3,  and  since  D0=Ga3  =  0  — 
by  (2.103)  and  (9.2)  —  and  8j  =  Uj)3  —  by  (4.22)  —  we  can  write 

TSj- y  (“iij  +  Ujii)  (9.7) 

By  (2.14)  and  (9.5)  the  non-vanishing  components  of  Ajj  are 

Ai'-A|=i  (9.8) 

As  a  result  we  have  the  following  expressions  for  given  in  (2.1 18)43 

Kop^  —  npio  +  Saip 

(9.9) 


BASE 
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k3o*^3»o 

The  components  of  director  displacement  6j  by  (4J2)  and  (9.5)  are 


x  1 

=  ”0.3-7  “a 

83  =  ”33 


(9.10) 


Now  we  proceed  to  calculate  the  covariant  components  of  ^  in  terms  of  the  covariant  com¬ 
ponents  of  the  displacement  vector  u  and  its  partial  derivatives.  From  (9.7)  we  can  write 


7  (”*  +  ”*>“  <1**1  “k 


(9.11) 


Using  (9.3)  we  get 


yn  =  uu  +  ru3 


Y12  =  7  (”u  +  ”2,1)  “  «*  ♦  Uj 
Y13 =  y(uu  +  u3,i)  ~  7  ui 

(9.12) 

Y22  =  ”2 3  +  sin^cos^Uj+r  sin2  $  u3 
Y23  =  7  (”23  +  u3,i)  ”  7  ”2 


Y»  =  ”33 

These  results  can  be  written  in  terms  of  the  physical  components  by  using  (8.4)  and  (9.2).  The 
appropriate  expressions  after  simplificiadon  are  as  follows 


1  d”* 

1’*=7(ir+u') 


BASE 
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v  1,1  **»,!*">  coc»  . 

*♦*  2  ^rsin*  36  r  8*  r  ^ 


v  =lA  +  i^L-^ 
2  (  dr  r  8*  r  * 


*»= 


1  Su^  .  cot 


r  sin  4> 


80  r  “  r 


J_  Sue  _  1  Su, 

Y*  2  dr  r  r  sin  $  80 ' 


(9.13) 


The  relations  for  and  given  in  (9.9),  after  using  (9.10)  are  reduced  to 

=  y  (up.a-Ua,p)-y  (pko)  uk  +  ua4p-{omp}unu3+{omp)  {mk3}Uk 
*3a  =  u3.3a“  (3ko)  ^3+{3ko)  U^)  uj 


In  expanded  form  these  relations  are  simplified  to 


Kn=tu3  +  u131+ru33 


*12  =  7  (“ll-U^  +  Ulf32  -  cot  *  uu 


*21  ”  7  («irun)  +  UW1  -  cot  4>  U23 


(9.14) 


*22  ss  sin2  $u3  +  U232  +  sin  cos  (j>  U13  +  r  sin*$  U33 


1  .  1 

*31  “  U3.31  -  7  U13  +  -jj  U! 


*32  =  U332  ”  “  u23  +  "jj  U2 


BASE 
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Using  (8.4)  and  (9.2),  we  can  write  the  physical  components  of  these  relative  kinematic  meas¬ 
ures.  The  results  after  simplification  are  as  follows 


**- H Wu'+-&)1 

_  _1_  due  +  1  d2**#  _  cot  $ 

i2  r  sin  $  d6dr  r  dr 


_  1  ^  .  1  ^6  cot  0 

^“r^sin*  30  r  3*dr  “  t2 


(9.15) 


,  .1  d,  l  ^ 

***  r  9r  sin  0  30  ^ 

The  equations  for  balance  of  mass  are  also  obtained  by  substituting  covariant  derivatives  of  dis¬ 
placement  vector  in  spherical  coordinates  in  the  expressions  (2.121)  and  (2.122).  The  simplified 
results  in  terms  of  physical  components  of  the  displacement  vector  and  its  derivatives  are  as  fol¬ 
lows: 


(9.16) 


P  =  PoO 


1  du^  i  due  du,  2  % 

- 


r  d<|>  r  sin  4>  36  dr 


BASE 
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9J  Linearized  field  equations 

In  order  to  derive  these  equations,  it  is  sufficient  to  substitute  from  (9.3)  and  (9.5)  in 
(2.123)-(2.125).  The  results  after  simplficiation  are  written  in  terms  of  the  physical  components 
of  various  tensors  by  using  (9.3)  and  are  recorded  below. 

Balance  of  linear  momentum: 


1  dlft,  .  v  cot*  .  .  1  . 

+  — :  + - «  .  -s—  + 


rsin$  dO  r  r 


i^sin^  dr 


spofii^  +  z1  8^) 


1  dl^e  ,  1  <^00  .  ^0r  .  cot<j>  ,  1  daz  L 

7"af  75n*  IT  T  “r^  ^  'ar+p°b9 


=  p0(iifl  +  z^5q) 


(9.17) 


1  d  ,  •  .  ,  ,  1  dx*  T^-Kee  i  dor 

iin<>  3<>  SU1  *  rsin$  36  r  ^sind)  dr  r 


spoCiiy  +  z1^) 


BASE 
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Balance  of  director  momentum: 


mT ^ (si” * ^ + 7^?  »  +  - + ^ + p»c* 


=  p0(z1ii^+  z2^ 


1  3s»e  1  dsee  .  cot  9  .  v  .  «er  .  qe 

r  9$  rsin$  96  r  s^®+s®^  r  ^^6  1,6 


p0(z1iie  +  z2^) 


(9.18) 


1  9  , .  A  v  1  4  %e  .  or 

- T—  -r-  (sin  <J>s*)  +  — —  — - +  -= - Xxr+Po 

rsin  <|>  9$  ^  rsin$  96  r  ^sin^ 

=  p0(z1'£I+z25r) 

Balance  of  moment  of  momentum: 


t^e  -  -  —  (s^-s^e) 


*1$  x4r®  r  V 


(9.19) 


^i©  ^©r =  r  Sftr 

Here  again  as  in  cylindrical  coordinates,  the  composite  stress  tensor  is  not  symmetric  and  the 
asymmetry  is  represented  in  equations  (9.19). 


BASE 


9J  Constitutive  relations 


We  follow  exactly  the  same  procedure  as  in  Section  8.3.  For  quantities  in  spherical  coordi¬ 
nates,  by  (9.6)  and  (9.3),  we  get 


and  by  (9.8) 


A  _  4r3sin2$  2 
2G  2i*sin2$  r 


(9.20) 


1 

r 


Using  these  results  in  (3.18)  and  (5.21)  we  obtain  the  following  contravariant  forms  of  the  con¬ 
stitutive  relations  in  spherical  coordinates.  These  constitutive  equations  are  obviously  written 
for  a  purely  mechanical  theory. 


T°  =  +  i- +  P<2^“  +  i- 


sop  _  p(l)apk/  +  2  jttlapk/  +  J_  +  p(2)c^r  +  2  j(3)oWr+  -L  i(4)a|J^K 

r  i2  r  t2  ** 


S«3  =  pOXOU  +  1  +  p(2)a^  +  ±  J?yaw ^ 


(9.22) 

(9.23) 

(9.24) 


The  covariant  components  of  the  relative  kinematic  measures  Jy  and  woe  calculated  and  are 
recorded  in  (9.12)  and  (9.14).  The  composite  constitutive  coefficients  I*k)pqr*  depend  on  the  con¬ 
stitution  of  the  laminates.  For  the  special  case  where  the  micro-structure  is  composed  of  n  iso¬ 
tropic  layers  with  different  elastic  constants  we  have 


Cr  =  \/)GPqG”  +  Mt/)(GprGq‘  +  Gp*Gqr) 


9-9 


(/= 

where  and  are  the  Lame  constants  of  each  layer  in  the  micro-structure.  By  (9.2)j  the 
non-vanishing  constitutive  coefficients  are 

C1111  =  (X+2n)/r* 

Cn22  =  X/r^sin2^ 

C1133  =  Vr4 

C2222  =  (X+^/^sin4^ 

C2233  =  X^sin2*  (9.25) 

C3333  =  X  +  2h 
C2323-^2* 
c1313=ma2 

C1212  =  p/^sin2^ 

where  the  subscript  (/)  is  dropped  for  brevity.  Of  course,  the  symmetries  of  the  constitutive 
coefficients  as  expressed  in  (8.28)  must  be  recalled  when  the  expressions  involving  the  con¬ 
stitutive  coefficients  are  to  be  expanded.  Substituting  from  (9.25)  in  (5.9)  and  using  the  results 
in  (9.21)-(9.24)  we  find  the  following  constitutive  relations  for  the  physical  components  of  the 
composite  stresses  and  composite  couples.  The  summations  in  these  relations  extend  over  the 
micro-structure  from  /  =  1  to  /  =  n. 


21} 

=  Z(X/ +  2p,)(Am,  +  Am2  +  -—£■  Am3) 


+  (Yee+Yir)  X  ^  Am,2  +  Am,3) 


BASE 
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+  k*Z  (X/f 2^X  y  W  +  y  W  +  W) 

+  K»ZMy  (9-26) 

v  243 

V  =  ^  X  M :W  +  —  Am^  +  y-  Am,3) 

1L  t2  t* 

+  (k^+k*)  £  My  W  +  y  Am,3  +  ~jj  W)  (9.27) 

V  =  ^  Z  WAm,  +  ^  Am?)  +  k*  £  Am?  +  -^  Am?)  (9.28) 


xe0  =  x«e 


(9.29) 


tee  =  (Y^+Yit)  X  XKAmi  +  -y  Am,2  +  Am,3) 

+  YeeL  (Xff2}i,)(Am,  +  -y  Am,2  +  y-  Am,3) 

+  **£^y  W+  y  Am,3+y-  Am,4) 

+  *»Z<V^AyAm?  +  ^Am?  +  |jAm?)  (9.30) 


%  =  2y*  Z  MrfAm,  +  y  Am?) 


+  «^Z)9<yAm?  +  ^Am,3) 


(9.31) 


BASE 


1 

1 
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■ 

1 

v  2t} 

V *  2Y*  I H*Am, +  yr  W +  -p-  A***3) 

1 

P  p2  ?3 

+  K*ZM*(y  Am,2+y  Am,3  +  ■—■  Am/*) 

(9.32) 

1 

1 

34n  2L? 

*i6  =  2Y.«Z  H/<Am,  +  —  Am3  +  —  Am,3) 

1 

P  p2  p3 

+  N«2M<(  2  Am,2  +  r  Am,3  +  ^  Am/*) 

(9.33) 

1 

1 

E_ 

V  =  (Y*  +  Yee)Z  tyAm,  +  —  Am,2) 

1 

1 

+  Yn  Z  (V*-2u,)(Am,  +  y  Am,2) 

P  2E3 

+  (K^+Kee)  Z  My  Am,2  +  -y  Am,3) 

(9.34) 

1 

1 

P  P^  P^ 

S*  =  Y*Z  (V2^Ky  Am,2  +  y  Am,3  +  -^  Am/*) 

1 

1 

1 

1 

+  (Tte+7n)  Z  ^Ky  Am,2  +  y  Am,3  +  -y  Am/*) 

p2  3p3  M 

+  z  (M-2fc)(y  Am,3  +  —  Am/*  +  y  Am,5) 

p2  3p3  M 

+  %z  Mj  Am,3  +  —  Am,4  +  y  Am,5) 

(9.35) 

1 

1 

S#e  =  2y^eZHKy  Am,2+y  Am,3  +  ~  Am,4) 

BASE 
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£3  t4 

+  (K*+K<*)  X  Mi(”  Am,3  +  —j—  Am,4  +  Am/5) 


f  21} 

S«r  -  2)^  X  M<(“2"  Am/2  +  Ami3) 
£2  E3 

+  **  X  MK-f- Am,3  +  Am/*) 


(9.37) 


S^arS^e 


(9.38) 


See  =  (Y*nr)  X  XK  i  Am,2  +  -2L  Am,3  +  Am,4) 


+  Yee  X 


£2  3^3  M 

+  £  V“J“  Am,3  +  -jj—  Am,4  +  Am,5) 

£2  363  e4 

+  *ee  X  (V2H|)(-j  Am,3  +  ^  Am,4  +  Am/5) 


(9.39) 


So. = 2*h  2  Mi(  Y  W + ^  W) 


£2  £3 

+%2M(<Y4,,,<3+'5-Am/4) 


(9.40) 


Using  (5.33),  (5.36)  and  (5.37),  die  composite  mass  density  p0  and  the  composite  mass 
moments  PqZ1  and  PqZ2  which  appear  in  the  equations  of  motion  are  also  calculated  for  an  ini¬ 
tially  spherical  laminate 


Po  =  X  PotijAm,  +  yX  Po/oAm,2 


(9.41) 


BASE 
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1  2£J? 

Po*1  =  y  4nl PotoW  fWW  (9.42) 

1  L3 

Po*2  -  J  1  Po(/yW  +  ^"Z  PoCoAm/1  (9.43) 

9.4  Energy  equations  and  constitutive  relations  of  linear  tbermodastidty  for  an  initially 
spherical  composite 

In  the  absence  of  heat  supply  or  heat  absorption,  the  energy  equations  (6.34)1i2  reduce  to 
the  following  forms  what  written  in  spherical  coordinates  as  defined  by  (9.1) 

ts*  [|r (sin  ♦  «•> + + 7  i (A) + p»<*>Wi>>  - 0  <9-44) 

1  d  dots 

Tita*  (s“  ♦  W  +  ~W]  +  h  -  % + p0(MrM>i%) = 0  (9.45) 

where  q$,  q@  and  q,  are  the  physical  components  of  the  heat  flux  vector  for  the  composite  lam¬ 
inate;  q^  and  q10  are  the  physical  components  of  the  composite  heat  flux  moment,  and  other 
quantities  have  the  same  meaning  as  section  6.  The  constitutive  relations  for  various  composite 
quantities  in  thermoelasticity  can  be  written  similar  to  what  was  done  in  section  8.4  for  initially 
cylindrical  laminates.  The  mechanical  parts  of  such  constitutive  equations  woe  deived  in  sec¬ 
tion  9.3.  Therefore  in  what  follows  we  record  only  the  thermal  parts  of  these  equations.  The 
complete  constitutive  relations  in  linear  thermoelasticity  are  obtained  by  adding  these  distinct 
parts.  Using  (7.10),  (9.8)  and  (9.20),  the  thermal  part  of  the  contravariant  stress  tensor  t'-’  is 

^(thermal)  =  -^{J^  +  j  J°*j  -*•  +  j-  J<2>*)} 

-(MJ^  +  J  J*2)ij  +  +  j-  F**)}  (9.46) 

or 

BAS 


t“  =  -^(J^1  +  f  J<w  +  ±  J™1)  -  +  +  i  J*3*1)  (9.47) 


t“  =  +Aj<^+4  J®*)  -  ^(J0)I2  +  2  J(2)i2  +  2  jOBj  (9>4g) 

*  r  r  r 


t° = -^0®“  +  -  P*3)  -  *i(*m  +  -  J®°) 


(9.49) 


r  r 

where  the  word  (thermal)  is  dropped  for  brevity.  Similarly  the  thermal  parts  of  the  composite 
couple  stress  are  written  by  using  (7.12),  (9.8)  and  (9.20) 


S«p  =  -4>0(J<1)a*  +  -  +  -j  J0^)  -  ^(J*2**  +  -  +  4  P**)  (9.50) 

r  r  *  r 


S'*3  =  -^(J*1*3  +  y  J*2**3)  -  4>i(J<2)a3  +  j  P*3) 


(9.51) 


The  constitutive  relations  for  composite  entropy  and  its  moments  T]m,  the  heat  flux  vector  q  and 
its  moment  qt  are  also  derived  by  substituting  from  (9.8)  and  (9.20)  in  (7.17),  (7.22)  and  (7.24) 


PoTl(m)  =  YijO^  +  j  J(nHl)ii)  +  tqp^1^  +  j  J<m+2)ip) 


+  <t»0(K(m)  +  y  K(m+1))  +  ^(K(m+1)  +  -  VP*2*)  (m  =  0,1,2)  (9.52) 

qi  =  -(L<0)#  + 1  -  (L<1»  +  j-  LWfyij  -  (L<°>°  +  j  L*1*3)^  (9.53) 

q«  «  -CL**#  + 1  L^XD^p  -«<*+!  L^14J  -  (L(I)a3  + 1  L®"*  (9.54) 

where  the  thermal  constitutive  coefficients  J®*,  K<k)  and  are  given  in  (7.6)-(7.8)  and  also  in 
(7.30)-(7.32). 


If  the  micro-structure  is  composed  of  isotropic  spherical  shells,  the  coefficients  of  thermal 
stress  and  thermal  conductivity  of  each  layer  can  be  represented  in  terms  of  only  one  constant. 
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Writing  (8.59)  and  (8.60)  in  spherical  coordinates  and  then  substituting  in  (7.30)  and  (7.32)  we 
get  die  following  results  for  the  non-vanishing  components  of  and  L***. 


■JEM"'*' 

(9.55) 

jW“-  >1. 

(9.56) 

(9.57) 

(9.58) 

L<B“=  W 

(9.59) 

L<lB3  *  kir  ^ 

(9.60) 

pm’s  and  k^’s  are  the  coefficients  of  thermal  stress  and  thermal  conductivity  of  each  layer  and 

the  summations  are  all  extended  over  the  micro-structure  from  l  =  1  to  /  =  n.  Substituting  from 

(9.55)-(9.57)  in  (9.47M9.49)  we  get  the  following  expressions  for  the  thermal  parts  of  die  physi¬ 
cal  components  of  the  composite  stress  tensor. 

■V - £  M  Am,  +  ^  Am,2  +  -|§-  Am,3) 

— <t>i  I  My  Am,2  +  -^  Am,3  +  Am,4) 

(9.61) 

o 

II 

$ 

II 

$ 

(9.62) 

© 

II 

II 

e* 

(9.63) 
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3?  2 El 

tee  =  -^oZ  P/(Am,+  —  Am,2  +  -p 


X  PX  Y  Am/2  + Am,3  +  p- Am/4)  * 

(9.64) 

vsv° 

(9.65) 

^=^oXPXAm/4  Am,2)-*,  X My  Am 2 4-  Am,3) 

(9.66) 

The  thermal  part  of  the  composite  stress  moment  is  also  similarly  calculated  by  substituting  from 

(9.55M9.56)  in  (9.50)  and  (9.51). 

s«  =  -$0P|(y  Am,2  +  Am,3  +  Am,4) 

I  ftOy-  M3 *  ^  W> 

(9.67) 

S^e  =  Se^=0 

(9.68) 

Jt 

ii 

$ 

CO 

(9.69) 

S#r=S0r-O 

(9.70) 

In  order  to  find  the  appropriate  forms  of  the  constitutive  equations  for  the  entropy  and  its 
moments,  heat  flux  and  its  moments  we  substitute  from  (9.55H9.57)  and  (9.58H9.60)  and  also 
(7.31)  in  (9.52)-(9.54)  and  obtain  the  following  results 

E_ 

Po^tO)  =  (Y(#+Yeer+'yrr)  X  PKAm,  +  —  Am/2) 

+  (**+*ee)  X  PKy  Am/2  +  Am,3) 
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p  p  2t* 

+  ♦o  2  (pc)j(AnJ|  +  —  Am/2)  +  $!  X  (pc)/(—  Am/2  +  —  Am/3)  (9.71) 

p«n(i) = (y#^B)  S  PK  y  W  +  W) 

+  (K^+Kee)  X  P/(  Y  Am/3  +  Y 

+  ♦o  2  (pc)K  Y  W  +  ~~  Am/3)  +  9]  X  (pc)/(  Y  Am/3  +  Am/4)  (9.72) 

P2  f3 

Po^a) =  (Y^+Yee'+Tir)  2  P/<  Y  +  Y 

F3  2tf 

+ <%«*•> I PK  Y  ^ + Y*  W) 

£*  *3  z3  2M 

+  0oZ  (pc)K  Y  W  +  Y-An^  +  9,X  (pc¥  Y  Am/4  +  Am/5)  (9.73) 


q4=_7'  ^^k/(Arn/+'7' 


^1  — ,  .  y.  2  ^  *  3. 

—  2k/(Am/2T  +  — Am/3) 


q9=-7^^£k^Am'+7l'ta'2) 

- —  X  kK—  Am,2  +  —  Am/3) 

rsin*  d0  ^  A  2  '  3r  " 


9r  =  “4i  2  k/(Am/  +  —  Am,2) 


(9.74) 


(9.75) 


(9.76) 
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1  frfro  ySa  A  2  ,  *  3\ 

1  d6t  L?  £* 

~7  Y  w+  £  W>  C9.77) 

1  ^9  ,^R  ! _ 2  .  A_,* 

qi9=~75Ki  lT2¥TAnl'  +-3TAn,') 

-^♦^•2,q<fto'3+¥A,n'4)  (978) 

It  should  be  recalled  that  in  the  above  relations,  ^  is  the  gradient  of  $D  in  the  03,  Le.,  r-direction. 


10.0  ANALYSIS  OF  COMPOSITE  LAMINATES  FOR  IN-PLANE  LOADING 

10.1  Introduction 

In  the  previous  sections  a  complete  thennomechanical  theory  of  composite  laminates  was 
developed.  In  this  section  the  results  of  stress  analysis  of  composite  laminates  with  traction-free 
edges  are  presented.  Composite  laminates  with  traction-free  edges  are  known  to  develop  inter¬ 
laminar  stress  concentrations  near  the  edge  region.  The  problem  of  a  finite  width,  symmetrically 
laminated  composite  plate  under  uniform  one-dimensional  stretch  has  been  studied  by  many 
authors,  see  Section  1.  Pagano  and  Pipes  (1970,  1973)  showed  that  free  edge  effects  on  the 
interlaminar  stresses  are  important  issues  in  determining  the  failure  and  the  strength  of  such  lam¬ 
inates.  Their  analytical  work  was  based  on  linear  elastic,  generalized  plane-strain  formulation 
and  numerical  solutions  were  obtained  using  a  finite  different  procedure.  Their  study  revealed 
that  certain  interlaminar  stresses  rise  in  magnitude  near  die  free -edge  region.  It  was  suggested 
that  a  possible  stress  singularity  exists  at  the  free  edge.  A.  S.  D.  Wang  (1977)  followed  a  finite 
element  scheme  to  investigate  the  same  problem  with  emphasis  placed  on  assessing  in  detail  the 
stress  field  closest  to  the  ply  interfaces  and  laminate’s  free  edge,  where  stress  singularity  is 
suspected.  In  this  section  the  analysis  of  the  same  problem  based  on  the  theory  developed  in 
previous  sections  is  presented.  A  finite  difference  scheme  was  adopted  for  the  solution  of 
governing  partial  differential  equations.  The  objective  of  this  numerical  modeling  study  was  to 
examine  the  three-dimensional  state  of  stress  at  the  free  edges  of  composite  laminates  and  to 
show  that  the  proposed  theory  is  in  agreement  with  recorded  experimental  data. 

10.2  Free  Edge  Boundary  Value  Problem 

Consider  a  prismatic  symmetric  laminate  rhown  in  Figure  10.1. 
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The  corresponding  boundary-value  problem  for  a  uniform  strain  field  in  x-direction  and 
traction-free  edges  at  y  =  ±b,  and  top  and  bottom  surfaces  (z  =  ±h)  were  derived  based  on  the 
linearized  field  equations  (2.126)  and  (2.127)  and  constitutive  relations  (3.24)  and  (3.31).  It  was 
assumed  that  representative  elements  are  made  of  orthotropic  plies.  The  stress-strain  relations  of 
each  ply  in  a  coordinate  system  with  major  axis  along  fiber  direction  is: 

x*  =  Cy  (10.1) 

where 

*  *  *  * 

T1  “''ll  ♦  x4  -t23 

*  *  *  * 

X2  -  T22  *  x5  “  t13 

*  *  *  * 

t3  -  X33  ,  X6  -  X12 

Yl  =  ul.l  .  Y4  =  U23  +  U3.2 
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T2  =  u2J  .  Ys*»U  +  u3.1 
Yj  =  u33  •  Yfc  =  uU  +  u2,l 

and  C  is  a  (6  x  6)  matrix  with  elements  related  to  the  nine  material  constants  of  each  ply,  i.e., 
extensional  elastic  moduli  in  three  directions,  shear  moduli  and  Poisson  ratios  in  corresponding 
directions  as  shown  in  [Whitney,  1989].  The  constitutive  relations  for  composite  stress  and 
composite  stress  couple  were  derived  by  rotation  of  coordinate  system  in  (10.1)  to  coincide  with 
the  direction  of  axial  loading  and  the  integration  of  these  relations  across  the  thickness  of  the 
representative  element  as  shown  in  equation  (3.24)  and  (3.31).  For  small  deformations  of  flat 
composites,  the  constraint  relation  (4.20)  was  employed,  i.e., 

6(0a,03,t)  =  — @°,e3)  (10.2) 

The  final  farm  of  constitutive  relations  for  the  composite  stress  and  the  composite  stress  couple 
assumed  the  following  presentation: 


x  =  Cy+  Dk 

(10.3) 

S  =  Dy+Fie 

where 

C  =  mC(1)  +  (l-m)D(2) 

D  =  y  [m2C(1)  +  (l-m2)C<2)]  (10.4) 

t2 

F=^-[m3C(1)  +  (l-m3)C(2>] 

K1  -  ul,13  »  *4  =  U3>23 
k2  =  u2,23  *  K5=u3t13 
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*3  =  0  ,  Kg  »  Uj^3  +  02,13 

In  (10.4),  C(d  and  C(2)  are  corresponding  C  matrices  fen-  each  constituent  (or  for  each  fiber  direc¬ 
tion)  present  in  die  representative  element  and  £2  is  the  thickness  of  representative  element  and 
m  is  defined  in  (3.15).  The  linearized  field  equations  (2.126)  and  (2.127)  for  a  static  loading  and 
in  the  absence  of  body  force  are: 


S^+oi-t^O 

Elimination  of  the  interlaminar  stress  vector  &  from  these  equations  resulted  in: 


(10.5) 


(10.6) 

For  an  axially  loaded  strip,  the  stress  and  stress  couple  components  were  taken  to  be  independent 
of  the  axial  direction.  Consequently  the  general  form  of  the  displacement  field  was  assumed  as: 


u1  =  ax1  +  U(x2pc3) 

u2  =  V(x2tx3)  (10.7) 

u3  =  W(x1pt3) 

Identifying  direction  1  with  x,  direction  2  with  y  and  direction  3  with  z,  the  field  equation  (10.6) 
reduces  to  the  following  set  of  partial  differential  equations: 


C«,U.y,  +  C55U^  +  CjjVyy  +  C„V  =  +  (Ctf+C^VVy, 


+  (1^45—1^36) W.yzz  “  ^toU.yyiz  “  ^26^,yyu ~  0 


C^U.yy  +  04511^  +  C22v,yy  +  CuV'iz  +  (C23+C44)Wiy2 


(10.8) 
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(Qj^CsfiKJyx  +  (C^f+C^Vyj  +  C44W>yy  +  C33W  ^ 

+  P36“D45)U.yzx  +  (D23-D44)Vtya 
"  F44W,yyzz  =  0 

where  Cy,  Dy  and  Fy  are  the  components  of  C,  D  and  F  matrices.  The  domain  and  the  boundary 
conditions  for  these  equations  are  shown  in  Figure  10.2. 


Figure  10.2 


Fori: 


1^  =  0 
V^O 
w=o 


Forll: 


Uv  =  0 

•wy.o 

V=0 

b 


x2i  =  0 


Forlll: 


*  X22  =  0 

X23  =  0 

b 


For  IV:  < 


Oj  =  0 

c2  =  0 
c3  =  0 


The  following  material  properties  were  used  for  the  analysis: 


Ex  =  48  x  104psi 
Ey  =  E2  =  5xl04psi 
GXy  =  Gyz  =  Gw  =  2  x  lC^psi 
VXy  —  Vjy  =  VXZ  =  0.21 
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A  rectangular  mesh  of  41  x  1 1  nodes  was  used  for  the  finite  difference  discretization  and  a  four 
layer  symmetry  laminate  under  a  uniform  axial  strain  yx  «  0.01  was  considered.  Each  layer  itself 
could  be  a  collection  of  thin  plies  repeated  in  a  consistent  pattern.  The  presented  results  are  for 
the  case  that  all  plies  in  the  top  and  bottom  layers  are  in  the  +6  direction  and  all  plies  in  the  two 
middle  layers  are  in  the  -8  direction  providing  a  [±6],  laminate.  Complete  stress  and  displace¬ 
ment  results  were  obtained  for  various  values  of  6.  These  results  are  presented  in  Figures  10.3 
through  10.13  and  summarized  in  the  following. 

10.3  Results  of  Finite  Difference  Simulation 

The  purpose  of  this  simulation  was  to  examine  the  response  of  composite  laminates  under 
uniaxial  extension  and  to  show  that  the  proposed  theory  reflects  the  complex  three-dimensional 
response  of  free  edge  problem  in  composite  laminates  as  recorded  in  the  literature.  Following 
this  verification,  a  systematic  discretization  technique  in  the  context  of  finite  element  method 
was  developed  and  extensive  analyses  simulating  various  flat  and  curved  composite  laminates 
under  in-plane  and  out-of-plane  loading  were  performed. 

In  Figures  10.3  through  10.9  various  components  of  stress  tensor  and  interlaminar  stress 
vector  for  [±30]„  [±45]4  and  [±60]ig  laminates  are  plotted  along  the  symmetry  line  of  the  lam¬ 
inate.  Figure  10.3  is  the  axial  stress  which  shows  a  decrease  at  the  free  edge.  For  [±30]4  lay-up 
this  decrease  is  about  50%  of  stress  at  the  centerline  y  *  0.  Figure  10.4  is  the  composite  shear 
stress  Txy  which  assumes  its  maximum  for  0  =  30°  and  approaches  zero  at  the  free  edge. 

Figure  10.5  is  the  normal  in-plane  composite  stress  in  the  y  direction,  perpendicular  to  the 
loading  axis.  The  value  of  this  stress  component  is  negligible  as  compared  to  the  axial  stress, 
about  0.3%  of  the  axial  stress. 

Figure  10.6  shows  the  interlaminar  normal  stress  in  z  direction.  The  value  of  this  stress 
component  assumes  its  maximum  on  the  free-edge  boundary. 
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Figures  10.7  and  10.8  are  the  interlaminar  shear  stress  components.  The  magnitude  of  the 
interlaminar  shear  stress  along  the  x-direction  is  in  the  order  as  normal  interlaminar  stress 
and  it  increases  with  a  high  gradient  as  it  approaches  the  free  edge.  The  magnimdc  of  the  yz 
component  of  shear  stress  at  the  centerline  is  negligible  compared  to  other  components. 

Figure  10.10  shows  the  variation  of  interlaminal  shear  stress  t u  for  various  values  of  fiber 
direction  0. 

Figures  10.11  through  10.13  present  the  various  stress  components  across  the  thickness  of 
the  laminate.  The  shear  stress  T^,  shows  very  small  variation  across  die  thickness.  Hie  normal 
interlaminar  stress  assumes  its  maximum  along  the  symmetry  line  and  it  approaches  zero  on  the 
top  and  bottom  surfaces  of  the  laminate.  The  normal  stress  Xyy  assumes  its  pick  on  the  top  and 
bottom  layers.  Its  value  at  the  symmetry  surface  is  not  zero  but  it  is  considerably  smaller  than 
those  values  at  top  and  bottom  surfaces. 

A  more  detailed  study  of  this  problem  is  presented  in  Chapter  14  based  on  a  finite  element 
scheme.  In  particular  it  is  discussed  that  even  for  a  symmetric  laminate,  the  problem  of  in-plane 
loading  of  a  finite  width  strip  is  a  three-dimensional  problem  and  should  be  modeled  accord¬ 
ingly. 
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Figure  10.5 

Extension  Analysis  —  Normal  Stress  (Tyy) 
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Figure  10.6 

Extension  Analysis  —  Interlaminar  Normal  Stress  (O3) 
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Figure  10.7 

Extension  Analysis  —  Interlaminar  Shear  (Oj) 


Interlaminar  Shear  Stress  yz  Component 
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Figure  10.8 

Extension  Analysis  —  Interlaminar  Shear  (02) 
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AxiaJ  Displacement  at  Composite  Surface 
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Figure  10.9 

Extension  Analysis  —  U(y,z) 


interlaminar  Shear  Stress  Vs  Theta 


Figure  10.10 
Effect  of  Fiber  Orientation 
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Composite  Shear  Stress  (xy  Component) 


Symmetry  Surface . -  Mid  Surface  .  Top  Surface 


Figure  10.11 

Variation  through-the-thickness  (xxy) 


Composite  Normal  Stress  (yy  Component) 


- Symmetry  Surface  .  Mid  Surface  .  Top  Surface 


Figure  10.12 

Variation  through-the-thickness  (tyy) 
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Composite  Normal  Stress  (zz  Component) 


- Symmetry  Surface - Mid  Surface  . Top  Surface 


Figure  10.13 

Variation  through-ther  thickness 
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11.0  WAVE  MOTIONS  IN  LAMINATED  FLAT  COMPOSITES 

Using  relations  (5.38H5.42),  the  linear  equations  of  motion  for  a  laminated  composite  with 
initially  flat  plies,  in  the  absence  of  body  forces,  can  be  written  as 


la&Utfa  +  +  Ojj  =  p0Uj  +  Po^iip  (H.l) 

+  la&Ufc#3  +  Oj  -  *  Po^Uj  +  Po*2^  (H-2) 

Eliminating  Oj  between  these  equations  we  get 

-  I<$Acja3  -  =  Poiij  -  Po*2^  (H-3) 

The  equations  (11.3)  are  the  differential  equations  for  the  displacement  vector  u  in  elasto- 
dynamical  problems.  These  equations  are  now  used  to  investigate  the  propagation  of  small 
amplitude  harmonic  waves  in  a  laminated  composte.  In  the  following  special  cases  that  we 
examine,  k  is  the  wave  number  and  c  the  phase  velocity  of  the  appropriate  wave 

(a)  Longitudinal  Waves  in  the  xj-direction 

For  waves  of  this  type  the  non-zero  displacement  component  is  Uj  and  we  have 

U!  =  Aj  exp[ik(xj-ct)]  ( 1 1 .4) 

where  Aj  is  the  constant  wave  amplitude  and  assumed  to  be  small.  Differentiating  (11.4)  with 
respect  to  xj  and  t  we  get 

fii—kVuj  ,  ultll  =-k2u1  (11.5) 

Substituting  (1 1.5)  in  (1 1.3)  we  obtain 


Po 


(11.6) 


For  a  composite  whose  micro-structure  is  composed  of  n  isotropic  layers,  by  (5.27)  and  (5.26)j 
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(5.26)i  we  have 

iffli  -  i  (\,)  +  -2M<,)Wnv  (H.7) 

(-1 

Substituting  this  result  together  with  (5.38)  in  (11.6),  the  wave  speed  c  for  a  longitudinal 
wave  would  be 


c2  =  •=!-; -  (11.8) 

Z 

1=1 


(b)  Horizontally  polarized  shear  waves  in  the  Xj-direction 

For  this  type  of  waves  the  non-zero  displacement  component  is  u2  and  we  have 


u2  =  ALexp[ik(xi-ct]  (1 1.9) 

Substituting  from  (11.9)  in  (11.3)  we  get  the  following  expression  for  the  wave  velocity 


Po 

which  for  the  special  case  of  isotropic  laminates  reduces  to 


(11.10) 


ZiVH 

C2^ -  (11.11) 

Z 

1=1 


(c)  Vertically  polarized  shear  waves  in  the  x1-direction 


The  non-zero  displacement  component  for  this  wave  is  U3  and  we  have 
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u3  =  A3  exptikfxpct)] 


(11.12) 


The  wave  velocity  for  this  case,  similar  to  the  above  cases,  is  found  and  we  obtain 


c*  = 


1^3 


Again  for  the  special  case  of  isotropic  laminates  we  get 


(11.13) 


C2*-^ - 

i  p 

i=i 


(11.14) 


(d)  Longitudinal  waves  in  the  x3-direction 

In  this  case  the  non-zero  displacement  component  is  u3  and  we  have 

u3  =  B3  exp[ik(x3-ct)]  (11.15) 

These  waves,  unlike  the  above  three  cases,  are  dispersive  and  die  wave  speed  depends  on  fre¬ 
quency.  The  non-zero  space  and  time  derivatives  of  (11.13)  which  are  relevant  to  (11.3)  are 


ii3  =  -k7c2u3  ,  u3i33  =  -k7u3  ,  ii3>33  =  kVu3 
Substituting  from  (11.16)  in  (11.3)  we  get 


(11.16) 


=  PoC2  +  PoZ2^ 

If  we  introduce  die  wave  frequency  0)  =  ck  in  (11.17)  we  get 


c*  = 


iffis-PoZ2®2 


Po 


(11.17) 


(11.18) 


For  the  case  of  a  composite  whose  micro-structure  is  composed  of  n  isotropic  layers,  this  rela¬ 
tion  reduces  to 
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c2 


»  (toL)2  *  M  X 

I(\r)+2mr))A“r - 5~  2  Po^W 

i*l  J  i*l _ 


£  Po*0^ 

i*i 


(11.19) 


(e)  Transverse  shear  waves  in  the  x3-dircction 

In  this  case  we  consider  a  transverse  shear  wave  propagating  normal  to  the  laminates  with 
its  amplitude  in  the  Xj-dirccnon.  Consequently  the  only  non-zero  displacement  component  will 
be  Uj  and  we  have 


Uj  =  B  1exp[ik(x3-ct)]  (11 .20) 

Here  again,  the  phase  velocity  c  is  obtained  similar  to  the  case  (d).  The  result  is 


c2_  Po^p2 

r  Po 


(11.21) 


which  shows  the  dependence  of  c  on  frequency  co.  For  the  case  of  isotropic  laminates  (1121) 
can  be  written  as 


n 


2P<rM- 

1=1 


(O2 

3 


£  powa-»,3 
1=1 


£  P^Am, 

1=1 


(11.22) 


It  should  be  noted  that  the  phase  velocity  of  a  shear  wave  propagating  in  the  x3-direction  with 
amplitude  in  the  Xj-direction  can  be  obtained  by  substituting  I^j  in  place  of  l|3]3  in  the  relation 
(11.21).  The  general  solution  of  a  transverse  shear  wave  propagating  normal  to  the  laminates  is 
the  sum  of  these  solutions. 
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12.0  WAVE  MOTIONS  IN  CYLINDRICAL  AND  SPHERICAL  LAMINATES 

The  results  of  sections  8  and  9  are  used  to  derive  equations  of  motion  in  terms  of  displace¬ 
ment  vector  u  for  the  cylindrical  and  spherical  laminates,  hi  each  of  the  two  cases  die  micro¬ 
structure  is  supposed  to  consist  of  n  isotropic  layers  with  different  elastic  constants.  Using  the 
same  notations  as  previous  sections  we  define  the  constitutive  coefficients  X(r)  and  p(r)  according 
to  the  following  relations: 


*(r)  =  5rMIVW  02.1) 

(r=  U,...) 

i*<r)=5r12><to>j  (i2.2) 

where  summations  extend  over  the  micro-structure  from  /  =  1  to  /  =  n.  For  future  reference  we 
further  introduce  the  quantities  p  ^  related  to  the  densities  of  different  layers  according  to  the 
following  relations 


P.w  -  y1  Z  <r  =  1 A-)  (12.3) 

where  summation  again  extends  over  the  micro-structure  from  /  =  1  to  /  =  n.  The  quantities 
X(r)i  p(r)  and  p*r)  defined  in  (12.1)-(12.3)  are  known  a  priori  for  each  composite  laminate. 

12.1  Governing  Equations  for  Cylindrical  Laminates 

Using  the  results  of  section  8,  the  equations  of  motion  are  derived  for  axial  symmetry.  In 
other  words  we  will  study  motions  which  are  independent  of  the  axial  coordinate  z  and  the  angu¬ 
lar  coordinate  0.  With  these  specifications  we  will  not  have  variations  in  01  =  0  and  03  =  z  direc¬ 
tions.  From  relations  (8.12)  and  (8.13)  we  calculate  the  physical  components  of  the  relative 
kinematic  measures: 


n  1  ,due  Uq 

Yee  =  —  •  Ye*  =  0  ,  Yer  =  j(-^ - — ) 
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(12.4) 

Kzz  =  0  ,  ^  =  -7  •  *n  =  0 

The  equation  of  motion  in  terms  of  the  physical  components  of  the  composite  stress  tensor  and 
the  composite  stress  couple  are  derived  using  relations  (8.15)  and  (8.16).  These  equations  are 
written  in  the  absence  of  body  force  and  body  couple. 


1 


du, 


Yrz  =  0  ’  2  *  y"=  & 


*06=  (^t)  .  *0**0  »  *18®° 


1  1  do*  ... 

7  %  +  7  77- =  Po(ii«  +  *  5e) 

7  ■^■-Po(“2  +  zlSI) 

"  ~  Tee  +  7  77  =  P0(«r  +  z%) 
78^+7  00-1,0= poCr'iie + z25e) 


(12.5) 


70z-''re=Po(zlUz+Z28z) 

-7500+70,-^  =  p^z1^  +  rS,) 

The  constitutive  relations  for  various  components  of  the  composite  stress  tensor  and  the  compo¬ 
site  stress  couple  are  derived  using  relations  (8.29)-(8.43)  along  with  definitions  (12.1)  and 
(12.2) 


tee  =  {k0>  +  2p(1)  +  7  (X(2)  +  2p(2))  +  -jj  (X05  +  2p(3))}yee 
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+  <*«+|x»+.^Lx»>fr 

+  { \ (X0)  +  v?*)  +  ~  <X0)  +  2*i0))  +  -4  (X<4)  +  2^(4>)}Kee  (12.6) 

2  3r  4r* 

t*  =  20l«>  +  -j-  H®  +  p  H®)*,  +  ( J  H®  +  £  (*®  +  p  H<4>)«*  (12-7) 

Tn  =  (X(,)  +  -i  X®**, + Hie)  +  <}  X®  ■ +  X0)»Cm  (12.8) 

V=2fti<'>  +  i(l0))r„  (!2.9) 

te,"201<»  +  i|i®)ftr+(|(i®+iti'3>)«rt  (12.10) 

To  =  20i(')+in®)rlr  (12.11) 

V={X<l|  +  2n<l)  + +  2)i<s)  )Y„ 

+  (X<»  +  i  X®)fo  +  (y  X®  +  i  X'3*)**  (12.12) 

s«e=  (y(X®  +  2n®)  +  ^:(X®  +  2Ii<3>)  +  ^  (Xm>  +  2^4>))%9 

+  (IX®  +  |-).®  +  ix'4>)y„ 

+  ( -r  (>.(3)  +  2n'31)  +  (X<4>  +  2n<41)  +  -jr-  (X®  +  V5*))^  (12.13) 

3  2r  5i 
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V = 2(|  K®  +  It'41)*,,  02.14) 

«*  “  (1 X<2’ +  ¥  X<3>>(1'® + «  +  <| X<3' +  ^  X<‘W  (1215) 

V'2(7t‘,!|+^%  (12,16) 


*iB  ~  ^  ^  *  *81  ”  ®  (12.17) 

It  should  be  mentioned  that  with  the  constitutive  relations  (12.6M12.17),  the  equations  (8.17)  for 
balance  of  the  angular  momentum  are  identically  satisfied. 

The  composite  mass  density  p0  and  the  composite  mass  moments  PoZ1  and  PqZ2  are  calcu¬ 
lated  using  relations  (8.44)-(8.46)  together  with  the  definitions  (12.3) 


«.-(*>+ £pP 

P„z'  =  y  Po®  +  Po0>  0218) 

wJ=jp»<31+ip»4> 

As  for  the  physical  components  of  the  director  vector  8,  using  (4.22)  in  conjunction  with  (8.1), 
(8.2),  (8.5)  and  (8.7)  we  find  the  following  results 


(12.19) 


In  order  to  substitute  the  constitutive  relations  (12.6)-(12.17)  in  the  equations  of  motion,  first  we 
eliminate  the  interlaminar  stress  components  oe,  o2  and  or  from  the  set  of  six  equations  (12.5). 
The  result  would  be  the  following  set  of  three  equations  of  motion.  In  these  equations,  as  previ- 


BASE 


ous  similar  equations,  double  dot  denotes  second  partial  derivative  with  respect  to  time  t 

-  'p"’  -  i  p°0)  i  -  <i  + i  p«4>)  ^ 

i  ^  i-(i  ***♦  i  ^  JlH  <1220> 

^^+7C«w-W+7|r*»- 

In  deriving  these  equations  we  have  also  used  the  results  (12.18)  and  (12.19).  Now  the  relevant 
constitutive  relations  are  substituted  in  equations  (12.20).  After  straightforward  calculations  and 
some  simplifications  we  get  the  following  results  for  the  left-hand  sides  of  these  equations: 

+  20i<»  +  ±  ,1®)  ^  +  (I  n®  +  „®>^  . 

±^±(^,+±±^1  (y»- 1  x®-  x®*. 

+  (X<»  +  2^<»  +  i  (X®  +  |1®)  +  X®  +  Ar  X«>)  4-  Y» 

r  3r  4r  dr 

-  j-  (2n(1)  +  ~  (X(2)  +  2p(2))  +  -j(*<3>  +  2p(3))  +  (X<4>  +  2p(4)}y90 


(12.21) 
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+  QP* + 24<5>)J,ii 

•  (bS0  -  j  fc®  (j  P.® + i  Pi4>)  ^ 02*) 

It  is  interesting  to  note  that  these  equations  are  uncoupled  in  terms  of  the  displacement  com¬ 
ponents  and  represent  three  types  of  transient  wave  motions  with  axial  symmetry.  Equation 
(12.22)  is  the  governing  equation  for  rotary  shear  motions.  Equation  (12.23)  represents  axial 
shear  waves,  and  finally  equation  (12.24)  is  the  dynamical  equation  for  radial  waves.  Since  the 
present  theory  is  designed  in  such  a  way  that  the  classical  theory  for  a  homogeneous  continuum 
can  be  derived  through  a  limiting  procedure  in  which  the  thickness  of  the  micro-structure 
approaches  zero,  we  expect  to  reconstruct  the  equations  for  wave  motions  with  axial  symmetry 
in  a  homogeneous  isotropic  medium  by  letting  ^  ->  0  in  equations  (12.22M12.24).  Doing  so, 
the  only  non-vanishing  constants  in  the  coefficients  of  these  equations  are  those  with  superscript 
(1)  and  we  obtain  the  following  results. 

Rotary  shear  waves: 

(12.25) 

(12.26) 

(12.27) 

BASE 


l  Sue  u$  _  i  92ue 
dr2  r  dr  r2  c*  dt2 

Axial  shear  waves: 

i  dug  i  d2ur 

dr2  r  dr  cy  dt2 

Radial  waves: 

1  ^  u,_  l  d2^ 

dr2  +  r  dr  “  r2  "  c£  dt2 

where 
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+  (X<»  +  -  (X<»  +  n«»)  +  \  <X<”  +  2n«>)  +  -L  (X«>  +  2n«)))  ±  -to 
r  3  r  4r  or 

-  -  {4®  +  -|-(X(3)  +  2h(3>)  +  (X(4)  +  2p(4>)  +  -2j  (X<5)  +  2{i(5))}K^e 

r  3r  4r  5r* 

+  (|x®+|:(X0>  +  |iP>)  +  -i-(X«>  +  2|»'4>)  +  -L  (X<» + 2|1<5>) 

Finally  we  substitute  for  the  relative  kinematic  measures  from  (12.4)  and  obtain  the  following 
equations  for  the  components  of  the  displacement  vector  u: 

o.®  -  -i.  p«)  ^  +  j  (p<’> + i  n®>  ^  i  n<% 

= ‘fc®  -  £  p®  i  -  <?  p",+ i ««  0^2) 

(4o) + J.  m(2>,  + 1  ^(1) 

=  tpW- i  P„0)-(}  pf  +  i  P<4>)  (12.23) 


{K(1)  +  2p(1)  +  ~  OX®  +  2ji®>  +  ~  (2X®  +  ^  +  —  (3X(4)  +  4p(4)) 
2r  3r  4x3 


+  -^-(X(5)  +  n(5))} 
5r4 


d2Uj 

a 7 


+  i  «X<»  +  2p<»  -  4  X®  -  4(X<4>  +  p<4>)  -  -4  (X<«  +  2p<!>»  % 
r  r  r  5r  or 


-  \  [X(l)  +  2p(1)  +  -  (3X(2)  +  4|i(2))  +  ~  (X(3)  +  2p(3))  +  (X(4)  +  2h(4)) 

r  r  3r  2r 
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.  .I-A 

Po  Po 

(12.28) 

and  we  have  deleted  die  superscript  (1)  for  simplicity.  These  are  the  familiar  equations  for  wave 

motions  with  axial  symmetry  in  cylindrical  coordinates.  Equations  (12.25)  and  (12.27)  are  of 

the  same  type  and  the  space  part  of  their  solutions  can  be  represented  in  terms  of  Hankel  func¬ 
tions  of  the  first  order.  The  space  part  of  the  general  solution  of  equation  (12.26)  can  be  written 

in  terms  of  Hankel  function  of  zero  order.  For  each  case  considering  a  solution  of  the  general 

form  F(r)e10*  we  get  the  following  results: 

Ufl(r.t)  =  (C,H|‘>(^-  r)  +  C,H,0)  (-£-  r))e'“ 

(12.29) 

u,(r,t)  =  {C3H»>(-!?-  r)  +  C.H®  (-£-  r))c™ 

Lj  Up 

(12.30) 

u,(r,t)  =  (CjHf^  r)  +  C^HP  (7?-  DJe1” 
t>r  Or 

(12.31) 

where  H^  and  H^2)  (n  =  0,1)  are  Hankel  functions  of  first  and  second  kind. 


12.2  Investigating  Wave  Motions  in  Cylindrical  Laminates 

Returning  back  to  the  equations  (12.22)-(12.24)  we  can  again  use  the  technique  of  separa¬ 
tion  of  variables  to  get  the  appropriate  differential  equations  for  the  spatial  part  of  the  wave 
motions. 

For  the  rotary  shear  motions  if  we  let 

Ue(r,t)  =  F(r)ei<ot  (12.32) 

in  the  equation  (12.22)  we  will  obtain  the  following  second  order  ordinary  differential  equation 
for  F(r) 

BASE 
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(J1(1)  -  \  <02po(3) -  J-  afp™  -  p(3)r,(r)  (12.33) 

3  4r  3r 

j  (p0)  -  j  afp®  +  j-  ptyFir)  +  (eo*p <*> -  -j  p°>)F(r)  =  0 
This  equation  can  be  written  in  the  following  form: 


where 


F"(r)  + 


qr+U(2) 

co^+prfy 


F(r)  + 


aPpVV-tiW 

otr^+Pr-py 


F(r)  =  0 


(12.34) 


a  =  |i(,|-|o1p®  ,  p  =  -^jp-  ,  y=-j-vn>  (12.35) 

The  origin  r  =  0  is  an  ordinary  point  for  this  differential  equation;  therefore,  die  solution  for  F(r) 
is  the  neighborhood  of  r  =  0  can  be  written  in  the  form  of  an  infinite  power  series  of  r.  The 
radius  of  convergence  of  this  series  depends  on  the  frequency  co  of  the  wave  motion  and  is 
approximately  equal  to  (lyl/a)1/2  for  small  frequencies.  At  the  critical  high  frequency 
toCT  =  (3p(1)/Po3))1/2»  a  =  0  and  the  equation  (12.33)  can  be  written  as 

*  u(1)o<4)  i  3u(1)o(1) 

(■§■  r  *  T  ~  -  (  '  a,  -  i2  -  H(1))F(r)  -  0  (12.36) 

4  Po  3  Po0) 

Substituting  F(r)  =  £  anifl  in  this  equation,  the  coefficients  a„  can  be  calculated.  This  is  a  con- 

n=0 

vergent  seires  and  its  radius  of  convergence  is  given  below: 


4  W 
9  Hd)p(4) 


For  the  axial  shear  waves  if  we  substitute 


(12.37) 


uz(r,t) = G(r)eie“ 


(12.38) 
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in  the  equation  (12.23)  we  will  get  the  following  differential  equation  for  G(r) 


Pb®  dvpv 

G"(r)  + - - - - - - G'(r)  + - - - 2 - - - G(r)  =  0 

(M(1>  ”  J  Cl^Po^V  -  j  C02p^4)  0i(,)  -  y  -  j  G^p  <4) 

(12.39) 

The  origin  r  =  0  is  again  an  ordinary  point  for  this  differential  equation  and  its  solution  in  the 
neighborhood  of  r  =  0  can  be  written  in  the  form  of  G(r)  =  £  */•  At  the  critical  high  ffe- 

n=0 

quency  CDa  =  (3p(1  Vpi3^1^  the  equation  (12.39)  adopts  the  simpler  form 


4p<!> 

G"<r> - “ST  fG(r)  =  0 

Po() 


(12.40) 


This  is  the  Airy  differential  equation  and  its  solution  can  be  represented  in  terms  of  the  Airy 
functions  of  the  first  and  second  kind 


(12.41) 


Due  to  our  initial  assumption  that  the  thickness  of  the  micro-structure  is  very  small  and  also  by 
definition  (12.3)  we  conclude  that  the  coefficient  of  r  in  the  arguments  of  the  Airy  functions 


(12.41),  i 


is  a  large  number.  Therefore,  we  can  employ  the  following  asymptotic 


representations  for  the  Airy  functions 


Ai(ar)  =  yj--  (ar)"1/4  e  3 


Bi(ar)  =  (ar)_1M  e 3 


t(«)m 


(12.42) 
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For  the  radial  waves  we  have  the  more  involved  equation  (12.24)  and  if  we  let 


u^r.t)  =  H(r)ei0*  (12.43) 

we  will  obtain  the  following  differential  equation  for  H(r) 

+ T + T + T + T’,H"(r) + 7 <a" + T' + 7  +  7')H'(r) 

+  4<“Vil)r!-|  uV™  -(0,  +  —  +  3-  +  T-  +  T-»H(r)  =  0  (12.44) 

r2  3  t  t*  r  r 

where 

a0  =  X<,>+2ti(,)-y«>2pj3) 

ttj  =  y  (2X<3)  +  p<3>)  (12.45) 

03=4x<4>  +  p<4> 

4 

a4=i-(X(5)  +  2H(5)) 

h=-^>  .  fc  =  -0-{4)  +  H(4>) 

(12.46) 

P4  =  -  -j  (X(5)  +  2|i(5))  =  -2a4 


Yj  =  3X(2)  +  4|i(2) 

Y2=t(X(3>+2^(3>) 
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•ft-j  (X(4)  +  2h<4)> 

Y4=  (X(5)  +  2n(S>)  =  404 

Equation  (12.44)  can  be  written  in  the  following  form 

H"(r)  +  ~  p(r)H'(r)  -  -j  q(r)H(r)  =  0 

where 


(12.49) 


(12.48) 


p(r)  =  f(r)/h(r)  ,  q(r)  =  g(r)/h(r)  (12.49) 

arc  rational  functions  of  r  and 


f (r)  =  P4  +  £3r+pir2  +  a0r4 

g(r)  =  Y4  +  +  y?2  +  Yii3  +  (etc  +  j  ~  (12.50) 

h(r)  =  04  +  cy  +  c^r2  +  ajr3  +  ey4 

Unlike  the  two  previous  cases  the  origin  r  =  0  is  not  an  ordinary  point  for  the  differential  equa¬ 
tion  of  the  radial  waves,  but  it  is  instead  a  regular  singular  point  However,  since  the  indicial 
equation  at  r  =  0  is 


r2+(p(0)-l)r-q(0)  =  0  (12.51) 

and  by  (12.49),  (12.50),  (12.46)3  and  (12.47)4;  p(0)  =  -2  and  q(0)  =  4  we  always  have  an  analyti¬ 
cal  solution  at  r  =  0.  The  other  independent  solution  has  a  singularity  at  the  origin  and  can  be 
found  by  the  familiar  techniques  of  series  solution.  The  nature  of  the  general  solution  at  origin, 
as  far  as  the  analyticity  and  singularity  of  the  solution  is  concerned,  is  independent  of  the  fre¬ 
quency  of  the  radial  waves. 
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12.3  Governing  Equations  for  Spherical  Laminates 

In  this  section  we  use  the  results  of  section  9  to  derive  the  differential  equation  of  wave 
motions  with  polar  symmetry.  With  such  an  assumption  the  relevant  field  variables  depend  on  r 
and  t  only.  Consequently,  we  do  not  have  variations  in  01  =  $  and  02  =  0  directions  and  the  only 
non-vanishing  component  of  the  displacement  vector  u  is  u,  which  we  will  denote  from  now  cm 
by  u.  The  physical  components  of  the  relative  kinematic  measures  are  calculated  from  (9.13) 
and  (9.15)  and  the  results  are  recorded  below 


u  „  chi 
Y<#-Yee  =  —  ♦  Yn-  ^ 

Yoe  =  Y$r '  Yer =  0 


(12.52) 


1  d  .  . 

K^Kee^-j-^fru) 


(12.53) 


The  constitutive  relations  for  various  components  of  the  composite  stress  tensor  and  the  compo¬ 
site  stress  couple  are  derived  using  relations  (9.26M9.40)  together  with  definitions  (12.1)  and 
(12.2).  We  also  substitute  in  these  constitutive  relations  the  results  (12.52)  and  (12.53)  for  the 
relative  kinematic  measures.  After  some  simplifications  we  have 


tw  =  (2(X<>  V") + -  <X°  V2)) + -4  ft0  Vs) + -T  ftwv4,»- 

r  3r*  t  r 

+  {X(1>  +-(4  X*2  V2*)  +  4  (T  X(3 v3))  +  -T  (X(4V4))}  02.54) 

r  2  r  3  r  *» 


tM=  |2(X<>V")+ 7  p  (X'W  j  (X'4V4)»7 
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+ (x<‘» + 1  <4  x°v«>) + 4  <4  xm+*i<j>) + 4  <x<4v4>»  #  (12.55) 

r  2  i2  3  r*  or 


A.  x«>+ A.  ^>)u 

+  {X(,)  +  2n<15  +  i  (4  X<2>  +  2h0>)  +-4x0)l|r  (12-56) 

r  2  3  r4  or 

‘t^Q  “  —  tj,  ~  X|^  =  XfQ  =  0  (12.57) 

s„  =  (X«V2)  +  T-  ft(3V3>)  +  4  <x(4)+n<4)>  +  4  (X<5  V”))- 

3r  r*  5r  r 

+  { 1  x(2)  +  (5X<3>+2h(3))  +  4r  (7X(4>+6p<4>)  +  — — r-  (X(5V5)) )  4r  (12.58) 

2  3r  4r  5r  or 

see  =  (X®V»  +  4  (X°V3>)  +  4  (X<4)V4>)  +  4  (X<!V5)»- 

3r  r1  5r  r 

+  { 1  ^(2)  +  J_  (5X<3>+2h(3))  +  (7X(4)+6ji(4>)  +  — -r  (X(5V5)»  4r  (12.59) 

2  3r  4r  Sr3  or 


s^e-s^-s^-Ser-O  (12.60) 

The  equations  for  balance  of  linear  momentum  and  director  momentum  are  written  in  terms  of 
the  relevant  components  of  the  composite  stress  tensor  and  the  composite  stress  couple  by  using 
relations  (9.17)  and  (9.18).  In  the  absence  of  body  force  and  body  couple  and  recalling  the 
assumption  stated  at  the  beginning  of  section  12. 1  for  polar  symmetry  we  have 


cot  $ 
r 


(**<,- tee)  + 


1 

x^sin^ 


1 

l^sin^ 


(12.61) 

(12.62) 

BAS 
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(12.63) 

See)+rssii.904'0 

(12.64) 

r^sin*  °*,  =  0 

(12.65) 

(12.66) 

In  writing  down  these  equations  we  have  also  taken  notice  of  the  fact  that  the  only  non-vanishing 

component  of  the  director  vector  is  83  =  5r  =  .  Hie  equations  (9.19)  for  balance  of  moment 

or 

of  momentum  are  identically  satisfied  due  to  the  constitutive  relations  (12.57)  uid  (12.60).  We 
also  notice  from  (12.54),  (12.55),  (12.58)  and  (12.59)  that 

=  ^  =  ^e  (12.67) 

Using  these  results  together  with  the  equations  of  motion  (12.61),  (12.62),  (12.64)  and  (12.65), 
we  conclude  that  the  first  two  components  of  interlaminar  stress  vector  are  zero,  namely 

o*  =  ae  =  0  (12.68) 

So  we  have  only  one  equation  of  motion  which  should  be  derived  by  eliminating  ar  between 

equations  (12.63)  and  (12.66).  In  order  to  write  down  this  equation,  first  we  derive  appropriate 

expressions  for  the  composite  mass  density  p0  and  the  composite  mass  moments  p^1  and  p0^ 
by  using  relations  (9.41)-(9.43)  and  (12.3).  The  results  are  as  follows: 

p.=pJ1)+7Po<2) 

Poz‘  =  jPo0)+-|-Po0)  <12.69) 
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9or 

Next  we  differentiate  (12.66)  with  respect  to  r  and  substitute  for  in  (12.63)  while  making 


use  of  (12.69).  The  resulting  equation  of  modem  in  terms  of  the  composite  stress  tensor  and  the 
composite  stress  couple  components  are  as  follows: 


7  ^  (r2t"+rsee)  "7%©  =  (Po1}  “  -~2  Po0))u 

-  (|r  p.°> '  ♦ £p^  p n  0  <12TO> 

Now  we  substitute  the  constitutive  relations  (12.55),  (12.56)  and  (12.59)  in  the  equation  of 
motion  (12.70)  to  derive  the  displacement  equation.  After  simplification  we  get  the  following 
result 


(X<»+2,l<»  + 1  (7X°>+2h0>)  +  i  (7X»)+6p<4>)  +  p-  0-,5)  +  P,!>» 

+  (2(X<‘W>)  +  JL  X.P)  -  ±  (\P>  +  |  n<5>)  -  -L  (7A.w+6h<4>) 

-  {2(JL(1)+2p(1>)  +  -  (X(2)+p{2))  +  -K  0lX(11f28|i(,))  +  4  a(4V4)) 
r  3r  r 

+ -4  (*(5v5)»  4 

5r  r 

-<fc®-^P.^-<|:  *»+£*«>§ 

-(}poO,+  ^:Pi4,)0  02-71) 
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Again,  as  in  the  case  of  the  cylindrical  laminates,  we  expect  to  recover  the  classical  theory  of 
wave  motions  with  polar  symmetry  through  a  limiting  procedure  in  which  the  thickness  of  the 
micro-structure,  approaches  zero.  Recalling  (12.1)-(12.3)  and  suppressing  the  superscript  (1) 
at  the  limit,  the  equation  (12.71)  reduces  to  the  following  form: 


9^  2  9u  2u  _  1  9*u 
dr2  r  9r  "  r2  =  9t* 


(12.72) 


where 


c?= 

Po 


(12.73) 


This  is  the  familar  displacement  equation  of  motion  for  waves  with  polar  symmetry.  In  order  to 
find  the  general  solution  of  this  equation  it  is  convenient  to  express  the  radial  displacement  u(r,t) 
in  terms  of  a  potential  function  $(r,t)  through  the  relation 


»  =  •& 
or 

Now  the  second  derivaqive  of  the  product  rtj)  with  respect  to  r  is 


(12.74) 


9 2  ...  92<fr  _  9<b  9u  ,  « 

If  we  rearrange  the  above  result  and  differentiate  again  with  respect  to  r,  we  get 


9  .  9u  2  .  92u  2  9u 

37(37  +  7“)=3?  +  7  37 


2u  3  .1  S2  .... 


which  reduces  by  (12.72)  to 


3  . 1  3J  ,  _  1  3*u  _  1  3  ,3V 

3r‘r  ar2  W)1'<g  St2  ~  Cl  *(^) 

Therefore,  we  conclude  that  if  the  product  nj>  satisfies  the  following  one-dimensional  wave 


BASE 
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equation 

il  (4)  =  JL  il  (r*)  (12.75) 

dr2  Cl  d t2 

then  u(r,t)  as  defined  by  (12.74)  will  satisfy  the  displacement  equation  (12.72). 

The  general  solution  of  (12.75)  is 

<*r,t)  =  |  f(t-r/CL)  +  ~  g(t+r/CL)  (12.76) 

where  all  the  two  terms  represent  waves  diverging  from  the  origin  r  =  0  and  converging  to  r  *  0, 
respectively. 

The  displacement  equation  (12.71)  for  waves  with  polar  symmetry  in  spherical  laminates 
can  be  investigated  by  substituting  a  separated  solution  of  the  form 

u(r,t)  =  F(r)ei‘“  (12.77) 

The  differential  equation  for  the  spatial  part  of  the  solution,  F(r),  can  be  written  in  the  following 
form 

F"(r)  + 1  P(r)F(r)  +  ~  Q(r)F(r)  =  0  (12.78) 

r  i2 

where  the  rational  functions  P(r)  and  Q(r)  are  given  by 

P(r)  =  N1(r)/D(r)  ,  Q(r)  =  N2(r)/D(r)  (12.79) 

and  Njfr),  N2(r)  and  D(r)  are  the  following  polynomials 

Nj(r)  =  2ct0r4  +  pji3  +  for2  -  a3r  -  a4 

N2(r)  =  a^p^V5  -  a',/  -  40^  -  ‘fer2  -  y3r  -  3ot4  (12.80) 
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D(r) *=  a,/*  +  ctjr3  +  a^r2*  <*$+  04 

where 

aD=X<1)+V1>-|o>2p0(3) 
a'0  =  2(0^  +  y  (&2p£h 
aj  =  2  (X(2)  +  -  ~  co 2p«<4) 

a2  =  l(7X(3)  +  2p(3)) 
ct3  =  j(7X(4)  +  6p(4)) 

(12.81) 

04=|^  +  ^5>) 

Pi-jX®-7«W 

p2  =  -2(X(3>+|p<3>) 

Y2  =  j(3a(3)  +  28p<3)) 

Yj  =  6(X(4)  +  p(4)) 

The  origin  r  =  0  is  a  regular  singular  point  for  the  equation  (12.78).  The  indicial  equation  at  this 
point  irrespective  of  the  material  constants  is 

r2  -  2r  -  3  =  0  (12.82) 

It  is  obvious  from  (12.82)  that  we  have  always  an  analytical  solution  at  r  =  0  and  the  other 
independent  solution  has  a  simple  pole  at  this  point.  By  (12.81),  the  radius  of  convergence  of 
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